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ABSTRACT 


The problem of allocating irregular shapes on a plane sheet to 
minimise the required area has been studied extensively in 
literature. The important parameter affecting the nesting problem, 
that of the shapes of the pieces and the allocation area has not 
received sufficient attention. The present work is concenti-ated on 
the study of the range of two dimensional shapes from l.he view of 
matching the shapes for efficient nesting. 

The aim of the work is to propose a comprehensive theory for a 
descriptive analysis of shapes, so as to facilitate shape 
matching. The shape description problem has been dissected in to 
subproblems of polygon decomposition, a classification of the 
decomposed primitives, the topological and geometrical inter- 
relations a.nong the primitives, and recognition of shape features 
for nesting. 

A complete polygon decomposition algorithm has been proposed which 
is aimed at giving a decomposition of the profile in to a set of 
convex polygons approximating the original profile; the objective 
of the decomposition is to provide a set of constituent polygons 
which would be useful for nesting oriented shape description. 

A complete convex polygon classification and description has been 
proposed in the present work. It is based on the number of sides 
and the shape f eatures such as prominent and sharp angles . 

The topological and geonmt.rical relations between the constituent 
convex polygons have been quantified, so as to provide simple 
nesting f eature r ecogn i 1; ion criteria. With this, the shape 





description is deemed to be complete. 

The proposed l.heory and algorithms rely heavily on artificial 
intelligence and fuzzy theoretic techniques. 

The present work is deemed to be a fundamental study of shapes, in 
addition to being a nesting problem solution. The concepts evolved 
in the thesis would be useful in the solution of allied 
geometrical problems such as the automatic machining sequencing 
problem in manufacturing. A conceptual level algorithm has been 
given in the appendix to illustrate the applicability the concepts 
to the allied problems. 



Chapter 1 


INTRODUCTION 

1 . Notivat ion of the Theaia 

Geometry problems are fundamental problems in mechanical 
engineering- They may appear as the inain problem in mechanical 
engineering research; or as subproblems in course of some other 
research problem. They are encountered in almost all fields of 
mechanical engineering. Some of the more prominent fields in which 
geometry problems are crucial, are listed below. 

Design and CAD : Any machine element design is finally a geometry 
sped f icat ion for manufacturing. The geometry aspect of many of 
machine design problems has often been relegated to the level of 
det erminat ion of size parametervS of some standard geometrical 
shapes. A full research on the possibi 1 i t i es offered by the 
'shapes' of the bodievS has rarely been done. Shape optimality of 
design is a field which would definitely be benefited by basic 
geometry investigations. 

Solid mechanics : A solid mechanics problem of stress and strain 
analysis of a body has geometry problems in the guise of the 
geometry specification of the body to which various force fields 
are appl i ed . 

Finite elements method : Apart from the solid mechanics 
aplications of FEK, the other problems such as heat transfer 
through solids, have geometry aspects in the form of the boundary 
constraints for the differential equations. Also, the automatic 
mesh generation for FEN, is a mainly geometry problem, 
nanufacturing process planning : Cutting processes can be viewed 
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as purely geometry problems > v/here in a given geometry of the 
stock is to be reduced the final geometry of the machined part. 

Any theoretical invest igat ion in to georaetry, clarifying the bavSic 
concepts, is likely to find applications in many of these fields. 

2 . The Meeting Probl am 

The nesting problera, in its purely geometrical form, finds 
applications in optimum sheet layout for sheet metal processing; 
and in similar situations such as garraent industry, leather 
industry. 

Stat ement of th© Problem 

The problem tackled in the present work is th© two dimensional 
nesting problem. 

A given sot of two dimensional profiles along with the number of 
pieces of each of the profiles, and the specifications of two 
dimensional rectangular sheet, the profiles are to be allocated on 
the sheet, without any overlap between them, so as to minimise 
the area of the sheet utilised for al locat icn.^ . The nesting 
efficiency is a measure of the area utilisation of the sheet. It 
may be defined as the ratio of the area of the profiles allocated 
to the utilised area of the sheet. Hence the objective function of 
the optimisation is the maximisation of the nest ing ef f ici ency . 

The profiles considered in the present work are non- 
self intersecting polygons, whose boundary is composed of straight 
line segments, without holes. 

3 . Basic Approach to the Meat ing Problem : 

The nesting problem can be viewed as an optimisation problem. The 
objective of the problem is to determine an allocation scheme to 
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minimise the cirea of the sheet. The constraints are the non- 
overlap of the allocated profiles with each other and with the 
exterior of the sheet. 

Due to the extreme diversity of the shapes of the profiles, a 
formulation for closed form .solutions seems irnpossible. Even 
application of classical search methods in optimisation would face 
very hi£h difficulties at the formulation stage it.self. Hence the 
possibilities of meaningful formulation of the problem as an 
optimisation problem and a successful application of a classical 
opt iraisat ion method seem bleak. 

The nesting problem is solved in industry with little difficulty, 
by human experts. The solutions reached may not be optimal; but 
they are sufficiently good for practical purposes. An attempt at 
emulation of the*, qualitative reasoning carried out by human beings 
would be the best approach to tackle the problem. 

The reasoning of a human expert, when solving a nesting problem is 
based on the ’'shape’ and ’size’ qualities for identification of 
the profile best suited for an allocation in the complimentary 
area of the sheet, and a .shape and size match for the eventual 
allocation. 

The reasoning of a human expert is essentially vague and inexact 
in nature. Fuzzy theory is a very good theory available for the 
modelling of such vague and inexact processes. 

A human expert bases his decisions on the experience gathered from 
the past attempts at nesting. The solution process is a search 
method with tentative decisions, backtracking. The decisions taken 
are based on heuristics. 

The above discussion suggests an approach based on shape analyais 
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and matching. The algorithins for shape analysis and matching would 
be fuzzy algorithms. The nature of the algorithms would foe of the 
type which may be classified as the heuristically guided search 
methods . 

In the present work, basic element.s of artificial intelligence and 
fuzzy tlieory have been tised to model the qualitative, inexact and 
the tentative aspects of the problem. The main thrust of the work 
is on shape analysis. A fundamental study of shapes from geometric 
perspective has been attempted in this thesis. 

4 . Discussion on the Problem and Solution Approach 2 . 

The basic idea behind the nesting, allocation algorithm is 
envisaged as follows. 

After allocation of some profiles, the utili,sed region of the 

sheet takes the form shown in figure 4.1 (a). The shaded area, is 

the complimentary area of the sheet, on. which the next allocation 

is to take place. The complimentary area has some shape features, 

projection P and bays B and B shown in the figure. The set of 
1 12 


Complimentary 

Area 




Allocated area 


(a) (b) 

Fig. 4.1 Basic Idea of Nesting. 



(c) 
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profiles yet to be allocated may contain profiles which may have 

similar shape features, A profile P , shown in figure 4.1Cb), has 

f 

a projection and a bay which may approximat ely match the 

projection P and bay B of the complimentary area. This profile, 
1 ' 1 

if placed in the complimentary area as shown in figure 4.1Cc), 
would maximise nesting efficiency by minimising the wasted area. 
This approach entails an invest igat ion of the following 
subprobl eras . 

1) A profile description based on shape features such as 
projections and bays and their int er-relat ions . The profile 
description would be an enumeration of the features and their 
inter-relations . 

2) Identification and devScription of the projection and bay 
features. 

A profile may have extremely complex shape rendering a shape 
analysis for identification of features and their inter-relations 
correspondingly di f f icult . The difficulty is sought to overcome by 
decomposing the profile in to a set of convex polygons which are 
simpler and hence more amenable to analysis. The topological and 
geometrical relations among the constituent convex polygons are 
then analysed. Projection identification and description is done 
on the basis of the information available at this stage. Bays are 
identified by taking a convex hull of the profile. Bays are the 
differential areas between the convex hull and the profile- Bay 
analysis is done by treating the bays as independent profiles. The 
bays, projections and the geometrical and topological relat ions 
between them form the prof i 1 e descript ion . 


5 



Thus the following subprobleras have been identified. 

1) Profile decomposition in to a set of convex polygons. 

2) Convex polygon description. 

3) Formulation of geometrical and topological inter-relations 
among a set of convex polygons and their conU^inat ions . 

4) Projection and bay identification and' analysis. A shape 
description as the projection and bay features and the remaining 
constituent convex polygons, and their inter-relations. 

5) Shape matching based on the profile description formulated in 
the steps given above. 

6) Allocation of a profile based on the shape matching, subject to 
the constraints of non-overlap of the profiles. 

The shape analysis and matching are the crux of the approach. The 
basic thrust of the present work is in shape analysis. The aim of 
the shape analysis is to give an abstract, qualitative shape 
description which is none the less complete. To give a rough 
analogy, the shape description gives a blurred picture. Some 
initial tentative decisions are taken on the basis of the blurred 
picture. As the solution proceeds, the picture is brought in to a 
sharper and sharper focus. The description is interpreted more and 
more exactly and quantitatively, as the decisions become more and 
more firm. 

The description of a shape shoxild necessarily have the following 
properties. The shape description should be unique to the profile 
being described. It should provide a one to one mapping from the 
space of all the shapes on to the space of descriptions. To 
preserve the abstraction and qualitative nature of the description 
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and at the same time sati.sfy the above criterion, the complete 
description should be composed of a number of sub-descriptions 
with varyin£ degrees of, what could be termed as 'closeness to the 
profile’ or 'goodness of description’. None of these sub- 
descriptions, by themselves, can provide one to one mapping. The 
conjunction of these descriptions with their individual 'goodness’ 
gives the description of the profile. It should be noted that the 
solution methodology should, hence, produce multiple shape 
descriptions . 

The shape matching is the matching of the descriptions. As opposed 
to a problem in machine vision, the shape matching here, is 
partial shape matching rather than full shape matching. A local 
region of the profile is to be matched against the complimentary 
space on the nesting sheet. Hence the local boundary features such 
as projections and bays .should form the basis of the shape 
descriptions. 

The shape analysis and matching should necessarily be general in 
nature. It should cover the whole range of possible two 
dimensional range. This is necessitated by the fact that though 
the shapes of profiles may be constrained in the problem 
definition. But the shapes generated in the complimentary area of 
the sheet during allocation of the profiles are beyond control. 
The complimentary area may assume any arbitrary shape even for the 
simplest of profiles. Hence any proposal or a study of nesting 
algorithms must be preceded by a complete study and analysis of 
arbitrary two dimensional shapes. 


7 



5 . Heuriat ic3 in the Solution Hethodology j_ 

The advocated solution process is the process of attainment of a 
goal state from an initial state. The initial state is the 
specification of the sheet and the profiles. The goal state is the 
final optimal allocation of the profiles on the sheet. The goal 
state may not be unique. 

The solution process is subdivided in to the subproblems stated in 
the previous section. Each of the subprobleins may be treated as an 
independent problem, with its initial state and goal state. For 
-many of the subproblems, the goal state evaluation is heuristic in 
nature. An attainment of multiple goal stat'es is desirable, 
especially for shape descriptions and shape matching. Hence the 
algorithm can not terminate on reaching a goal state. Multiple 
paths need to be explored in the solution tree, to obtain multiple 
goal states. 

The traversal from the initial state to the goal state may be 
modelled as a tree traversal, or a tree search problem for an 
optimal goal state. At every node of the tree, a number of 
solutions are generated, which are the children of the current 
node. The generation of the children at each of the nodes is done 
on the basis of the conceptual framework developed in the thesis. 
The solution tree may be very large. Hence an exhaustive tree 
search is not advisable. The children nodes are evaluated using 
■heuristic functions developed for the subproblem. The best of the 
children is explored first. On reaching a dead end or a goal 
state, the solution path is backtracked and multiple goal states 
are generated. The extent of backtracking i.e. up to which of the 
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nodes in the existing solution path should the backtracked, and 
the number of goal states, are determined again by heuristics. 

The Solution method proposed, thus, may be called a heurist i cal 1 y 
controlled generative space search. The heuristic functions are 
applied for evaluation of the solutions and goal state.s, and for 
controlling the backtracking [28]. 

6 ■ The Fuzzy Theory Application in the solution ± 

The fuzzy theory is a good tool to model ’vague’, ’qualitative’, 
’inexact’ reasoning. A good development of fuzzy theory has been 
given in standard textbooks on fuzzy theory [13, 14, 19 & 20]. In 
the present work, basic elements of fuzzy theory have been used to 
handle quantifications of the qualitative aspects of the shape 
analysis . 

If the solution proces.s is modeled as a decision tree, then at 
every decision node a number of decisions are possible. The 
heuristic functions evaluate the possible decisions i.e. the 
children of the current node. The quantity returned by the 
heuristic functions is taken to be the ’desirability’ or 
’goodness’’ of the children. This is illustrated in the Fig. 6.1. 



Fig. 6.1 Heuristic Decision Tree or Fuzzy Spawning Process The 
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node N has children N , W , N and W . Their goodness measures 
i j k 1 m 

are given in the figure as 0.1, 0.4, 0.8 and 0.7, respectively. 

The node N and W have a low goodness. Hence they should be 
j k 

pruned. Nodes N and N have comparable and high goodness. They 

1 ffl 

are both expanded, to explore the solution space. At each decision 
node more than one decisions may be followed, thus spawning a 
number of subprocesses. This type of proce-sses have been termed as 
fuzzy spawning processes in the present work. 

The goodne.ss of the achieved goal may be evaluated independently 
by a heuristic function or it may be evaluated as a function of 
the goodness of the decision nodes in the path. This function, 
called the mapping function in literature, has to satisfy 
necessary constraints [19]. Sufficient conditions for the mapping 
functions have not been formulated. Some common mapping functions 
have been given in literatvire [13, 14, 19 & 20]. Both the-se 

approaches, for evalviation of goodness of the goal, have been 
advocated for different subproblems in the present work. 

The heuristic functions which are the same as the fuzzy membership 
functions, are based on some qualities in the problem domain. They 
quantify the qualitative aspects. They are also termed fuzzy 
qualitative predicates as they detect the presence of a quality 
and return a number indicating the measure of the quality in the 
solution or decision alternative. On detection of a quality in the 
domain, a corresponding action is recommended by a rule. This 
action gives rise to a decision alternative. 

The solution process may be viewed as a forward chaining of the 
fuzzy qualitative predicates, with .fuzzy inf erencing and with 
backtracking which is controlled by some fuzzy qualitative 
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predicates. The fuzzy spawning model fits in to this paradigm as a 
backtracking done to generate multiple goal.s. This paradigm is 
more convenient for implementation. The research done on fuzzy or 
possibilistic inferencing [14] will be useful for a thorough study 
of the problems through experimentation. The generalised modus 
ponens propounded in [14] is suggested as a fuzzy inferencing 
likely to give good results. 

The general method followed for quantification of the qualitative 
aspects for the formulation of the fuzzy qualitative predicates is 
as follows. The rough quantity indicating a measure of the quality 
is taken, e.g. for largeness of geometrical .size, a length or a 
radius may be chosen as the quantity. It is compared with a 
particular quantity of the sa.me type characteristic to the profile 
under consideration, e.g. for a length, it may be longest distance 
between two boundary point.s of a profile. The normalised parameter 
obtained after comparison, is subjected to a fuzzy membership 
distribution for an evaluation of the qualitative predicate.- 
Consider the predicate ’Long Length’ for an example. The length to 
be evaluated under this predicate is normalised by taking a ratio 
with a characteristic length. Let the length ratio be L . The 

V 

fuzzy membership distribution for ’Long Length’ is shown in Fig. 
6.2. The membership value is the value corresponding to the value 

of L . The fuzzy membership distribution functions embody the 

r ■ . • ' 

heuristics; they need to be refined by experimentation over a 
range of cases. For the beginning, a linear distribution of the 
type shown in Fig. 6.2 is proposed. This process is called fuzzy 
normalisation. In the thesis, a generic function ’ f uzzy_normal ise ’ 
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Mefnb ei'-ahip 


Length ratio > L.^ 

Fig. 6.2 A Fuzzy Membership Function. 

has been used to indicate this process. The reference quantity, 
the characteristic length in the above example, is called the 
normalisation and comparison parameter. 

7 . Scope of Work ± 

The present work is limited to a development of the basic concepts 
required for a complete theoretical system for shape analysis 
AlgorithmvS have been proposed; but they are more in nature of 
first proposals demonstrating the important niche of the concepts 
which have been developed. They should be treated as initial 
proposals providing fundamental guidelines for further work in the' 
field- The algorithms have not been studied for complexity and 
implementation aspects. The attempt here has been to reduce the 
problems to the level of standard problems in the fields, such as 
path traver-sal in a graph, for which standard solutions are known 
to be available. Computational complexity has not been a 
consideration in any of the problems. Stress has been laid on the 
Gompleteness of the theory, and correctness of approach, which 
appears to be missing in the research work in the field of shape 
analysis 
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The present work is meant to be an exploration of a fundamentally 
new approach to geometry and shape handling in mechanical 
engineering problems. The concepts developed in the present work 
have a powerful potential for giving complete solutions of 
hitherto unsolved geometrical problems in mechanical engineering. 
As an example of the applicability of this approach, an algorithm, 
at a conceptual level, has been given for the solution of a 
process planning problem in the appendix. 

The statement of problem has been given in section 2 . A non- 
self intersecting polygon, constrained by linear boundary segments, 
may have very short line segments in the boundary approximating a 
curved boundary. The curved boundary segments are identified and 
treated as curved rather than as being composed of line segments. 
Hence an arbitrary two dimensional shape can also be analysed by 
the proposed method. 

The basic thrust of the present work is on shape analysis. The 
underlying philosophy of the approach has been to obtain an 
analysis consistent with human perception. The subproblem 
definitions have been proposed on thi.s basis. 

Shape matching for nesting should be based on local feature match 
such as the projections, bays and strings composed of projections 
and bays along the boundary. It is a partial and local shape match 
as against a complete and global shape match required in machine 
vision and part recognition sy,steras. A complete and global shape 
match is not much relevant in nesting. A partial and local match 
is encountered more often; and hence description of the local 
boundary features and their inter-relations should form the basis 
of the shape description. 


13 



A large nuBiber of techniques may be devised for nesting and 
allocations. An explorative study of nestixig problem and the 
possible algorithms would be very vast. It has not been attempted 
in the present work. A basic approach of shape matching, which is 
likely to be extreniely fruitful, has been pursued. 

The shape analysis is fuzzy in nature. The solution niethodologi 
applied, are based on artificial intelligence techniques. The 
special problems of fuzzy algorithmvS and fuzzy inferencing have 
not been studied in depth. It has been assumed' that mechanisms 
suitable for the proposed algorithms, may be vselected from thosc-^ 
available in literature. Some suggestions for the mechanisms have 
been given from time to time. 

The artificial intelligence aspects of the problem are also not 
dominant. A very basic technique of heuristically controlled 
search has been applied. In this aspect also, an in depth study of 
heuristically controlled search has not been done. The proposed 
heuristics need to be tested and refined through wide ranging 
experimentation. 

8 . Comparison of the Approach in Present Uork with some Standard 
Approaches 2. 

The skeleton of a two dimensional profile [8] has been selected to 
be the basis of the whole shape analysis. The reasons underlying 
this selection have been laid out in section 1 of chapter 2. 
Diverse approaches have been tried for shape analysis by various 
researchers. The details of these approaches are given in the 
literature review. The approach followed in the present work Is an 
’area analysis’ approach, with profile decomposition and the 
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relations between the primitives formed by the 
decomposition, rather than a boundary analysis approach. A detailed 
di.gcussion of this issue is given in the chapter 2. 

Another rough classification of the approaches for shape analysi.s 
is between ’syntactic pattern recognition’ and ’graph theoretic’ 
approaches. The approach followed in the present work could be 
classified to be under the graph theoretic approaches. 

One of the major problems with the syntactic approach is that it 
is, if not impossible, extremely difficult to develop 
computational grammars, which are applicable to the whole range of 
shapes. This problem has been mentioned in [18], The most basic 
di f f icul t i ej3 in this approach are the problems of discretisation 
of a continuum. For example, in a syntactic boundary analysis, the 
alphabet for the grammar is composed of short, directed vectors, 
of unit length. For a profile, in which the lengths of the linear 
boundary segments have a large range, say varying from one to 
fifty, the parse tree becomes extremely large. Attributed entities 
need to be used to overcome this problem, e.g. for boundary 
analysis, the length of the vector is the attribute of the vector. 
The next, and even more fundamental problem is that, what should 
be the orientation of the vectors. Obviously, unit vectors in all 
the directions can not included in the alphabet; the alphabet size 
would become infinite. The approach followed so far in the 
literature is that, either consider only polygons composed of 
boundary segments which are oriented in certain fixed directions 
only [17], or consider the boundary segments oriented in some 
range, say 0 < 0 <90, to be exactly the same [27]. None of these 
approaches seem satisfactory. 
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9 , Lit erat ur e Survey ^ 

The literature review consists of review of the nesting algorithms 
for irregular shapes and the review of the subprobleras and their 
solutions . 

9,1. Nesting ^ 

An exact formulation of nesting problem has been given by Bailleul 
et al [6]. It can be seen that the problem can not be solved 
analytically to get a closed form solution. 

A nesting problem with just one shape fallvS under the category of 
lattice formation. An approach to nesting .solution which 
constrains the solutions to single row, double row or in general a 
layout strictly forming strips has been explored in a lot of 
papers. An approach where in the profilevS are enveloped in 
rectangles and the rectangles are allocated on the sheet has also 
been quite popular due to its relative simplicity. Such 
constrained nesting layouts, in general, would not give good 
nesting efficiency in all cases. Since the approach followed ih 
the present work is extremely different from these approaches, 
they are not reviewed here. A good review of such approaches is 
given in [ 2 2 ] . 

Albano and Adamowicz [3] proposed an algorithm which clusters 
irregular shapes in to rectangular modulevS of minimum area, bavSed 
on some heuristics, such as the best pairwise clustering is 
obtained by placing two pieces of the same profile next to each 
other along the width, aligning the lengths, and rotated through 
180 degrees. The proposed rectangular module formation is an 
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exteiivsion of the pairwise clustering hei 4 ristic al£orithro. Albano 
and Sappupo [4] gave a heuristic algorithm to allocate these 
rectangular modules on a sheet based on an OvSt imation of the waste 
that will be generated by the future al locat ioxis . Kaheshwari 
combined parts from both these algorithiftvS and modified the 
heuristics to give a nesting algorithm which short listed profiles 
based on equality, of lengtiis, giving priority to longer profiles 
and laying the short listed profiles from bottom left upwards 
[22], Albano reported an interactive implementation of the 
algorithms in [3] [4], which allowed the user to control and 
modify the allocations proposed by these algorithms [2]. Bailleul 
and others followed more or levSs the same approach, i.e. creating 
strips by clustering profiles and allocating the strips on the 
sheet [6]. They have implemented an interactive program which 
presents allocation schemes of the user and allows him to select 
the one which is consistent with the anisotropic properties of the 
sheet and the profiles. 

These nesting algorithms do not take in to consideration the basic 
property that dominates the nesting problem, that is the shape of 
the profiles and the shape of the complimentary area of the sheet. 
The first such attempt was done by Cai Yuzu et al [10]. They 
implemented an expert system for the shapes encountered in the 
garment industry, matching the projections and the bays of the 
pieces for allocation. The expert system does not analyse shapes 
but the shape knowledge, in the form of bays and projections, and 
allocation heuristics knowledge for previously known types of 
profiles is available, is encoded in the rules in the knowledge 
base. The expert fxinction executes on the standard paradigm of 
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rule application and limited backtracking. 

A nesting solution methodology should necessarily be based on 
shape analysis, as the property of ’’shape’’ forms the basis of the 
nesting problem. For a nesting solution to be considered a fully 
automated solution, it should be able to analyse and allocate 
diverse shapes and not a limited subset. 

9.2. Shape Analysis Techniques 

Various techniques have been employed for shape analysivS, such as 
thinning algorithms, statistical feature recognition and analyvSis 
techniques, Fourier Transforms of shapes, Syntactic analysis of 
the boundary and area, shape decompositions as area decompositions 
and boundary decompositions, and shape descriptions. The basic 
pattern recognitions methods applied in shape analysis are given 
in texts [16 & 25]. The special techniques evolved for shape 
analysis have been surveyed by Pavlividis in [24], and by Shapiro 
in [31]. The shape analysis work done in pattern recognition is 
mostly aimed at machine vision and automatic part recognition 
systems. Hence most of the research work ivS done on the assumption 
of closed world, i.e. a fixed set of well known profiles, or a 

tf 

class of profiles are analysed and matched for identification. 
Such work, where it is felt that it would be difficult to extend 
the techniques to diverse shapes, has not been included in the 
literature review. 

Skeleton, alternatively called Hedial Axis Transform (HAT), 
Symmetric Axis Transform (SAT) in the literature, is one of the 
basic techniques which have been quite popular with researchers in 
shape analysis. It is reviewed in the next section. 
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9.3. Skel eton j_ 

The skeleton was proposed by Blum as a transform of two 
dimensional .shapes [8]. Since then it has attracted a lot of 
attention from, researchers in the shape analysis field. The 
skeleton and some of its basic properties can be found in .standard 
texts on pattern recognition [15, 16 & 25]. Blum and Nagel 
investigated its properties and the possibility of its 
applications in shape analy.sis [9]. Some of its features .such as 
the forking points, skeletal interaction zones ha.ve been mooted by 
Seng - Beng Ho and Charles Dyer [29]. ria.ssone and I1oras.so have 
implemented a two dimensional shape editor based on a skeletal 
representation of shapes structuring a symbolic description of 
shapes and a shape topology [23]. 

As most of the work in shape analysis is aimed at machine vision, 
most of the algorithms for computation of skeleton work on 
digitised images. Digitisation, though it allows construction of 
simple algorithms for skeleton computation, produces its own 
problems, the foremost being that the skeleton is computed in the 
form of a set of pixels. A recovery of a line diagram description 
of the skeleton from the set of pixels is difficult and 
computationally complex. D. T. Lee has constructed a fast 
algorithm foi' ‘Computation of skeletons from direct boundary 
representation [.21]. 

As the skeleton gives a good structurisation to multif erously 
branched shapes which is consistent with an intuitive concept, it 
has often been used to decompose profiles using the lines forming 
the skeleton, in to basic shapes and projections. Arcelli and di 
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Baja applied succevssive threv^holding on the skeletal radiUvS avS the 
method of segmentation [5]. Pigier et al also used vsuccessive 
resolution reduction on object boundary resolution, by blurring 
the intensity extr-erna on a gray digitised image [26]. At every 
stage of the resolution, the skeleton is computed and a hierarchy 
is imposed on the lines forming the skeleton. 

9.4. Polygon Decomposition 

Some of the polygon decomposition methods which directly use the 
skeleton have been mentioned in the above section. They are 
different from the standard polygon decomposition problem in the 
sense that they do not explicitly decompose the polygon in to a 
set of convex polygons. Even a convex polygon may be decomposed 
where as a resultant of the decomposition may not be a convex 
polygon. Hence these methods are not included in this section. 

The polygon decomposition methods may be divided in to area 
decomposition methods and boundary decomposition methods. The 
first area decomposition method reported in literature is that of 
Pavlividis [25]. It follows the boundary of the polygon and' 
eliminates concave vertices sequentially. The concave vertex 
elimination is restricted to cuts joining the concave vertex with 
some other vertex of the polygon. 

Ahuja [1] has presented a decomposition method based on the Vornoi 
diagrams of the vertices of the profile. The concave vertices are 
eliminated by drawing a line (or lines) to the nearest Vornoi 
neighbour. Hence this algorithm can decompose the profile only in 
to convex polygons which are formed by Joining the vertices. 
Chazelle [11] has formulated an algorithm which again follows the 
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boundary of the profile and eliminates the concave vertices. The 
concave vertices are eliminated by cutting the profile by a line 
lying within the reflex angle of the convex vertex, and joining 
the concave vertex with another vertex of the profile. The 
possibility of intersection of cuts has been considered. The 
intersecting cuts are modified so as to ntinimise the number of 
constituent convex polygon. The cuts may also be between a concave 
vertex and some points on the boundary, which are not the 
vertices. The results are very good for most of the examples given 
in the thesis. But as in all other methods, the factor of the 
id ent f i abi 1 i ty of the constituent polygons ha,s not been 
considered. Hence, though the decomposition may be optimal from 
the minimality of the number of convex polygons, it may not give 
results consistent with human perception. 

The boundary decomposition method was first propounded by Shapiro 
and Haralick [32]. They have developed a graph theoretic 
clustering algorithm based on the mutual visibility of the 
vertices across the interior of the polygon, to decompose the 
boundary and then form the constituent polygons form the 
combinations of the decomposed boundary. This algorithm does not 
guarantee a decomposition in to convex polygons [24], 

Bjorklund and Pavlividis [7] extended upon the boundary 
decomposition method by using a graph model Tor relations among 
the boundary elements in addition to the mutual visibility. The 
links of the graph, termed the similarity graph, are the adjacency 
and similarity relations between boundary elements. The complexity 
of the graph is curtailed by considering only the links with the k 
closest neighbours and deleting the rest of the links. The set of 
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matching algorithm, for the relation graph matching. A structural 
description formalism has been proposed, which is used in the 
shape matching algorithm. Both, the description formalism and the 
shape matching algorithm give very Interesting leads which may be 
useful in a complete theoretical formulation. 

9.7. General Theoretical Background _j_ 

The application of the fuzzy theory and artificial intelligence 
methodologies is of a basic level in the present work. Hence a 
comprehensive literature review of these fields has not been 
included in the thesis. Basics of the fuzzy theory, used in the 
the present work, may be obtained from some standard texts [13, 
14, 19 & 20]. The fuzzy classification of convex polygons is 
similar to the fuzzy clustering given in [20]. The inexact 
neighbourhood relations developed in the chapter 5, are similar to 
fuzzy graphs [13]. Possibilistic or inexact reasoning has been 
developed in [14]. The rudiments of heuristically controlled 
search methods may be referred from [28]. 

Graph theoretic concepts have been used liberally throughout the 
present work. These concepts can be referred from a standard text 
book on graph theory [12]. 

10 . Organisation of the Thesis Report and Conventions ; 

10.1 Organisation x_ 

The chapter one is the introductory chapter. The basic approach, 
literature review of similar problems are given in it. In the 
second chapter, the skeleton of a two dimensional shape, which 
forms the basis of the present work is elaborated. An algorithm 
for profile decomposition is proposed in the third chapter. An 
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analysis leading to a complete classification and description of 
general convex polygons is given in the chapter four. The 
relations between the constituent convex polygons of a profile are 
developed in the fifth chapter. Chapter six concludes the thesis 
by discussing the present work, stating the scope for further work 
on this problem. A solution of a process planning problem has been 
given in the appendix. 

The sections in the chapters are numbered, starting from one, in 
each chapter. This is done as every chapter is essentially 
independent; also, to avoid long numbering, as the hierarchy of 
numbering is too deep. 

10.2 Diagrams 

Since the present work is geometrical, it has been illustrated by 

ample diagrams.. In the diagrams, very often, the skeleton of the 

profile is drawn along with the profile. The profile edges are 

labeled by the vertices V . In some diagrams, the profile boundary 

i 

and its interior are distinguished from other lines by short 
hatching. The skeleton is drawn by a dot-dash line (centre line). 
The cuts and approximations are denoted by a section line. In 
chapter 5, the N-relations are indicated by a double arrowed line 
segment. 

The diagrams are essentially illustrations of the principles 
explained in the respective sections. To maintain the link between 
chapter sections and diagrams, and to facilitate cross 
referencing, the diagram numbering ia not contiguous. It follows 
the number of the section in which the figure appears. 
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Chapt er 


SKELETON 

1 ■ Requirementa of a Basis f or Shape Analysis ^ 

Shapes of polygons can be extreraely diverse. A precise analysivS of 
shapes is a highly complex l:ask. A feature based cluster analysivS 
approach may be possible for a limited number of shapes. But when 
it comes to an analysis of all possible shapes, such an approach 
would produce a very high dimensionality; and also the generation 
of rigid and mutually exclusive partitions in the feature space 
denoting classes of shapes, as required in a crisp or non-fuzzy 
cluster analysis, would be impossible- Rather than a pattern of 
clusters of points, a continuous vspectrum would be obtained. The 
selection of features, for such an analysis is extremely 
difficult. An approximate analysis using fuzzy partitioning in 
some form, is necessary. 

Even before an approximate analysis is attempted, the shapes need 
to be transformed into a form more amenable to analysis, reducing 
the dimensionality and complexity- The transformation should 
necessarily ensure that there is no loss of information . 

The features which form the basis of the shape • analysis should be 
global rather than local. Features based on boundary analyvSis are 
likely to be local features influencing a small local region. 
Boundary analyvsis by definition follows the boundary of the 
profile. It can build only the relat ionships between parts of the 
boundary along the boundary. It can not build cross relations 
between parts of boundary across the interior of the profile. Such 



cross relations can be extremely importarrt . For example a 
boundary analysis may not differentiate between cases vShown in 
Fig. 1.1. In part (a), the boundary points A and B are quite far 
apart from each other. There is no significant relation between 
them. In part (b), A and B are close by; they denote a 
constriction in the shape. Such relations can not be directly 
detected in boundary analysis. It does not have a mechanism to 
store such, across the interior, relations. 

A transform, proposed by Blum on two dimensional shapes, variously 
called Skeleton, Iledian Axis Transf ormCMAT) , or Symmetric Axis 
Trans f orm(SAT ) has been reviewed in literature survey (ref. 
chapter 1 section 9.3). This transform is well suited for shape 
analysis. It fulfills all the requirements above reasonably well. 
It preserve.^ interaction acroSvS the interior along , with 
interaction along the boundary. It reduces a shape into a more 
familiar domain of line diagrams which may be easier to analyse. 

2 . Skeleton Def init ion J. 

2 . 1 Def init ion of the Skel eton of a Profile 2 . 

Compute the minimum distance of every point in the interior of the 

profile from the boundary of the profile. With every interior 

point, at least one point on the boundary is associated, which is 

the foot of the line of minimum distance. The skeleton is defined 

as the set of interior points which are associated with more than 

one point on the boundary. In Fig. 2.1.1, the points P , P and P 

1 2 3 

are associated with one point on boundary, points F , F , F 

12 3 

respect ively. Point P i associated with points F and F 


skeleton 





as distance F -P is the same as distance F -P 

4 skeleton 5 skeleton 

2 . 2 Pi f f er ent perspectives on the skeleton 

The points forming the set are physically contiguous. The skeleton 
is a continuous line diagram. The skeleton, along with the 
distance associated vith the constituent points, is completely 
equivalent to the profile. Hence the transformation is without any 
loss of information [9, 15 &. 25]. 


The minimum distances of the interior points can be thought to to 

form a field in the interior. The ridge points of the distance 

field form the skeleton [25]. In Fig. 2.2.1, the distance of point 

P is the largest among the distances of points P , P , P 

skel eton 12 3 

and P . P is a ridge point of the di.stanc6 field, 

skeleton skeleton 

The skeleton can alternatively be thought of as a ’quench line of 

prairie fire’. If a field congruent to the profile is .set afire at 

the boundary simultaneously, fire fronts advance from the boundary 

at the same rate. Wherever two fire fronts meet, a quench line is 

formed. The set of these, quench lines is the skeleton [15]. In 

Fig. 2.2.2, the rectangle V -V -V -V is the profile,. Fire fronts 

12 3 4 

start from the rectangle V -V -V -V at time t = 0 . At t = t’ , 

12 3 4 

fire fronts reach the position V ’-V ' -V ’-V ’, forming the quench 

12 3 4 

lines V -V ’ , V -V ’ V -V ’ and V -V ’ . The fire fronts advance 
1 1 2 2 3 3 4 4 

through time t’’ till they are completely quenched forming the 

skeletal lines V -A, V -B, V -B, V -A. This perspective on 

1 2 3 4 

skeleton gives rise to an important parameter of the skeleton, the 
quench velocity. It is the rate of formation of the quench line 
with respect to the distance(q) through which the fire front 
advances. 
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Thus the quench velocity, V = cls/dq = 1/Cdq/ds); where a is the 

<5 

length of the formed skeleton, i.e. distance along the skeletal 
lines. 

The shape can be deemed to be generated by a circular disc swept 

along the skeletal lines; the radiusCr) of the circle being the 

minimum boundary distance associated with the skeletal points 

[15]. In Fig. 2.2.3, circle C is swept along the skeleton. At 

point 0 on the skeleton, which is the centre of the circle (C ), 

1 1 
the circle has radius r . As it is swept to 0 , the radius 

1 2 

increases to r . This view of the skeleton is particularly useful 
2 

in shape analysi.s. 

3 . Discussion on Skeleton j_ 

The skeletons of some common polygons are illustrated in Fig. 3.1. 
Skeleton of a polygon consists of mainly lines emanating from the 
convex vertices of the polygon. As the line is traversed, the 

radius of the circle increases. At some point the circle grows so 

large that it comes into contact with another boundary segment. At 
this point the skeletal line forks, giving forth skeletal lines 
produced by combinations of the boundary segments. 

The concave vertices of the polygon may interact with a boundary 
edge, giving rise to curved skeletal segment.s. If two concave 
vertices int'eract, a linear skeletal segment is produced. The 
skeleton of a convex polygon consists of only straight lines. 

A non-self intersecting polygon without holes has a ’tree’ skeleton 
[15,25]. The tree is geometrically located in the plane of the 
profile. The nodes of the geometrical tree are the forking points 
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of the skeleton and the vertices of the polygon. The vertices form 
the leaf nodes of the tree. Its edges are the skeletal lines. The 
skeleton of a polygon with a hole is a, planar graph [ 9 ] , [ 1 5 ] , [ 2 5 ] . 

Skeleton is sensitive to sharp m i n o r 1> o xx n d a r y v a r,' 1, ri t .1 o n s . A s in a 1 1 

notch in the periphery can distort the skeleton of a rectangle 

drastically. In Fig. 3.2Cb), the notch V -V -V in the rectangle 

4 5 6 

V -V “V “V , has disfigured the skeleton of the rectangle. The 

12 3 7 

skeleton of the rectangle is shown in part (a) of the figure. A 
roinor projection produces extra long linevS in the skeleton due to 
interaiction between the concave vertices as' shown in Fig. 3.3. 
Hence the skeleton can not be directly used for .shape analysis. An 
approxiinat ion luethod is needed to get rid of the disturbances due 
to minor variations. 

Skeleton is based on the interactions of the boundary across the 
interior of the polygon. It does not consider interactions of 

boundary across the exterior of the profile. In Fig. 3.4, the- 
points A and B are far apart in part (a). In part (b), they are 
quite close by. These relations are not detected in the skeleton. 
An elaborate analysis is needed to encapsulate these interactions. 
This analysis can not be based on the skeleton. It is based on an 
analysis of the convex hull of the profile and its differential 
area between the convex hull and the profile. It is the bay 
analysis proposed in the complete prof i 1 e description , in chapter 
5 , s action 2 . 

The concept of 'Visibility of Boundary Elements across the 

Interior of a Profile/ has been expounded in 1 i t erature [ 1 , 11&32 ] . ^ 


32 





If the line joining two vertices of a profile lies fully within 

the profile then the two vertices are said to be visible to each 

other. This definition can be extended to any pair of boundary 

points. Skeletal interaction between a pair o£ boundary points is 

a more strict version of visibility. If a pair of boundary point.s 

interact to give rise to a point on the ' skeleton, they must 

be visible to each other. In Fig. 3.5, the boundary 

points P and P interact with each other to produce a skeletal 

point s/xhe skolotal circle of interaction is C. cantered at S. 

.the circle C lies fully within the profile, the chord P^-P^ 

must also lie within the profile. Hence the points P and P^ are 

1 ^ 

also visible each other. 

4 ■ Skeletal interaction 

Skeletal interactions in a polygon can be of three type.s : 

1) between two edges of the polygon, 

2) between a concave vertex and an edge, 

3) between two concave vertices. 

4.1 Skeletal interaction between two edg e s j_ 

The skeleton produced by two interacting edges i.s a straight line 
segment. It is the bisector of the angle between the two \ edges. 
The radius distribution along the skeletal line is linear. In Fig. 
4.1.1, along segment PQ of the .skeletal interaction between edges 
V -V and V ’-V The radius variation along PQ is linear from rl 

at pLnt P to r at point Q. The quench velocity is constant. 

V = S / r = sinC0/21. 

n w d^Toritv reduces with increase in the angle between edges. 

Quench velocity 


3 4 





As 0 -> 'tT, V -> 0, and 

q 

as 0 -> 0, V -> infinity, 
q 

The quench velocity for skeletal line of two parallel edges is 
infinite, l.e. the whole line forma simultaneously as the fire 
front advances from both the edges. 


4.2 Skeletal interaction between a concave vertex and an edge ^ 


The skeletal interaction between a concave vertex v and an edge e 

is a parabolic arc as depicted in Fig. 4.2.1. The radius increases 
* 

as the arc is traversed away from the concave vertex. The quench 
velocity reduces along the arc. 


For the parabolic arc shown in the Fig. 4.2.1, 
y**2 = 4ax ; 

q = X + a = y**2 / C4a) + a ; 

ds dx dy 

V = = aqrt( ( )**2 + ( — -)**2 ), 

q 

dq dq dq 

V = sqrtC^ 1 + Ca/(r-a))) . 

q 

V >1, and as s -> infinity, r -> infinity, and V -> 1. 

q q 

4 . 3 Skeletal interaction between two concave vertices ; 


The skeletal interaction is the perpendicular bisector of the line 
joining the vertices (Fig. 4.3.1). The radius distribution has a 
minimum at the midpoint of the line joining the vertices. The 
quench velocity has a maximum at the same point. 
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9 = sqrt Cq**2 - l**2/4) ; 


ds 1**2 / 4 + S**2 

V = = sqrt C ) ; 

q 

dq S**2 

V >1, and as s -> infinity, V -> 1 . 
q q 

4 . 4 Interaction strings and zones 

A skeletally active polygon boundary element, an edge or a concave 
vertex, interacts with one or more other elements. These 
interactions are sequential i.e. if we start from an extreme of 
the interaction field of an element, the element interacts with 
only one element at a time. As the skeletal line is traversed, a 
point is reached where the circle comes into contact with another 
element. As we sweep past this point the contact with the earlier 
element is lost and the element starts interaction with the next 
element. This continues till the other extreme is reached. It is 
illustrated in Fig. 4.4.1 with explanatory notes for the sequence 
of interactions. Thus the total skeletal interaction of an element 
is a ’string of interactions’. The skeletal tree is composed of 
individual strings of interactions of all elements. 


The 

skeletal string 

of 

segments in 

a string of 


int eract ions 

is 

a.£ain composed of the 

distinct zones 

of each interaction pair. 

In 

Fig . 

4.4.1b, s egment 

AB 

is the interaction 

zone 

of 

elements V 

1 

and 

e 

on skeletal line. 

BC 

between V 

and e , 

and 

so 

JL 

on . Segment 

PQ 

1 



1 

2 





is 

the interaction z 

one 

on edge e , 

for e 

and 

V 

. Angle AV B 

is 




1 

1 


1 

1 


the 

interaction zone 

at 

V , for V and e . 






1 1 1 

The area bounded by the interaction zone and the skeletal line 
associated with the interaction zone is the skeletal interactibn 
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ht A : 
B : 
C ; 
D 

E : 
F : 


Interaction 
int eract ion 
interaction 
interaction 
interaction 
int eract i on 


starts 
Q.-e.y ends; 
e - ends ; 

Q ends; 

(O ends; 

€ “fis ends ; 


int eract ion 
i nt eract i on 
int enaction 
int eract ion 
int eract ion 


(0 - starts 
starts 
starts 
^ -1S3 starts 
starts 


Ca) String of interaction of edge ’e’. 



At A : 
B : 
C : 
D : 
E : 
F : 


Interaction 
interaction 
int enaction 
int enaction 
int enaction 
interaction 


Vl-el starts 
Vl-e2 ends; 
Vl-e3 ends; 
VI -e4 ends; 
VI -e5 ends ; 
Vl-e6 ends; 


int enaction 
int eract i on 
interaction 
interaction 
interaction 


VI -e2 starts 
VI - 6.3 starts 
VI -V4 starts 
VI -e6 starts 
er-e6 starts 


(b) String of interaction of vertex V, 


Fig. 4.4.1 Strings of Interaction. 



area of zone. In E^ig. 4.4.1Cb), the area ABPQ is the skeletal 

interaction area for the interaction zone of the V -e interaction 

1 1 

for the edge e . Similarly, the area AV B is the skeletal 

1 1 

interaction area for V -e interaction for vertex V . The total 

11 1 

area V APQB is the interaction area for the skeletal line AB. The 
1 

sum of all the skeletal areas asvsociated with an edge is the 
skeletal interaction area of the edge. Similarly skeletal 
interaction area is defined for a vertex. 

5 > Concavit i ea in polygon and their ef f ects on the skeleton 2. 

It can be easily seen that the skeletal lines produced by 
interactions of the concave vertices are different from those 
produced by two polygon edges. The skeletal lines of interaction 
between two concave vertices are geometrically straight lines but 
the radius distributions and quench velocity are different from 
those for skeletal lines of interaction between two polygon edges. 

The skeletal tree can be deemed to be composed of subtreevS of 
different types, the convex skeletal interactions and the concave 
skeletal interactions. The convex skeletal interactions are 
homogeneous. They form strictly linear geometrical trees. In 
contrast the concave skeletal interactions are not homogeneous. 
The interactions between two concave vertices are different from 
those between an edge and a concave vertex. The concave 
interaction subtrees are composed of concave interact ion strings * 

The properties, such as quench velocity, radius distribution , 
geometrical size, of the convex and concave subtrees and strings, 
and their distribution through the skeletal tree are the important 
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features for shape analysis. 


6 . Properties of skeleton ; 

The features of skeleton which are important for shape analysis 
are defined below. 

6 . 1 Characteristic lengths t_ 

Various characteristic lengths are used in shape analysis. They 
are, 

1) naximum of the distances between two vertices of the polygon 

CL ), 

1 

2) Maximum radius in the skeleton CL ), 

2 

3) Maximum length in the skeleton + maximum diameter in the 

skeleton CL ). 

3 


6 . 2 Properties of concave interactions j_ 

6.2.1 Concavity measure t_ 

The exterior angle of a concave vertex is a good measure of the 

concavity of the vertex. In Fig. 6. 2. 1.1, a part of a profile with 

edges e, e , e and a concave vertex V is shown. 0 is the 
1 r e 

exterior angle between e and e . 0 is the angle between 

1 r p 

perpendiculars to e and e . It is the angle of interaction for 

1 r 

the skeletal interaction between V and e. 

In Fig. 6.2. 1.1, 0 = "XT - 0 ; 

P e 

concavity measure C = Fuzzy_normalis6C0 , TT) . Cref. chapter 

cv P 

1 section 6 for the function f uzzy_normalise) . 

6.2. 2 Peripherality j_ 

A profile may have a lot of minor concavities giving it a jagged 
boundary. Such jagged boundaries must be smoothened before a shape 
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analysis as they do not contribute much to shape. On the contrary, 
they unnecessarily complicate the solution space and hence the 
shape analysis process. 


An essential feature of the jagged edges is that the edges 

adjoining the concave vertex are small. But, this alone can not be 

the criterion for detection of jaggedness or ' per ipheral i ty ' of a 

concavity. A smooth and curved boundary segment may be composed of 

a large ‘number of small lines. The peripheral concavities are 

dist inguished by their long interaction strings and pre.gence of 

large radius sections interacting with boundary elements straight 

across the interior. In the Fig. 6,2.2. 1, the notches V -V -V in 

1 2 3 

part (a) and V -V in part (b) have geometrically long interaction 
2 3 

strings. The interaction starts with small radius r . The radius 

1 

increases till it becomes much larger than r , upto the R in 

1 max 

the figure; then it reduces again to r , which is again quite 

2 

small compared to R . Presence of long interaction strings is 

max 

also not a reliable criterion. Only the long range interactions 
across vast distances of interior areas give a strong indication 
of peripheral concavities. Some per'ipheral concavities may be 
hidden from the interior by presence of another peripheral 
concavity. In this case, the overlying peripheral concavity is 
identified as such first. The concept of per ipheral ity is being 
used in the profile decomposition algorithm, developed in the next 
chapter. In the algorithm, the decompoBit ion is achieved by 
deleting concavities ( concave vert ices ) one by one. Hence, when 
this concavity is deleted, the other concavity gets identified as 
a peripheral concavity. In Fig. 6. 2. 2. 2, the per ipheral concavity 
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V -V -V IS occluded by the peripheral concavity V -V --V . In this 
3 4 5 1 2 3 

case, V -V -V will be identified as a peripheral concavity. After 


12 5 

it is deleted, V -V -V will 

3 4 5 

concavity and deleted. 

Peripheral ity may be defined as 
the edges adjoining the concave 
radius section. Determine the 
interaction 'string. In Fig. 6 
radius , then 

peripheral ity , P = Cvery_sraall 

pr 

large(R ). 


be identified as a peripheral 

the conjunction of smallness of 

vertex and existence of a large 

argest radius in the concavity 

2.2.1, if R is the largest 

max 

edge ) OR very_smal 1 ( edg e )) AND 
1 2 


max 

6.2.3 Importance of a concavity j_ 

A concavity which penetrates the interior and divide-s the profile 

into two clear parts is an important concavity. On the contrary, 

minor peripheral concavities may interact with a large section of 

the boundary, with large radius, but they do not cleave the 

interior of the profile. Such a concavity V -V-V is depicted in 

1 ■ 2 

the Fig. 6. 2. 3.1. The general trend of the radius(r) versus 
distance along skeletonfs) characteristics of concaviti-es within 
the interaction zone give an indication of cleavage. The rC.s) 
curve is concave upwards. In Fig. 6.2.3.1b, the rfs) curve of the 
interaction of the concave vertex V, which cleaves the interior of 
the profile, is illustrated. Sharper the concave nature of the 
rCs) curve, clearer the cleavage. But this function fails to take 
care of situations where a smooth and curved concavity may be 
composed of a sequence of short line segments (Fig. 6.2 . 3.2a) , and 
where a thin bridge connects two large polygons (Fig. 6 . 2 . 3 . 2b) . 
In both these cases, though the interior is cleaved, the rCs) 
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curve, as shown in Fig. 6. 2. 3. 2, is not sharply concave upwards. 
This is due to the local nature of the interaction string of a 
concavity. In situations given in Fig. 6. 2.3. 2, the interaction 
string of the concavity reflect the situation in the local 
vicinity of the concave vertex. Cleavage is a global property. A 
global r(s) cxirve, of concave natxire, truly reflects the cleavage. 
A global rfs) function may be extremely varied. Its study is would 
be a daunting task. Hence a scheme is proposed for importance 
evaluation which is not exact but it is very simple. 

Cut the skeleton into two at the point of the lowest radius in the 

concave interaction string. Evaluate the largest radius on both 

part.s of the skeleton. Evaluate the sequence of largest radius to 

the left CR the smallest radius R , and the largest radius to 
1 min 

the right (R for a concave variation i.e. ’large-small-large’ 
r 

variation, as a conjunction of fuzzy predicates ’large’ and 
’small’ applied to their respective arguments. The resulting fuzzy 
set membership grade gives the 'importance’ of the concavity; that 
is , 

Importance = LargefR ) AND smallCR ) AND largeCR ). 

1 min r 

It is suggested that the fuzzy conjunction be implemented as the 
average of the three fuzzy membership value. 

6.2.4 Proximity of two concavities 

Proximity of two concavities is NOT of the relative distance 

between them, with respect to characteristic length L . The 

3 

relative distance between the two concavities is the comparison of 
the sum of their radii at the closest points in their interaction 
strings and the distance along the skeleton between their 
interaction strings. For example, in Fig. 6. 2.4.1, the distance 
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between V and V is the distance ABCDEF. BCDE is the skeletal 
12 

line, whereas AB is the radial distance along the right extreme of 

the interaction zone of vertex V with the opposite edge e. 

1 

Proximity P = NOT ( Fuzsy_norinal is e Cd,L )). (ref. chapter 1 
X 3 

section 6 for the function fuzzy_normalise) . 

Note that distance along skeleton is the distance between the 
nodes of a tree and hence can be easily determined. 

6 . 3 Properties of convex int eract ions 


6.3.1 Forks j_ 

The skeleton, as observed earlier, is a tree coraposed of 

geometrical lines. The nodes of the tree are the forking points of 

the line diagram. The.9e nodes will be referred to as the 'forks' 

of skeleton in the further work. In Fig. 6.3.1, f and f are the 

1 2 

forking points of the skeleton. They play an important role in the 
processing of the convex interactions. 


Vi 



Fig. 6.3,1 Convex Skeletal Interaction Lines. 
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The ’Size’ ‘function in the last formula may be the 
’ Local_Si ze ( f ork ) ’ or the ’ Global__SizeC f ork.) ’ , in accordance with 
whether the largeness is to be local or global. 

6.3.4 Proximity of Forks j_ 

Proximity of two neighbouring forks is quantified as the length of 
the NACSIL connecting them compared to the locally largest radius, 
i.e. the larger of the radii of the two fork, or the globally 
largest radius in skeleton. 

If 1 is the length of the NACSIL between forks f and f , and 

1 2 

Rg and R1 are the global and the local maximum radii 

max max 

respectively, then 

Local_Proximity C.f ,f ) = Smal 1 ( Fussy_normal is e C 1 ► R1 ))» 

1 2 max 

Global_Proximi ty C f ,f ) = Smal 1 ( Fussy_norraal i s e ( 1 , Rg )). 

1 2 max 

6.3.5 Intensity of CSI L 

Intensity of the CSILs is a property indicating the importance of 
the CSIL in the skeleton. It combines the radius and length of the 
CSIL. It also depends, although indirectly, on the quench 
velocity. 



Fig . 6 .3 .5 Intensity and Skeletal Interaction Area, 
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Intensity is 


defined a.s the area associated with the skeletal 


interaction which generated the CSIL, con^pared to the total 

of the polygon. In Fig. 6.3.5, the area V -P -f -P is 

1112 

interaction area for the ACSIL V — f ; and the area P -P ~P ~P 

1 1 5 3 2 < 

the interaction area for the NACSIL £ -f . 

1 2 

The intensity ’I’ is given by 

I = Fus 2 y_norinaliseC interaction a.rea, area of 
polygon) . 

6.3.6 Int ens ity of forks Inten.sity of a fork is the sum of 
intensities of the intensitie.s of all the CSILs converging at 


area 
the 
i 3 

the 

the 

the 


fork . 



Chapter 3 


PROFILE DECOMPOSITION 

1 . Ob J»ctiv 

The proposed nestine method is based on shape description and 
matchine of the shapes on the basis of description, for profile 
allocation. A description scheme encompassine all possible shapes 
requires some structurinfi and 'classification of the shapes. 
Classification of all shapes is a highly complex task. The 
possible shapes are too numerous, too diverse and they have too 
many refiions of complex ranges amonfi the extreme cases. Hence as a 
first step in a structurins process, the profiles are decomposed 
into constituent convex polyeons. Convex polyfions are more 
amenable to imposition of classification. Shapes are described in 
terms of the constituent convex polygons (primitives), and their 
topological and geometrical relations. 

The basic requirement of the decomposition scheme is that the 
decompositions should be consistent with human perception of the 
shape. Various decomposition schemes have been proposed in 
literature. All the schemes have been briefly reviewed in the 
literature survey (ref. section 9, chapter 1). The results of some 
representative decomposition algorithms Pavlividis £25], Ahuja £1] 
and Chazelle £11], are discussed here., The objective of these 
schemes is to decompose the profile into a minimum number of 
convex polygons. These schemes are exact in nature. Hence small 
disturbances in the boundary lead to large and undesired 
variations in the decompositions. It is illustrated in Fig. 1.1. 
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Fifi. 1.1 Need for Smoothenine before Polyfion Decomposition. 



(a) Pavlividis alfiorithin. 



Ca) Ahuja aleorithm. (b) Chazelle alftorithm. 

Fig. 1.2 Results of Some Decomposition Algorithms. 
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Part (a) shows a rsctancle with a very siaall notch in the 

boundary. The small notch induces decompositions as shown by the 

dashed lines. A similar effect is seen in the part (h) also. The 

small notches should be smoothened out to maintain consistency 

with perception. In the decomposition algorithms mentioned above, 

for some simple profiles, the decompositions obtained are not 

consistent with human perceptions. Some such profiles are shown in 

Fig. 1.2. Part (a) depict® the result of the application of the 

Pavlividis algorithm [25]; Some of the constituent convex polygons 

are shaded. It can be easily seen that the decomposition does not 

yield much information about the shape features. Especially the 

decomposition in the left half of the profile, producing the 

quadrilateral V -V -V -V seems too artificial. In part (b), the 

3 4 11 12 

results of application of Ahuja decomposition algorithm to the 

profile V -V -V -V --V are shown [1]-. The decomposition into three 
1 2 3 4 5 

triangles again is not natural. A decomposition in to two 

quadrilaterals, by dropping a perpendicular from vertex V on to 

3 

the edge V ~V , would be much better. The results of the Chazelle 
5 1 

algorithm are t in t’la (c) [11]. The presence of a long 

cut running through th. interior of the profile again belies the 
perception. The drawbacks in the Pavlividis and Ahuja algorithms 
are mainly due to the fact that only vertex to vertex cuts are 
used for decompositions. In the Chazelle algorithm, some of the 
drawbacks of the other two algorithms .have been overcome. But 

still, the decompositions may seem artificial in some cases. This 
drawback arises due to a strict adherence to the objective i^e. 
decomposition into minimum number of polygons. 



2 . Important Concepta for Decoropoaition Problem j. 

Ue propose modifications in the formulation of polygon 
decomposition problem. It is not sufficient to obtain solutions 
which will decompose the profile into a minimum number of convex 
polygons. The optimisation criterion should also ensure that the 
constituent polygons are such that they are easily identifiable in 
the profile. It' can be seen from Fig.s 1.1 and 1.2 that the 
decomposition schemes do not guarantee it. The resultant 
decompositions which might be minimal in terms of the number of 
polygons. But a shape description based on these constituent 
polygons and their interrelations, will not be useful for nesting. 

The ’ease of identif iability’ of constituent polygons depends upon 

the way the profile is 'cut up’ in decomposition. The ’cuts’ are 

extraneous to the profile. If the cuts are kept to a minimum it 

will increase the identif iability of the constituent polygons. It 

is illustrated in Fig. 2.1. The cut A-V -B produces three convex 

3 

polygons which are more easily identifiable than the convex 

polygons produced by the other cuts. The length of the cut A-V -B 

3 

is smaller than that of the other cuts. The convex polygons 

generated by the cut A-V -B have boundaries which are more exposed 

3 

to the exterior of the profile than the boundaries of those 
generated by the other cuts. It may be inferred from these 
examples that for a decomposition to produce easily identifiable 
constituent convex polygons, the boundary of the constituent 
polygons should be largely composed of the boundary of the 
profile. To express it in different words, the constituent 
polygons should be ’maximally exposed to the outside ' ; or. 
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Fifi. 2.1 Identiflability o 


Decompositions based on the 
identifiability of consituent 
polyfions . 

Decompositions based on the 
minimum number of constituent 
polyeons . 

Constituent Polygons. 



(c) Approximate decomposition (d) Exact decomposition 

* ■ V 

Fifi. 2.2 External Visibility - 'Lagoon’ 
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expandine further, the constituent polyfions should be 'maximally 
visible’ form far, in the plane of the profile. 

The concept of 'maximum external visibility' alone with ease of 
identif lability., also incorporates another concept. A profile may 
have a laree 'laeoon' almost closine the profile onto itself (ref. 
Fifi. 2.2Ca)). Thoufih it is possible that this laeoon may be used 
in nestine. it is quite unlikely that there will be another 
profile with a narrow 'neck' with a laree body attached to the 
neck, fittine the laeoon. The laeoon is most likely to be wasted, 
especially so if its boundary is irreeular. For example, the 
probability of another profile which may fit in to the laeoon of 
the profile in Fie 2.2Cb) is quite low. Since the proposed shape 
description does not consider the effects of the cross relations 
of profiles the possibility of nestine of the laeoon can be 
neelected. Hence an approximate decomposition scheme closine the 
lagoon and neglecting its presence would be better suited for 
nestine (ref. Fig. 2.2c). The exact decomposition given in part 
(d) of the fieure is unnecessarily complex. It should be noted 
that the external visibility of the lagoon, from a large distance, 
is low. 

The quantification of the concept of 'external visibility from 
far’ would require probabilistic or possibilistic analysis 
considerine the random possibility of nesting of some other 
profiles into the concavities. It is not attempted in the present 
work. A simplified version 'external visibility from close 
quarters’ is used. It neglects the 'lagoon cases’ and their 
offects altogether. It embodies the ’ease of Identifiability’ in 
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terma of 'cut mininiality’ . All hitherto references to external 
visibility should be taken to mean this concept. 

Cut minimality is represented by the minimality of the total 
lenfith of the cuts . External visibility can be represented by the 
ratio of total lensth of the line diafiram formed by the profile 

and the cuts to the total leneth of the orieinal line diaeram i.e. 

the perimeter of the profile. 

Another essential feature of the decomposition into convex 

polycons is that some minor variations in the boundary should not 

induce larfie discontinuous chanfies in decomposition. In Fi£. 

2.3Ca)» the edees V -V and V -V are collinear. A line joining V 

2 3 5 6 3 

and V completes the decomposition. In part (b) , the two edaes are 
5 

slightly offset. Two cuts are needed to decompose the profile. as 

shown. Thus a very minor variation has induced a major chanse in 

the decomposition. In a decomposition algorithm, the minor 
variations should be smoothened out and approximated. They may not 
be neglected altogether. But their effects should be reduced 
ensuring continuous variations in decompositions. This is 
especially important for small concavities on the boundary. Small 
concavities such as very shallow bays and very .sharp notches, 
which are not likely to be affect the nesting efficiency much, may 
bo covered by 'straightening up' the profile boundaries. In Fig 
2.4, a sharp notch is covered up by an approximating edge. Also, a 
sequence of minor concavities, if not approximated, would give 
rise to a large number of small constituent polygons, thereby 
destroying usefulness of the decomposition of the profile. 
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V2-V3-V5-V6 are colllnear. V2-V3 & V5-V6 are sllfihtly 

offset. 

(a) (b) 

Fifi. 2.3 Minor variations in boundary affectinfi decomposition. 



Fifi. 2.4 Strai£ht enins Up / Approximation. 

It is important to note that if small projections are simply 
ignored, it could be fatal for neat ins. Approximations should be 
done only on concavities. This implies that the approximations 
will strictly be external envelopes of the profile. 

The extraneous areas introduced by ’straishtenins up’ could be 
taken to be a measure of approximation. As the approximation Is 
always an external envelope of the profile, the possibility of 
spurious ’sood approximations’ due to cancellation of positive and 
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nefiative approximations, is absent. Hence ratio of orifiinal area 
to area of approximation £ives a aood measure of the accuracy of 
the approximation. 


3 . Formal Statement of the Problem x. 

The profile decomposition problem can be stated as follows : Given 

an arbitrary polyaon P approximate it by an aaa^eaate of convex 

o 

polyaons P , i = l...n, so that the objective function, 

i 

’aoodness’, GCv,a) is maximised, where ’v’ is a measure of the 

external visibility of the constituent polyaons, and 'a’ is a 

measure of the accuracy with which P is approximated. Note that 

o 

the approximation is by external envelopes. 


External visibility measure 'v' is taken to be 


PerimeterCP ) 


Perimeter(P ) + L 
o c 

where L = total lenath of cuts, 
c 

Accuracy measure 'a' is taken to be 

Area(P ) 

^ 1 + — — — — ^ — — ) . 

n 

y~" AreaCP ) 
i = i i 

The objective function, ’aoodness’, G = v * a . 

Explanation of Measure of Accuracy Formulation t, 

Durina shape description, if the extraneous area introduced by 
approximation is described as a bay then there exists a 
possibility that this area could be utilised in nestina. Thus the 
extraneous area may not bo completely wasted. If it is utilised 
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fully, the approximation le not really an approximation and the 
accuracy measure should be one. If the area is not utilised at 
all, then the approximation accuracy could be taken to be the area 
ratio. In absence of any information on the possibility of 
utilisation of this area, the chances of both the cases are taken 
to be the same and an averaee value is taken. 

If on the contrary, the extraneous area is not described as a bay, 
then it will not be utilised. In this case the approximation 
accuracy should be taken as just the area ratio. 

4 . Exploration of the Solution Space 

This problem could be viewed as a problem of deletion of the 
concavities in a profile by a judicious combination of cuts and 
approximations. In absence of approximations, atleast one cut 
should originate at or pass through a concave vertex. In figure 
4.1, some cuts required to delete the concavity at the vertex V 
are shown. More than one cutting line will be required if the 
cutting line does not fall within the reflex angle of the vertex. 
In such a simple case, a procedure of concavity deletion one at a 
time, with selection of the locally best solution at each deletion 
could be sufficient. The possibility of intersection of two cuts 
increases the complexity of the solution space and thus rules out 
such simple solution methods. 

Approximations add another dimension to the complexity of the 
solution space. The boundary lines of constituent polygons now 
need not pass through, the concave vertices. In figure 4.2, a 
portion of the boundary of a profile is shown as the edges 

connecting vertices V to V . The approximate decomposition line, 

• ’ 1 ^ 8 
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Fifi. 4.1 Cuts at a Concave Vertex. 



j- 


Fifi. 5.1 Ran«e o 


i Possibilities at a Concave Vertex. 



shown dashed in the fieure, does not pass throufih the concave 
vertices V4 and V6. Hence the solution space, with this provision, 
becomes the space of all possible geometrical graphs ~ geometrical 
line diagrams, produced by the combinations of the line segments 
in the plane, which either lie within the profile or terminate on 
the profile boundary. 

5 . Method of Solution 

A meaningful mathematical formulation of the problem for 
application of the classical optimisation methods seems 
impossible. A combination of heuristic and combinatorial search 
methods is a promising avenue for optimisation. The search over 
the solution space could be divided into search over the cuts and 
approximations for each concavity for the concavity deletion, and 
a search over the sequence of deletion of concavity. The advocated 
algorithm for decomposition is based on some ordering of the 
concavities for deletion and deletion of one concavity at a time. 
With each deletion, the profile is subdivided in to two or more 
subprofiles which may or may not be convex. The decomposition is 

\ -i 

complete when every subprof il^ is a convex polygon. 

The search space in each concavity deletion is uncountably 

i ‘ ■ 

infinite. In figure 5.1, a portion of the boundary of a profile is 

shown as edges from V to V . A number of approximate cuts are 

1 8 

shown, . all of which may delete the concave vertex V6.It has 
discrete regions wherein the problem characteristics change 
disGontinuously from region to region. But they vary continuously 
within each region. The proposed solution method is to locate 
’most likely’ cuts and approximations in each region using 
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heuristics. Then the problem could be solved by heuristically 
fiuided combinatorial search method. 

Different concavity deletions may have cross effects fiivinfi rise 

to the problems such as cut intersections, and interaction of cuts 

and approximations, and interaction between two approximations. 

They are illustrated in figure 5.2. In part (a), cuts cutl and 

cut2 intersect. The problem is which of them is to terminated at 

the intersection point, and which one should be extended to the 

boundary. In part (b), the approximation is done for the new 

concave vertex at V5 generated by the cut, which has been applied 

earlier. The cut now becomes redundant. In part (c), a chain of 

concave vertices V to V is being straightened up by 

2 7 

approximations at various combinations of the vertices. These 
approximations intersect and create a few more possible 

approximations. These complications are taken care of by heuristic 
rules; and the different deletions are effectively isolated. 

The process can now be modelled as a generative search method. The 
search nodes at each step are generated by rules for cuts, 
approximations and their cross relations. The solutions are 
evaluated by heuristic evaluation criteria. The heuristic 
evaluation criterion is the 'usefulness' of the solution. The best 
solution node is explored first. 

If a dead end is reached or the goal state is not sufficiently 
good enough, then the solution tree is backtracked. A dead end is 
a solution with a low objective function 'goodness 'value. ^ 

For complex shapes, likely to be encountered as cotRplementary area 
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in nestinfi, the solution tree may become very large . Large amount 
of backtracking may become necessary to obtain a good solution. If 
the solution is quite bad then the process may spend time 

unnecessarily searching and backtracking over similar solutions 

' \ 

which may not offer any aleni|lcant Improvement. To avoid such a 
condition it is necessary to identify the critical decisions, 
which when backtracked, would offer significantly different 
solutions. The factors affecting the 'criticality’ of a decision 
are the importance of a concavity, degree of difference of the 
solutions, and the local 'goodness’ of the possible solutions. 

Multiple shape descriptions are necessary for nesting. Also the 

simplifications in the optimisation problem erode the confidence 

! 

in the final solution. The goodness of the solution depends upon 

.V 

the heuristics, which can’t be guaranteed to work equally well 
across the whole range of shapes. Hence it is necessary to obtain 
multiple solutions in the decomposition process. 

The multiple decompositions should be sufficiently different from 
one another. They should be comparably 'good'. This can be 
achieved by developing all of the most critical solution nodes. 
Then finally the best comparable few solutions can be reported as 
the decompositions. 

In the fuzzy spawning paradigm (ref. chapter 1, section 6), the 
pruning of the decision alternatives for decomposition is based on 
the high usefulness. Additional pruning is introduced in the form 
of the backtracking control. The goal state evaluation is the 
heuristic function goodness of decomposition. The backtracking 
control is based on the difference criticality. The heuristic 
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functions or the fuzzy qualitative predicates are the coodness, 
usefulness, criticality etc. 

6 . Digcussiona on Concavity Spectrum and Concavity Deletion j. 
Concavities in a profile may assume diverse shapes and sizes. 
These concavities also interact with the non-concave regions of 
the profile and with other concavities in the profile. These 
features give rise to a complex range of profiles. An attempt at 
spanning the range of these features is presented below. 

6 . 1 Shape and Size of a Concavity ^ 

Here, only an isolated concavity with smooth monotone borders is 
considered, i.e. there are no intervening local projections or 
second order concavities on the border. 

1) General shape of the concavity ; A sharp notch (Fig. 6. 1.1 (a)), 
a gentle bay (Fig. 6.1.1Cb)), a vary shallow depression (Fig, 
6.1.1CC)). 

2) Nature of the boundary : Its boundaries may be smooth (Fig. 
6.1.2Ca)), or composed of long straight line segments(Fig. 
6.1.2Cb)). 

3) Size with respect to the other profile features : The concavity 
may be dominating - large, or negligible - small. 

6 . 2 Multiple Concavities 

1) Distribution of features listed in section 6.1, over the 
concavities. 

2) Consecutive concavities along a border which may give rise to a 
concave curve (Fig. 6.2.2). 

3) Concavities interspersed with various different convex regions. 


65 




6 . 3 Gcowatrical Arrangeinentg of Concavitieg in the Border ± 

1) The sequences described in section 6.2 may be restricted in one 
region of the profile boundary. 

2) Sequences scattered over the profile boundary. 

6 . 4 Geometrical Arrangements of Concavities Across the Interior 

1) Skeletal interaction among concavities giving rise to necks 
(Fig. 6.4.1). 

2) Skeletal interaction between concavities and convex regions of 
the profile (Fig. 6.4.2). 

6 . 5 Concavity Deletion ± 

The concavity deletion process should be capable of handling the 
whole spectrum of concavity possibilities. We propose two 
disparate methods for the concavity deletion. They are, ’Cuts' and 
’Approximations ’ . 

ye propose to delete concavities one at a time. These concavities 
are the local concavities between adjacent edges. The cases of 
concavities with smooth curved boundaries (ref. Fig. 6.1.2) are 
handled by sequences of successive applications of the two 
deletion methods. The sequences are generated by the successive 
nodes of the generative space search, that is, each concave vertex 
in the profile is deleted independently. 

The vertices of profile with interior angle greater than TT are 
the concave vertices. Each of these concave vertices is deleted by 
application of a ’Cut’ or an ’Approximation' . 
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are the local concavities between adjacent edges. The cases of 
concavities with smooth curved boundaries (ref. Fig. 6.1.2) are 
handled by sequences of successive applications of the two 
deletion methods. The sequences are generated by the successive 
nodes of the generative space search. That is, each concave vertex 
in the profile is deleted independently. 

The vertices of profile with interior angle greater than TT are 
the concave vertices. Each of these concave vertices is deleted by 
application of a 'Cut' or an 'Approximation'. 
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Fia. 6.4.1 Interaction between two Concavities 



Fie- 6.4.2 Interaction between Concavity and a Convex Reeion 



Approximation of concave vertex V 
Fie- 6. 5.1.1 Definition of an Approximation. 
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6.5.1 Approxiroatione j. 


It can be seen from discussion eiven above that profiles can have 
ereat variations and may bo extremely complex. A profile may have 
a large number of major and minor variations in the boundary. The 
minor and not so minor disturbances would tend to clutter up a 
decomposition solution with a large number of constituents. It 
will reduce the usefulness of the shape description. This 
situation can be avoided by a careful use of the approximations. 
The minor concavities are straightened up or deleted by 
’enveloping’ by approximations. 

Approximations, as a method of concavity deletion, are in the form 

of ’straightening’. In Fig. 6. 5. 1.1, the vertices V and V , 

1 r 

adjacent to the concave vertex being deleted, V, are joined by a 
straight line and thus the concavity is deleted. A sequence of 
cuts and ’straightening up’ approximations in the decompositioning 
process, illustrated in Fig. 6. 5. 1.2, may give other types of 
desirable approximations. 

Approximations of different concavities may overlap or otherwise 
interact with each other (ref. Fig. 6. 5. 1.3). The concavity 
deletion order becomes critical when the approximation interaction 
are feasible. 

6.5.2 Cuts i 

Cuts are taken to originate at the concave vertices. A cut, if it 
lies in the reflex angle of the concavity, consists of only one 
line. Else it consists of two lines (Fig. 6.5.2.1). The cuts are 
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Approximation 



Fifi. 6. 5. 1.2 Sequence of Approximations and Cuts. 



Fifi. 6. 5. 1.3 Overlapping and Intersecting Approximations. 



Fig. 6.5. 2.1 A Cut at a Concave Vertex. 



generated by application of the heuristic rules stated in section 
7 . 

A cut could theoretically, lie anywhere in the interior angle. The 

heuristics generating the cut should locate the most probable 

regions of the interior angle which would give maximal external 

visibility i.e. minimal length cuts. We propose that the skeletal 

interactions give a good basis to locate these regions. These are 

the skeletal interactions of the concave vertex being deleted. In 

another words, the concave skeletal interaction strings of the 

skeleton are being replaced by convex skeletal interaction linos. 

If a ray is passed from the concave vertex till it intersects the 

boundary of profile, and scan the whole range of the interior 

angle, the length of the ray and its variation with the angle 

would yield a lot of Information about the possible cuts. The cut 

should be located in the minimum length regions or in the regions 

where the characteristics of the length variation undergo major 

changes. These regions are roughly located by the skeletal 

interactions. In Fig. 6. 5. 2. 2 the skeletal interaction of the 

concave vertex V is shown. A cut should be located in the zone' 

0 

delimited by the angles 0 to 9 . With in this zone, 03 completely 

2 5 

and regions of 04 are ruled out if the cut is to end on the 

. y . 

skeletal interaction regions of vertex V on its interacting 

0 

edges, i.e. segments AB on edge V -V and segment CD on edge Y - 

,12 5 

V . Thus, it can be seen that the skeletal interaction regions of 
6 

the concave vertex, on the edges interacting with the concave 
vertex, are the regions where the foot of the minimum length cut 
are likely to lie. 


71 





Apart from passlnc rays at discrete an£ular interval , checkinfi the 
intersections and scanning the interior angle, there is no 
alternative to skeletal interaction for exploration and study of 
cuts. The process of passing rays would be computationally too 
expensive. The skeletal interaction inherently takes care of whole 
diverse concavity spectrum. Thus, to minimise cut length, cuts 
should be based on the geometrical regions of the concave skeletal 
interaction strings. 

There may exist situations where an interaction between the 
concave vertex and a boundary segment is just missed; e.g. vertex 
V and edge e in Fig. d.5.2.3. A good and viable decomposition 
could result from these missed interactions, as is apparent in the 
figure. The identification and an appropriate handling of these 
cases should be done on the basis of a fuzzy or possibilistic 
definition the skeleton. It is beyond the purview of the present 
work. 

Exploration of the Skeletal Interaction Regions and Formulation of 
Cut Rules ±_ 

A concave vertex interacts with a sequence of profile edges and 
other concave vertices. These separate interactions produce 
disparate regions with differing characteristics. Each of these 
regions should contribute atleast one cut to the set of solution 
nodes. In each interaction region, cuts connecting the concave 
vertex with the extremes and the centre of the interaction region, 
and the minimum distance cut should be generated. 

If the cut lies outside the reflex angle of the concave vertex. 
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then another cut line on the other side of the reflex ancle is 
generated by application of the same heuristic rules (ref. Fig. 
6 . 5 . 2 . 1 ) . 

6.5.3 Cut Interactions and Cut Modifications 

Cut Interaction Cuts taken during a particular concavity can not 

interact with each other. They are exclusive of each other. In 

Fig. 6. 5. 3.1, the skeletal interaction and the interaction regions 

of the concave vertex V are shown. The cuts C , C and C lie in 

12 3 

the angular interaction regions of V with edges V -V , V -V and 

12 2 3 

V -V respectively. Hence they can not intersect with other. After 
3 4 

a concavity is deleted, the profile changes. It may be replaced by 

an envelope, or it may be decomposed into multiple polygons. The 

skeleton of the new profile is determined afresh. In the new 

skeleton, the interaction zones are also different. Thus a cut may 

interact with cuts taken in a previous concavity deletion; e.g. 

cut taken in deletion of vertex V and cut taken in deletion of 
1 12 
vertex V in Fig. 6. 5. 3. 2. 

2 

Cut Modifications ± A cut modification results from an interplay 
between a cut and a subsequent approximation across the interior 
edges created by the cut (Fig. 6. 5. 1.2). (Interior edges are 
those which lie inside the original prof i 1 e boundary .) In such a 
case the earlier cut is obliterated and the new cut resulting from 
the approximation is retained. 
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6 . 5 . 4 Hauriat ic Evalmat ion of tha Intarinadiat a Solutions j. 


The ’fioodnass' measure of the intermediate solutions fisnerated by 
the heuristics could be used as the heuristic evaluation function 
'usefulness'. This would turn the decomposition process Into a 
hill climbina process with its attendant pitfalls like foothill 
climbina< To overcome these drawbacks, the 'aoodness' measure is 
modified by the 'concavity reduction measure’. This introduces the 
possibility of not followina of the locally best path every time. 
Hence, with a eood heuristic intermediate solution evaluation 
function, foot hills may be avoided. 

As the name sueeasts, the concavity reduction measure is a measure 
of how many concavities are baina reduced by the cut and to what 
extent. It may be noted that a cut might only reduce a concavity, 
and might not delete it. Also, more than one concave vertices may 
be affected by one cut, as the cut C in Fig. 6. 5. 4.1. This factor 
should be designed to increase the usefulness of a cut joining two 
concave vertices. It should also enhance the usefulness of the 
cuts which reduce the concavity to a large extent and those which 
completely delete the concavity. The exact evaluation of the 
concavity reduction measure is given later, in section 7.1. 

6.5.5 Criticality of a Concavity Deletion 2 . 

The criticality of a deletion is an indication of how different 
the possible decomposition results are from each other. It also 
depends on the importance of the children solutions. The useful 
solutions, which may be developed further, if they are drastically 
different from one another, may give rise to highly differ eht 
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decompositions. Hence these solutions are critical for the 
generation of multiple and different decomposition solutions. 
Multiplicity of decompositions is necessary to generate multiple 
and varied descriptions for a profile. The success of shape 
matching approach to nesting depends upon the ability to generate 
multiple descriptions. 

During the decomposition process, at a particular concavity 
deletion, from the possible solutions, the ’most useful’ solutions 
are shortlisted. All the most useful solutions can bo checked 
against the solution with highest usefulness for difference. The 
criticality of the deletion is the largest difference and the 
importance of the concavity. 

The difference between two deletion solutions can be a difference 
between a cut and an approximation or difference between two cuts. 
A cut and an approximation are absolutely different. The 
difference between two cuts depends upon whether the two cuts are 
1) vertex to edge cuts or vertex to vertex cuts, 2) one line cuts 
or two line cuts, 3) angular separation between the cuts. The 
first two sources of difference are combinatorial in nature, 
whereas the last one is continuous. A heuristic difference 
evaluation function combining these sources of difference and 
returning a difference measure needs to be formulated. The 
combinations of first difference can be given weightages. The 
difference is the weighted sura of the normalised angular 
separation. The exact formulation of the heuristic difference 
evaluation functions is given in section 7.3. 


77 



6 . 6 Sequence of. Concavity Deletion 

Sequence of concavity deletion becomes extremely important with 
the possibility of cut Interaction, cut modifications and 
approximation interaction, hence forward generically called 
deletion interaction. If there is no likelihood of either of 
deletion interaction, then the order of deletion is immaterial. 

As cuts are restricted to the skeletal interaction zone, cut 
interactions are possible only for the concavities which are 
visible to each other across the interior of the profile (ref. 
section 3, chapterZ). Approximation and cut modifications are 
possible for concavities proximal alone the boundary of the 
profile. Since cuts are based on the skeletal interaction, the 
proximity of the concave skeletal interaction strine reeions 
increases the possibility of cut interaction. Due to the nature of 
the skeleton, proximity of the concavities alone boundary is bound 
to reflect in proximity of the skeletal lines. Thus deletion 
ordering ia highly important for concavities with skeletal lines 
in proximity. 

As noted earlier, a complex profile may have large number of minor 
concavities. The deletion interaction of such minor concavities is 
not likely to produce significantly different decompositions. 
Accordingly, the concavities present in the polygon have been 
given properties of ’peripherality’ and ’importance' (ref. section 
6.2, chapter 2). These two properties of concavities can be used 
in combination to decide the deletion order. 

The peripheral and not important concavities are merely 
disturbances on polygon boundary. These sh^^ldy**^ deleted first to 
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fiive a smoother boundary; otherwise they will unnecessarily 
interfere in the deletion process of important concavities, 
complicatinfi the solution space. 

Among the non-peripheral concavities, the order of deletion of 
’not important’ concavities does not substantially affect the 
final decomposition solutions. The important concavities are 
deleted first and then finally the non-peripheral and unimportant 
concavities . 

The properties of the concavities are evaluated with respect to 
the particular polygon to which they belong (ref. section 5). This 
necessitates an ordering of the constituent polygons. This problem 
is circumvented by taking the global reference values for the 
property evaluation, o.g. for peripherality and importance 
evaluations, the reference length with respect to which the 
’smallness’ or ’largeness’ of radius is to be evaluated, is taken 
with respect ' to the whole skeleton and not the skeleton of the 
particular profile. This ensures that large polygons are processed 
more carefully than the small polygons. 

For the high peripherality and low importance concavities, the 
order is not much important. For non-peripheral concavities, the 
deletion will in general be in a decreasing order of importance. 
When multiple concavities are seen to be equally important, all 
are presented for deletion. The ’usefulness’ of deletions then 
decide which one will be explored before others. If any of these 
concavities are proximal then the decision is mark d as critical. 
The critical decisions are used for generation of multiple 
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decomposition schemes, by explorine all the useful alternatives at 
the critical decision. 

Criticality of Deletion Order 

The deletion order is critical in the high importance phase of the 
decomposition process. In this phase, a deletion decision which 
has multiple candidates is a critical decision. The proximity and 
importance of the candidates for deletion decide the criticality 
of the deletion order decision. 

6 . 7 . Process Control Strategy s_ 

The decomposition process is designed for production of multiple 
decompositions of equal goodness. Heuristic function 'usefulness’ 
is used to guide the process in local decisions. Another heuristic 
function 'criticality' is used to guide the backtracking and to 
produce multiple goal states. 

At each decision node, concavities are short listed for deletion 
on the basis of the ordering heuristics. The criticality of 
deletion order is evaluated. Deletion heuristics are applied to 
each of the candidate concavities. The resulting solution nodes 
are evaluated for usefulness. The criticality contribution by the 
deletion is evaluated. The total criticality of the decision, 
function of both deletion criticality and order criticality, is 
evaluated. The most useful solution is followed. 

In case of a dead end, the decision tree is backtracked to the 
previous critical decision. After a satisfactory decomposition is 
achieved, the most critical decision nodes are < plored further to 
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produce more decompositions of equal goodness. 


7 . Complete Formulation of Heuristic Rules and Procedures j_ 

The proposed heuristic rules and functions carry out the 
quantification of the qualitative concepts discussed so far. The 
heuristics are based on study of some situations. They are quite 
ad hoc in nature as they represent a tentative embodiment of 
qualitative principles guiding the development. Carefully planned 
experimentation, over a wide range of shapes, is needed to refine 
and improve the heuristics. 

7 . 1 Rules for Cut 2 , 

More than one cuts are generated at the deletion of each concave 
vertex . 

1) A cut is generated for all the profile boundary elements 
skeletally interacting with the concave vertex being deleted. 

2) If the element is a concave vertex then the cut is the 
straight line connecting the two vertices. 

3) If the element is a profile boundary edge then generate four 
cuts by connecting the concave vertex with : 

a) centre of the interaction zone. 

b) the extremes of the interaction zone (two cuts). 

c) the minimum distance line between the concave vertex and 
the profile boundary within the interaction zone, if it is not 
already generated in th'e above cuts. 

4) If the cut lies outside the reflex angle then, 

produce cuts over all other interacting profile boundary 
elements. Reject those lying inside the refle angle or on the 

same side of the ref lex angle as the first cut line. 
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5) Usefulness evaluations of cut : 

Usefulness is the 'Goodness' G, modified by concavity reduction 

factor, i.e. 

U = G * C ; 

u 

where G = V * A , 

and C is the concavity reduction factor of tla cut. 
u m 

C = H CW * C ) ; 
u i ui 

’m’ is the number of concave vertices being affected by the 

cut , 

W is the important weightage of the i'th concave vertex, 
i 

C is the concavity reduction of the i’th concave vertex, 
ui 

Note that a cut joining two concave vertices reduces the concavity 

at both the vortices. 

Evaluation of C j, 

ui 

a) For the concave vertex being deleted : 

C la a function of the final angle at the vertex being deleted, 
ui 

This function is plotted in Fig. 7.1. 0 is the exterior angle of 

e 

the profile at the concave vertex V. 0 is the final interior 

f 

angle at the vertex V after the cut i.e. angle between the two 

lines forming the cut. 0 is the compliment of 0 . If the cut lies 

c e 

In the reflex angle of the vertex, then the final angle is zero. 

This factor allows some concavity to remain after the deletion 

action. In the Fig. 7.1, 

0 = TT - 8 ; 

c e 

© = + 0 . 

0 is the limit on 0 above which the C is zero. Value of 8 is 

, f : ; ui d;: ; v; ; 

proposed to be 0 / 2. 

c 
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b) For other concave vertices : A concave vertex other than the 

one being deleted may be affected by a cut. For these vertices, 

taking complement of the exterior angle to be the measure of the 

concavity (0 = TT - 0 ; ref. Fig. 7.1), 

c e 

C = C(© before_cut) - (0 after_cut)) / TT; 

ui c c 

if the concavity is deleted then (0 after_cut) =0. 

c 

7 . 2 Rules for Approximations (ref. Fig. 6. 5. 1.1) 

1) Straightening up : 

Draw the approximation line by connecting the two vertices 
flanking the vertex being deleted. 

2) Cut modification ; 

If any of the edges being eliminated by the approximation is an 
interior edge then the constituent polygon with the interior edge 
is modified to the new edge formed by the approximation. 

3) Approximation interaction : 

The last approximation overrides the previous approximations. 

4) Usefulness : 

The usefulness of an approximation is the goodness value G (ref. 
section 3) of the decomposition. 

7 . 3 Deletion Criticality 

Let C be the concavity being deleted. 

Let D , i = 1,... ,N , be the most useful deletion solutions, 
i • d 

Let HUD = D , 1 <= h < N , be the deletion with the highest 
hi i d 

usefulness value. 

Let D , i /= h , be the difference measure of t» ^ HUD with D . 

i ^ ■ ■ 


fi 



7.3.1 Evaluation of D 


ii 

The difference evaluation ia commutative, i.e. D CD ,D ) 

f i j 

D CD .D ). 

£ j i 

1) If the two deletlbrfa are a cut and an approximation then 

= 1 . 


D 

f 


2 ) Difference weifihtages for types of cut, D : 

W 

Cref. section 6.5.5 for sources of differences and their 
combinations . ) 

Vertex to vertex cut : D = 0.8. 

U 

Vertex to edfie cut : D =0.6. 

W 

This reflects the observation that a vertex to vertex cut is more 

different from other cuts than a vertex to edfie. The numbers 

proposed for D are ad hoc; they are adjusted so as to lie within 
U 

the lower and upper bounds of zero and one respectively. The 
bounds are the same as those normally used in fuzzy membership 
functions. 

Combinations of two cuts may be, 

a) one line and one line, 

b) one line and two lines, 

c) two lines and two lines. 

3) Difference evaluation between two one line cuts : 

1 • ^ 

The situation is illustrated in Fifi. 7. 3. 1.3. The concave vertex V 

is beinfi deleted from the profile. 0 is the interior anfile at the 

i 

vertex, and 0 is the anele between the two one line cuts beinfi 
d , ' ' 

compared for difference evaluation. 

Cut C has difference weifihtage D and C has D . 

1 U1 2 

D = D * D * C0 /0 ). 
f U1 W2 d i 
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D ia bounded, 0 < D <1. 
f f 

4) Difference between a one line cut and a two line cut : 

The case is illustrated in Fifi. 7. 3. 1.4. V is the concave vertex 

being deleted from the profile. The line L and L form the two 

2 3 

line cut cut2 The cut is the one line cut. 

1 

Cut C has difference weightage D , 

1 U1 
Line L of cut C : D , 

2 2 U2 

Line L of cut C : D . 

3 2 U3 

D = ( D * D * (0 /TT) + D * D * (8 / TT) ). 

f W1 W2 12 m U3 13 

D is bounded. Minimum D = 0. 

f f 

For maximum D , 
f 

0 = TT, 

12 

0 =0 (ref. section 7.1, number 5) 

13 d 

with 0 -> 0, 0 -> TT / 2 , 

e 13 

substituting, maximum D = 0.96 < 1. 

f 

5) Difference between two two line' cuts : 

The case ia illustrated in Fig. 7. 3. 1.5. V is the concave vertex 

being deleted from the profile. The line L and L form the two 

13 

line cut cut . The cut is formed by lines L and L . 

1 2 3 4 

Cut C is a two line cut comprised of lines L and L . 

1 12 

Cut C is a two line cut comprised of lines L and L . 

2 3 4 

Line L of cut C has difference weightage D , 

1 1 W1 

Line L of cut C has difference weightage D , 

2 2 U2 

Line L of cut C has difference weightage D , 

3 3 W3 

Line L of cut C has difference weightage 'D . 

4 4 W4 
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Reflex ancle 



C< 7. 3 .1.4 Difference between a One Line Cut and a Two line Cut 


Cut 



Fic. 7. 3. 1.5 Difference between Two Line Cuts. 


9 

1 

U = ( ) * D * D 

f U1 U3 

TT - 9 


e 


9 

2 


TT - 9 


) * D * D . 
U2 U4 


e 

D is bounded. Ifinimum D =0. 

f f 

For maximum D 

f 

9 -> 0, 9 +9 -> TT + 9 (ref. section 7.1, number 5), 

e 12 d 

i.e.9+9->3*'TT/2, 

1 2 

substituting, maximum D = 0.96 < 1. 

f 


7.3.2 Criticality Evaluation for D 

i 

The deletion criticality C (D ), between pairs of the HUD and D 

r i i 

is the fuzzy logical conjunction of importance of the concavity 

Imp(C), difference D , and equality of usefulness of HUD and D 

fi i 

i.e. the fuzzy predicate 

equal (usefulness (HUD) , usefulness (D )). 

i 

C (D ) = Imp(C) AND D AND equal (usefulness (HUD) , usefulness(D ))- 
r i f i . i 

7.3.3 Deletion Criticality Evaluation j. 


It is the fuzzy logical disjunction over all the D . 

i 

C = OR C (D ), i = 1,...,N , i /= h . 
r r i d i 

7.4 Rules for Deletion Order i. 


1) Classification of concavities into four types : 
a) High peripherality and low importance. 
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b) Hifih importance and hifih peripherality . 

c) High importance and low peripherality. 

d) Low importance and low peripherality. 

It ie the application of the fuzzy logical functions; 

type = highCperipherality) AND low(importance) ; 

1 

type “ high( importance) AND high(peripherality) ; 

2 

type = low(peripherality) AND high(importance) ; 

3 

type = low(peripherality) AND low( importance) . 

4 

2) Check existence of concavities of type one by taking an 

alpha cut [13] [19,20] on the set of membership grades of type 
one. The suggested value of the alpha cut is 0.2. If there are no 
type one concavities then process the concavities for type two, 
three and four as given in the following parts 3 and 4; else 
select the concavity with highest type one membership. 


3) Shortlist the important concavities from the profile, 

IC , i = 1,2, .. .N . 

1 ic 

4) Select the most important concavity (concavities) from them. 

All of the most important concavities are presented for deletion. 

MIC , j = 1,2, . . . ,N 
j mic 

5) Evaluate the proximity of the most important concavity with 

the important and the most important concavities. Deletion order 

criticality of a concavity pair is the fuzzy logical conjunction 

of proximity and importance; i.e. 

C CC ,C ) = P (C ,C ) AND IrapCC ,C ). 
r i J r i j i J 

The order criticality of the deletion order decision is the fuzzy 

disjunction of all the pairs; i.e. 

C = OR C (C ,C ). 

: r r . 1 j 
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The concavities for which the order of deletion is critical, are 
noted, alona with the criticality. The order of deletion is taken 
to be the descendine order of importance of the concavities. The 
deletion order criticality la utilised in deciding upon the 
backtracking (ref. section 5). 

7 . 5 Rules for Process Control 

1) Criticality of a decision : 

It is the fuzzy logical disjunction of the deletion criticality 

C and the order criticality C , C = C OR C . 
rd ro r rd ro 

2) Judgment on solution state : 

Dead end recognition : It is the fuzzy predicate ’low' applied to 
the local goodness value; i.e. lowCgoodness (HUD) ) . 

Goal state recognition ; The necessary condition is that all the 
concavities are deleted. It is the fuzzy logical ’satisfactory’ 
applied to the goodness of the final decomposition solution, 
sat la factory (goodness (Solution) ) . 

3) Critical decision recognition and backtracking ; 

Critical decision recognition is the application of the fuzzy 
predicate ’high’ to the decision criticality list. Backtracking is 
done upto the last ’high’ criticality decision. 

8. Complete Algorithm j. 

1) Evaluation of the skeleton and its properties. 

2) Generation of Concavity deletion solutions ; 

a) Deletion order analysis to determine Concavities for 
deletion. 

b) Deletion of concavity. 

c) Usefulness evaluation of each solution- 
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d) Criticality of the decision. 

The roost useful explored first. 

The decision stored on stack. 

3) Goal state / dead end / end of process recognition. 

4) If none- of the three go back to skeleton evaluation. 

5) if goal state, 

a) Determine and- analyse the bays formed by 
approximations . 

b) store the final results, 

c) Get the list of most critical decisions, 

d) backtrack to the last of them, 

e) go to skeleton evaluation. 

6) If dead end, 

a) Get the list of roost critical decisions, 

b) backtrack to the last of them, 

c) go to skeleton evaluation. 

7) if it is the end of the process then quit. 

9 . Example of Decomposition 

The decomposition algorithm is illustrated by an example. Consider 

the profile given in the Fig. 9.1. 

The concavities are V , V , V . 

4 7 9 

Importance order (descending) of concavities (ref. section 6.2.3, 
chapter 2) : 

V and V have approximately same importance. V has much lower 

7 9 ^ 

importance. 

V > V > V . 

7 9 4 

Per ipheral ity of all the concavities is low (ref. section 6.2.2, 
chaptor 2). 

91, ' 



Candidates for deletion are V and V . 

7 9 

V and V have good proxiinality (ref. section 6.2.4, chapter 2); 

7 9 

hence the order of deletion of V and V is critical. 

7 9 

1) Deletion of V in First Step : 

7 

V interacts with edge E and vertex V . Hence two cuts are 

7 4 9 

possible; (Note that the three cuts between V and E (ref. 

7 4 

section 7.1) are reduced to one cut at the centre of the zone of 

interaction.) Cut between V -E and Cut V -V . Approximations have 

7 4 7 9 

low usefulness and hence are not considered. 

Usefulness order (ref. section 7.1) ; 

The usefulness of the two cuts is almost the same. But, 

Cut V -V > Cut V -E . 

7 9 7 4 

It may be noted that the cut V -V does not eliminate the 

7 9 

concavity at V ; but it only reduces the concavity. 

9 

Hence V may be deleted by the cut V -V or by cut V -E . 

7 7 9 7 4 

2) Deletion of V by cut V -V : 

7 7 9 

Now the profile is decomposed in to two profiles, as shown in Fig. 

9.2. The subprofile P , has concavities V , V , V . All have low 

2 4 9 7 

importance and low peripherality . The order is also not critical. 

Hence let the order be V -V -V . 

9 7 4 

For V deletion, Cut V -V is the most useful (ref. Fig. 9.2). 

9 9 4 

Now only V is the concavity. The approximation V -V is the most 
7 9 6 

useful for deletion of V . The resulting final decomposition is 

7 

shown in Fig. 9.3. 
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7 7 4 

The profile is decomposed in to two subprofiles as shown in Fig. 

9.4. The subprofile P has concavities V , V . Both the 

2 9 4 

concavities have similar low importance and low per ipheral ity . The 

order of deletion is not critical. Let V be deleted first. Of the 

9 

possible cuts and approximations, cut V -V has highest 

9 4 

usefulness, outstripping the remaining alternatives by a large 

margin. Hence cut V -V is taken. The resulting decomposition is 

9 4 

shown in Fig. 9.5. 

4) Deletion of V in First Step : 

9 

As the deletion order of V and V is critical, the solution is 

7 9 

backtracked to the first step and V is deleted in the first step. 

9 

V interacts with edge E , vertices V and V . The possible cuts 

9 7 4 7 

are V -E , V -V , V -V . Cut V -V is the same as the deletion of 
9 7 9 4 9 7 9 7 

V by cut V -V , in the first step. It will lead to the same 

7 ' 7 9 

solutions. 

Usefulness order : Cut V -V > Cut V -E . The usefulness of cut 

9 4 9 7 

V -V is appreciably larger than that of cut V -E as the cut V - 

9 4 • 9 7 9 

V deletes concavities at both V and V . 

4 9 4 

Hence V is deleted by cut V -V . The profile is now decomposed in 
9 9 4 

to two subprofiles, P and P (ref. Fig. 9.6). Subprofile P has 

1 2 2 
concavity V . The concavity interacts with edges E , and E . The 
7 9-4 4 

cut V -E produces a decomposition which has been produced 
7 4 

earlier. Hence the concavity V is deleted by the cut V -E . The 

7 7 9-4 

resulting decomposition is shown in Fig. 9.7. 

Thus the given profile is decomposed in to three different 
decomposition schemes, given Fig.s 9.3, 9.5, 9.7. In this case, 
the three decompositions are not much different from one another. 




Chapter 4 


CONVEX POLYGON DESCRIPTION 

1 . Introduction 

Convex polygon description is a stage in the shape description of 
a profile. After the profile is decomposed into a set of 
constituent convex polygons, a description of the constituent 
polygons is obtained based on its qualitative aspects. In case of 
a convex profile this description itself becomes the shape 
description. For profiles with concavities, description of 
constituent polygons is useful in projection matching during the 
allocation. 

The' description of a convex polygon can not be only qualitative. 
It must combine both, the qualitative as well as the precise 
quantitative geometric features of the polygon. Polygon 
description is a complex, structured geometrical knowledge. An 
approximate global description neglecting minor features, may be 
used as a criterion to pare down and delimit the feasible solution 
space for allocation. Then the exact description is used for final 
allocation in the reduced solution space. 

1.1 Objective of Convex Polygon Description j. 

A polygon description should be independent of translation, 
rotation and mirror imaging transformations. It should identify 
features of the polygon which are likely to be helpful in nesting. 
It should also be capable of representing the geometrical 
relations between the features. The features and their geometrical 
interrelations together constitute a complete polygon description. 
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1.2 Shape and Size of a Polygon j_ 

A polygon is characterized by its shape and size. Shape of a 
polygon may be specified by the shape features and their 
interrelations. Shape is a necessarily qualitative property. Size 
is a quantitative property. Size can not be defined as a single 
global quantity. Size specification must be done with reference to 
the shape. Thus size and shape are intricately and inextricably 
woven together. They together define the polygon. We propose to 
define size of a polygon as the size of its shape features. The 
size of the dominant features and their combinations will be used 
as the global polygon size. 

1 . 3 Basic Approach to Description 

Convex polygon description is based on classification of the 
polygons into a set of well defined shape classes. The proposed 
shape classification is hierarchical in nature . It can be modeled 
as a decision tree based on the features detected in the convex 
polygon. The proposed shape classes completely cover the whole 
shape space. They are based on shape features and their 
geometrical interrelations. The features and interrelations form 
the basis of the shape space. They span the whole shape space. The 
coordinate directions of the shape space are divided into 
different regions. The shape space is classified into different 
groups formed by combinations of these regions. The division and 
the subsequent classification are based on qualitative rules . The 
rules are in turn based on a systematic shape space exploration 
carried out by studying features and their interrelation regions 
and their combinations. 

All the divisions and the classifications are fuzzy in nature. 
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This preserves the continuity of shape variation from one class to 
another by allotting grades of membership of classes to the 
polygons. A polygon may have nonzero membership grades for many 
classes. Thus the complex relationships among classes of shapes 
can be modeled in this approach, which may not be possible in a 
crisp approach. 

This approach is best illustrated by an example. The polygon in 
Fig. 1.3.1 could be described as just an approximate regular 
hexagon or a regular hexagon with a projection. Either of the two 
would be an approximate description. Neither of them alone would 
be sufficiently accurate. In the present work, both the 
descriptions are maintained; and the polygon is given different 
grades of membership or different confidence values for both the 
descriptions. For example, it may be described as a regular 
hexagon with say, 0.7 confidence and a regular hexagon with a 
projection with 0.75 confidence value. 

An extremely simplified illustration of a general case is given in 

Fig. 1.3.2. The feature space is just two dimensional i.e. only 

two features are considered. X , X ,...,X are the zones of the 

12 m 

feature ’X’ along x axis. Y , Y ,...,Y are the zones of the 

1 2 n 

features of the feature ’Y' along the y axis. Combinations of 

these zones give us classes C,C,... . A profile P lying 

12 1 

halfway between C and C has equal membership grades for C and 

2 5 2 

C and zero membership grades for the other classes. The profile 
5 

P2 may have same membership grades for C and G , but it has 

2 5 

nonzero membership grades for C and C also. The whole set of 

1 ' 4 

class descriptions and the membership grades together describe the 
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profiles. 

If the description process is modeled as a decision tree, at each 
decision node some alternatives may be ruled out altogether. All 
the remaining alternatives need to be explored. The decision 
process spawns one or more subprocesses at each of the nodes. In 
each subprocess, the decision taken at the parent node to spawn 
the process is firm. In Fig. 1.3.1, at a decision as to’ whether 
the profile is to be a hexagon or a hexagon with a projection, 
both the decisions will be taken, each one spawning its own 
subprocess. In the subprocess* spawned by the decision of hexagon 
with a projection, in all further analysis, the profile will be 
treated as a hexagon with a projection and similarly for decision 
of the profile being just a hexagon. If the projection were to be 
very minor then the alternative of hexagon with a projection 
would be neglected altogether. Only when the profile is such that 
it lies somewhere halfway in between two classes, both the classes 
would be considered. This type of analytical process will be 
referred to as a fuzzy spawning process in further work. 

1 . 4 Shape Features 

A convex polygon may be deemed to be formed by the process of 
distortion of a regular convex polygon with the same number of 
sides. The distortions take the form of different stretching 
operations on the regular polygon. The stretching operations may 
result in formation of long sharp projections, mild projections, 
or elongations of the polygon. The stretching operations on a 
regular polygon are illustrated in Fig. 1.4.1. A regular polygon 
has only one fork in its skeleton. All the ACS I Ls converge to that 
fork. NACSILs are absent. The fork and the skeletal circle of the 
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^2. Xg ... Xm ^ 


Fig. 1.3.2 An Illustration of a Feature Based Classification 
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regular polygon are drawn in the figure. Only the sides of the 
polygon undergoing the stretch operation are shown to avoid 
cluttering. 

In Figure 1.4.1a, vertex V is pulled out pushing back vertices V 

2 1 
and V along the skeletal circle. The length of sides V V and 

3 • 12 

V V increases at the expense of the other sides of the polygon. 
2 3 

The angle 4> is reduced to a small value. The symmetry of the 

skeleton of the regular polygon is disturbed. The angles and the 

lengths of the skeletal lines are no longer equal. A long sharp 

projection, a beak, emanates at the vertex V . In Figure 1.4.1b, a 

2 

similar but weaker distortion results in a mild projection, a 
bulge. This type of stretch operations is referred to as type one 
stretch operation henceforth. 

In Figure 1.4.1c, the edge e is pulled out maintaining its 

2 

orientation. The interior angles of the regular polygon remain the 

same. But an extra skeletal line is formed due to interaction 

between edges E and E . It results in an elongation of the 

13 

polygon. It also generates additional forks and the NACSILs in the 
skeleton. This type of stretch operations is referred to as type 
two stretch operation henceforth. 

The two stretch operations mentioned above can produce, on 
repeated applications in different order combinations, all the 
possible variations in the skeleton of a convex polygon. Hence 
they can span the whole shape space. 

The resultants of stretching operations, i.e. beaks, bulges and 
elongations, are important shape features. Very weak stretch 
operations result in minor distortions. Minor projections and 
minor elongations do not greatly affect the shape. They can be 
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approximated out to simplify the polygon and its skeleton. On the 
other hand, if the stretch operation is very strong, then the 
other features pale into comparative insignificance. The stretch 
operations dominate the shape. It is the middle level resultants 
of the stretch operations, beaks, bulges and elongations as listed 
above, are most important shape features. These important shape 
features and their interrelations form the basis of the shape 
space for explorations. 

Of the shape features, beaks and bulges are the result of the same 
type of stretching operations, whereas elongation represents an 
entirely different operation altogether. It is contended that the 

A 

sequence of application of these two operations can be separated 
into two distinct sequences. This contention is more rigorously 
justified later on in this chapter, in section 2.6.12. These two 
operations and their sequences can be explored separately. They 
are classified separately and finally combinations of these 
classes are used to classify the convex polygons. 

For regular polygons, as the number of sides increases, the 
polygon becomes more and more similar to a circle. The difference 
in shape of a polygon, produced by addition of an extra side, is 
negligibly small if the number of sides is large. Regular polygons 
with five or more sides are treated as circles. This introduces 
extra approximation. But it is necessary to curtail the complexity 
of the convex polygon classification problem. Triangle and 
quadrilaterals are treated as such. 

1.5 Polygon Simplification t_ 

The minor features need to be approximated out of the polygon to 
airoplify it and thereby make it amenable for analysis. 
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Approximations are done as a prelude to any attempt at shape 
analysis. A long sequence of minor stretch operations, applied in 
different directions, may distort the shape in vastly different 
ways. A thorough analysis of the effect of such sequences would be 
quite complicated. It is not attempted in this work. Long 
sequences of minor stretch operations are approximated by short 
sequences of major stretch operations. 

Minor type one stretch operations give rise to irregular polygons 

which are very similar to the regular polygons (ref. Fig. 1.5.1). 

In Fig. 1.5.1, polygon in (a) has been distorted to (b), by an 

operation near vertex V . Vertices V and V are pulled towards 

6 15 

vertex V . In the process, V and V are also affected. V and V 
6 2 4 1 6 

maintain their positions. The ACSILs of the skeleton are shifted 

from their positions and the angular symmetry is disturbed. It can 

be seen that the shape of the polygon is hardly affected. The 

resultants of such operations are neglected altogether in 

classification. Minor type one stretch operations, when coupled 

with some major stretch operations, or just application of some 

major stretch operations may give rise to rounded corners (ref. 

Fig. 1.5.2). In Fig. 1.5.2, two major type one stretch operations 

have produced two beaks. These two operations have generated a 

rounded corner V -V -V . 

12 3 

Minor type two stretch operations produce small elongations. Again 
these stretch operations hardly affect the shape of the polygon 
(ref. Fig. 1 . 5 . 3 ). Extremely strong type two operations are 
illustrated in Fig. 1.5.4. It can be easily seen that they 
generate rounded corners. A sequence of minor type two operations 
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Fig. 1.5.1 Minor type 1 Stretch Operation 



Fig. 1.5.2 Two Major Type 1 Stretch Operations. 





Minor elongation / 
type 2 stretch 
operation 


Fig. 1.5.3 Minor Type 2 Stretch Operation 



(ref. Fig. 1.5.5) generates a boundary pattern which closely 
resembles a curved boundary. 

2 . Skeletal Basis for Description Process 

The skeleton and skeletal properties provide a very good basis for 
detection and replacement of features to be approximated. 
Identification of shape features and their interrelations is 
facilitated by skeletal characteristics. 

2 . 1 Description of Convex Polygon Skeleton j_ 

The skeleton of a convex polygon is devoid of any concave skeletal 
interaction. It is composed of only convex skeletal interaction 
lines, i.e. CSIL as they are referred to in the chapter 2. The 
skeleton is again a geometrical line diagram. It is a geometrical 
tree. The nodes of the tree are the forks, and the end points of 
the adjacent edge convex skeletal interaction lines i.e., the 
ACSILs. The edges are the CSILs. The non-ad jacent convex skeletal 
interaction lines i.e. the NACSILs are the edges between the forks 
and the ACSILs are the edges between forks and the ACSIL 
endpoints. 

2 . 2 Observations on the Properties of Convex Polygon Skeleton ± 

The skeletal circle radius(r) as a function of the distance(s) 
along the skeletal lines has some important properties. 

2.2.1 The Uniroodality of r(8) i_ 

Along a path in the tree from a leaf node to another leaf node, 
rCs) is strictly unimodal ; it is increasing - decreasing. It is 
not concave upwards . In Fig. 2.2.1, as the skeletal path from V 

1 

to V is traversed, the r increases from zero till it reaches its 
, 5 

maximum at fork f and starts decreasing again till it reaches 

1 

zero again at V . 

5 


106 



-rounded corner 



operation 

Fig. 1.5.4 Extremely Strong Type 2 Stretch Operation 
Curved boundary 



Fic. 2.2.1 Uniraodality of r(») 



Along any other path, r(3) is either monotone increasing, monotone 
decreasing, or unimodal and convex upwards. 

2.2.2 The Convexity of r(el ^ 

The dr/ds is also unimodal, convex upwards along a path from a 

leaf node to leaf node. The area underneath the rCs) distribution 

is convex. In Fig. 2.2.2, any edges of the convex polygon, lying 

to the right of the V -V line, must necessarily lie within the 

2 4 

shaded region. The angles between the edges must be larger than 

the angle (|> between e and e . As the dr/ds is given by 

12 34 

dr 1 1 

:= ( ) = ( ) sin(e/2) ; 

ds ds/dr quench velocity 

(ref. section 4.1, chapter 2, for quench velocity formulation.) 

as the path is traversed to the right of C , i.e. s increases, the 

1 

magnitude of dr/ds will increase as (0/2) > (<|>/2). But as r is 
decreasing with increase in s, dr/ds < 0 ; and hence as s 
increases dr/ds decreases. This is the situation to the right of 
the location of the maximum r. To the left of maximum r, the 
situation is reversed, and the dr/ds increases with increase in s. 
Thus the dr/ds is also unimodal. This implies that the area under 
the r(s) curve is convex. 

Along any other path, dr/ds is either monotone increasing, 
monotone decreasing, or unimodal, convex upwards. 

2.2.3 Continuity of NACSILs i. 

The NACSILs in the skeletal tree are continuous. They are not 
fragmented into physically disparate units. If only the NACSILs 
are considered neglecting the ACSILs, then the resulting line 
diagram is also a tree. In Fig. 2.2.3, it can be seen that all the 
ACSILs start from a fork and terminate in a vertex, or a leaf node 
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of the skeletal tree. An NACSIL could not exist at the vertex end 
of an ACSIL. It can occur only at the fork end, where it would be 
connected to the other NACSILs at the fork. 

2.3 Fork - NACSIL Tree ^ 

If the ACSILs converging at a fork and the fork are considered to 
be one unit and the NACSILs are considered to be the different 
entities, the effects of the two types of stretch operations can 
be discerned separately. The result is again a tree, with the 
units of fork with its ACSILs as the nodes and the NACSILs as the 
edges. This proposition is based on observation given in section 
2.2.3. This tree has some important properties useful for convex 
polygon classification. The skeletal tree forks with their ACSILs, 
are the nodes of this tree. They are henceforth referred to as the 
forks in this chapter. The forks of the skeletal tree and these 
forks are almost the same; the only difference being the ACSILs. 
Hence the same terminology is retained, rather than introducing 
new terms. 

The tree is referred to as the fork - NACSIL tree. A fork - NACSIL 

tree is illustrated in Fig. 2.3.1. Fig. 2.3.1(a) shows the 

skeletal tree of a convex polygon. The fork - NACSIL tree of the 

polygon is given in Fig. 2.3.1(b). The ACSILs of fork fl in part 

(a) of the figure, i.e. ACSILs f -v , f -v etc., are merged into 

1112 

fork f in part (b). 

1 

2 . 4 Properties of the Fork 2 . NACSIL Tree j. 

In the fork - NACSIL tree, the fork size and proximity definitions 


remains 

the same as 

in 

the skeletal tree. 

The 

intensity 

of 

the 

fork is 

defined as 

the 

sum of intensities 

of 

the ACSILs 

at 

the 


fork . 
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V±o 



(a) Skeletal tree of polygon 





(b) Fork - NACSIL tree of the polygon in (a) 
Fig. 2.3.1 Skeleton and Fork - NACSIL Tree 
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A path in the fork - NACSIL constitutes a string of forks. A 
string of forks or, a fork string, may be geometrically long or 
geometrically short, depending on the physical distances in the 
string. The physical distances are the lengths of the NACSILs. A 
fork string may be called long or short depending on the number of 
forks in it. In a dense fork string, the neighbouring forks are 
proximal to each other. In a sparse fork string neighbouring forks 
have a low proximality i.e. they are far flung apart. 

2 . 5 Symbols Used in Representation of Fork - NACSIL Tree ± 

The following set of symbols is used to represent the fork - 
NACSIL tree and its various subtrees. 

1) A locally large fork is represented by LL. 2) A globally large 
fork is represented by GL. 3) A locally small fork is represented 
by LS. 4) A globally small fork is represented by GS. 5) Proximity 
between a pair of forks is denoted by a short line joining the two 
forks. Distant pair is denoted by a long line joining the two 
forks . 

For example, 

LL LS represents a distant or not proximal pair of locally 

large forks . 

GL — GL represents a proximal pair of globally large forks. 

The same set of symbols may be used with the skeletal tree without 
any conflicts as the two trees are essentially the same; they are 
different representations of the same tree. 

2.6 Oba er vat i ono on the Properties of Fork ^ NACS I L Tree ; 

The observations on the fork - NACSIL tree are equally valid on 
the skeletal tree. All the observations essentially follow from 
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the observations in section 2.2. The observations form the basis 
of the approximation and classification process. Detection of 
situations for approximation, the actual approximations and the 
ordering of the approximation process is done on the basis of 
observations from sections 2.6.1 to 2.6.10. The classification and 
description is based on the rest of the observations. 

The observations are in the form of logic propositional rules, 
i.e. 

if 'A’ then ’B’ 

where A is a feature of the fork - NACSIL tree or skeletal tree, 
and B is a shape or boundary feature. The tree features are 
qualitative fuzzy predicates such as 'a string of locally large 
forks’. As a result, the feature detected are not mutually 
exclusive. The same tree feature may satisfy more than one 
qualitative predicates with different confidence values. 

The approximation and classification process is by and large a 
fuzzy spawning process. In case of border line cases of the 
qualitative predicate values, all the alternatives ensuing thereof 
are followed up and expanded. For example, for a fork pair which 
la ambiguous predicate values of say, around 0.5 and 0.6 for all 

of LL--LL and LL- LL and LL LS, all the three solutions 

will be followed. Some of these solutions may in later analysis 
give very low confidence values for all decisions thereof. These 
will be discarded; but the rest will be followed. In the final 
description, all the descriptions of not very low confidence 
values, spawned by these three predicates will be present. This 
strategy is in conformance with the classification scheme 
expounded in section 1.3. It ensures that the approximations will 
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not be too damaging to the description. 

The features of the fork - NACSIL tree and the skeleton tree 
analysed in the following observations completely cover all the 
possible ranges of the features of the tree. The fuzzy approach to 
the feature categories ensures completeness. 

2.6.1 Location of Small Forks 

Small forks are always located at the’ extremes of the strings. 
This observation follows directly from sections 2.2.1 and 2.2.2.' 
For a string from a leaf node to leaf node, the larger forks are 
in the middle of the string. The large forks dominate the shape of 
the polygon. The small forks can be approximated out by a pair of 
polygon edges and the ACSIL between them. 

It can be seen from Fig, 2.6.1, that the small forks must lie at 

the extremes, else the convexity of the polygon will be violated. 

In the figure, dur to convexity of rCs). r < r , and 

1 2 

r > r > r . 

2 3 4 

2.6.2 Relation between Size. Proximity and Angle between 
Interacting Edges j_ 


In the Fig. 2.6.2, 

edges e 

12 

between vertices V ad V , 

1 2 

and 

e 

34 

between vertices 

V and V , 

3 4 

give rise to a fork f at 

1 

point 

0. 

Their interaction 

continues 

till point 0’ where the 

edge 

e 

23 

between vertices V 

2 

and V stops their interaction, giving rise 
3 

to 


another fork f at O’. The radii of the skeletal circles at forks 
2 

f and f are R and r respectively. The distance between the forks 
1 2 

is d . In the figure, 

KOQ) = R / sin(4>); 

1(PQ) = r / sinC4*); 
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By similarity of the triangles QOA and QPB, 


(R / sin((j>)) - d r 

CR / sinC(|>)) R 

i . e . 

R - d * sin(<|>) = r ; 


Proximity ratio P is defined asP =d/R, 0<P < infinity ; 

r r r 

Largeness ratio L is defined as L = r / R, 0 < L <= 1 ; 

^ r r 

As <j> is half the edge angle, 0 < (j> < 90 ; 

P * sin(<j>) + L =1. ( 2.6.2.) 

r r 


The proximity ratio and the largeness ratio are in conformance 
with the definitions of proximity and size given in the chapter 2. 
Hence the relation 2.6.2 can be used as an indicative of the 
relations between size, proximity and quench velocity. 

2.6.3 Uneven Size Proximal Fork paira i. 

In a pair of forks, for a locally large fork to be proximal to a 

locally small fork, the angle ()» must be quite large. For a locally 

large and small pair of forks, L is low. For a low L , P * 

r r r 

sinC(l>) must be high, from relation 2.6.2. For a proximal pair, P 

r 

should be low. Note that P lies between 0 and infinity. Low P 

r r 

does not mean that it should be low compared to one. A low P may 

r 

be close to one (ref. section 7.2, definitions of the fuzzy 

qualitative predicates). ^ should also be quite close to one, to 

satisfy the relation 2.6.2. Hence the angle should be comparable 

to 90. Such a situation is illustrated in Fig. 2.6.3. It may be 

treated as a rounded corner. The edges e and e may be extended 

1 2 

till they meet, to approximate and simplify the polygon, as the 
area approximations involved will be quite small (ref. section 7, 
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for the exact area approximations). 

2 ■ 6.4 Uneyen Size Distant Fork Pairs 

In a LL LS fork pair , the small fork can be nefflected without 

any appreciable loss of accuracy. 

The locally small forks indicate rounded corners and they can be 

approximated out by a pair of edges and an ACSIL between the 

edges. In the Fig. 2.6.4, the large fork f is a major feature of 

1 

the shape. The small fork f lies very close to the tip, and is 

2 

the result of a cut off tip. 

2.6.5 Proximal Fork Pairs 

A proximal fork pair is the result of a minor type two stretch 
operation. It must necessarily be LL--LL. The proposition in 
section 2.6.3 takes care of the LL — LS fork pair. The two proximal 
forks can be treated as one fork without appreciable loss of 
accuracy. The two forks can be merged into one fork. All the 
skeletal lines are preserved, but they, instead of emanating from 
the two original fork points, emanate from one fork point lying 
somewhere halfway between them. The basic angle properties of the 
skeletal lines such as ’skeletal lines bisect angle between the 
edges’ and so on, remain approximately true. 

In Fig. 2.6.5, the two forks f and f are large and are very 

1 2 

close to each other. The overall shape of the profile is not 

disturbed much by the minor type two stretch operation. It 

resembles a regular hexagon. 

2.6.6 Global Size of Forks i_ 

Globally small forks are indicators of very strong type two 

stretch operations. The stretch operation is so extreme that the 
edges stretched out become very small in comparison with the rest 
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small fork 


= large fork, F = 

2 

2,6.3 Uneven Size Proximal Forks 



= large fork, F = small fork 

2 

2.6.4 Uneven Size Distant Forks 



F and F are comparable. 

1 2 

ig. 2.6.5 Proximal Fork Pair 
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of the polygon. They represent rounded corners in the polygon and 
hence can be approximated as the rounded corners detected in the 
earlier observations. This observation is in line with the one 
given in the section 2.6.4. 

2.6.7 Subtree of Forks 

It may be noted that the propositions outlined in sections 2.6.1 

to 2.6.6 may be, used to simplify fork pairs. As a result, only LL- 

LL fork pairs are left in the skeleton. Concatenation of 

these pairs may give rise to fork string or fork trees. 

LL LL pairs have large values of the proximity ratioCP ) and 

r 

largeness ratioCL ). It can be inferred from relation 2.6.2 that <J> 

r 

should be quite small. The fuzzy qualitative predicates are so 

defined that the LL LL pair exists only at (|> below 50 (ref. 

section 2.7). The number of CSILs with low emanating from a fork 
is limited. It can be seen from Fig. 2. 6. 7.1 that the angles 0’ 
must sum upto 360. Hence only three CSILs can maximally be present 

in a fork with 4> of about 60. As LL LL fork pairs can not 

exist at <J> equal to 60, a subtree of large trees can not exist in 
skeleton. 

\ 

A subtree of large forks can occur only if they are proximal. 

Subtrees of forks can be reduced to one fork by repeated 

application of proposition in section 2.6.5. A subtree of large 

forks, f , f , f , and f , is illustrated in Fig. 2. 6. 7. 2. 

12 3 4 

A subtree of large forks can not have distant forks. Thus, if a 
compact subtree is taken care of by approximations, the NACSILs 
resulting from type two stretch operations will be in the form of 
a set of strings rather than forming an arbitrary tree. This 
argument is again similar to that in sections 2 . 6 . 4 and 2 . 6 . 6 . 
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2.6.8 Curved Boundary 

A numerically short dense string indicates the presence of a small 
curved boundary segment which can be treated as a rounded corner. 
The definition of the numerical shortness is that the string 
should have less than five forks. The rounded corner may be 
approximated by extension. 

A numerically long dense string indicates the presence of a curved 
boundary. In Fig. 2.6.8.1(a), the curved boundary from vertex V 

2 

till vertex V ' produces a skeletal interaction wherein a large 
2 

number of locally large forks are densely packed in a small length 
of the skeleton. This is the result of a type two stretch 
operation on an approximately circular boundary. 

The ACSILs of the forks are of very low intensity, as the angles 
between them are close to 180. A curved boundary on both the sides 
of the skeleton produces closely packed low intensity ACSILs, 
emanating from almost every fork, on both the sides. For a curve 
only on one side, such ACSILs are restricted only on the side of 
the curve. Both the cases are illustrated in Fig. 2. 6. 8.1. 
Replacement of a string will be complete when all the curves it 
generates are approximated. 

A dense long string of large forks may be a part of larger 

substrings (ref. Fig. 2. 6. 8. 2). The substring f -f is part of 

1 2 

strings V -V , V -V and so on. In such a case, it is necessary to 
13 1 4 

identify the curve limits on the boundary and the string which 

gives rise to the curve. Ue propose to identify the string and the 

curve by determining the minimum of the areas enclosed by the 

string and the curved boundary segment associated with it. In the 

Fig.2.6.8.2, the curve V -V will be associated with the upper 

2 3 
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Sparse strings interspersed with dense strings are reduced to only 
sparse string by proposition in section 2.6.8. 

After all the approximations outlined in sections 2.6.1 to 2.6.8 
are carried out, the dense strings, all the small forks and hence 
the small fork strings are deleted from the fork - NACSIL tree. 
The tree is composed of only sparse strings. Since the sparse 
strings can not be branched, more than one sparse strings can not 
exist in the tree. The fork - NACSIL tree is composed of only one 
sparse string. 

2 . 6 .10 Order o_f Application of Rules Formulated from Observations 

m 

Repeated applications of the propositions stated in sections 2.6.5 
and 2.6.6 to a dense string can produce quite erroneous results. 
Hence dense strings need to be detected and replaced before any 
other approximations are applied. 

2.6.11 The ACS I La in a Fork 

A sequence of large number of very low intensity ACSILs within a 
fork represents a circular curve in the boundary. 

If a circular curve is delimited on either sides by an ACSIL 

showing a sharp increase in intensity and the angle between the 

extremities of the curve is low, then it is a rounded corner. In 

Fig. 2.6.11, the boundary segment V -V ’ is a circular curve and 

1 1 

V -V ’ is a rounded corner. 

4 4 

2.6.12 Separation of the Elongationa from the Polygon i_ 

Type two stretch operations preserve the angular relations between 
ACSILs and the polygon edges. Their effect is restricted to 
generation of additional forks and NACSILs joining the forks. The 
change in lengths of the polygon in these operations is a result 
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of these changes. It is not the fundamental distortion in shape; 
rather it is a secondary distortion. 

The type one stretch operations basically distort the angular 
relations. The length changes can be treated as the results of the 
angular distortions. For a given fork radius, the angle 
specifications are sufficient to completely reproduce the profile. 
Hence the type one stretch operations can be treated as basically 
angular distortions. 

A separation of elongations by merging all the forks together, 
preserves all the angular relations. Feature analysis for beak and 
bulge features can be carried out on the merged fork as well, as 
it could be on the unseparated tree. These two feature analyses 
would be exactly equivalent. The separation of elongations does 
not affect the feature analysis. 

separation thus cleanly sifts two distinctly different phenomena. 
It helps in simplifying the analysis preceding the classification 
and description of convex polygons. 

By means of separation, the tree and hence the polygon is is 
divided into a one fork polygon and the elongation line diagram. 
These two are explored and analysed separately. The subtrees of 
the fork - NACSIL tree, are merged into the single fork. They are 
thus taken care of. 

The separation is illustrated by two examples CFie^res 2.6.12.1 
and 2.6.12.2). In Fig. 2.6.12.1, the original polygon and the 
skeleton are given in part (a). For simplicity, it does not 
contain small forks and dense strings. The fork - NACSIL tree is 
separated into unifork C^J) snd elongation (c). Fig. 2.6.12.2 runs 
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(a) Original polygon and 

its fork - NACSIL tree 



■ NS-^. 



Elongation 



Unifork polygon and its skeleton 
(b) Seperation 

Fig. 2.6.12. 1 Separation of Uni fork and Elongation - I 
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Elongation 



^7 V6 


Uni fork polygon and its skeleton^^^^ 

(b) Separation 

2 . 6 . 12 . 2 Separation of Uni fork and Elongation 


L2 



along similar lines. Its elongation is along only one direction 
instead of two directions as in section 2.6.12.1. 

2.6.13 Unifork j_ 

The single fork into which the forks of the fork - NACSIL tree are 
merged, and the polygon it represents, is referred to as the 
unifork in this chapter. The exact meaning will be clear from the 
context of usage. 

As observed above, all the angular relations of the polygon are 
preserved in the unifork. For a given radius of the unifork, the 
angles uniquely determine the edge locations. This implies that 
the beaks and bulge features of the original polygon are 
completely preserved in the unifork. 

The unifork analysis and classification requires extensive work 
and is expounded in section 3. 

2.6.14 Elongations 

As observed earlier, the elongations i.e. the NACSILs are 
continuous. The line diagram of elongations is a tree. 
Observations given in sections 2.6.1 through to 2.6.9 and their 
resulting approximations imply that the elongation tree will not 
be excessively branched. After the approximations only one sparse 
string is left in the fork - NACSIL tree. The approximation 
process being a fuzzy spawning process, fork subtrees will either 
give rise to rounded corner approximations or they will be merged 
together as one fork before separation. Due to this processing, 
the elongation tree is extremely simplified. It has only the major 
elongation corresponding to the sparse string, which txy 
observation 2.6.9, can not be branched. Hence elongation tree is 
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only a string of NACSILs. 

Elongations are illustrated in Fig. 2.6.12.1 and 2.6.12.2. At the 
center, corresponding to the unifork location, the elongation 
string has at the most two major branches as evident in the 
figures. 

The angles between successive NACSILs in an elongation strings, 

are very large, almost close to 180. It can be observed in Fig.s 

2.6.12.1 and 2.6.12.2. This observation relies on the propositions 

given above. The justification for the proposition is as follows. 

In Fig. 2.6.14.1, edges e and e form an NACSIL P -0. It 

1 2 1 

terminates at the skeletal circle centered at 0, of radius R, in 

fork f . The next NACSIL in the elongation string 0-P , starts 
1 2 
from 0. The edges e and e form the NACSIL 0-P . 0-P ends in 

3 4 2 2 

fork f at P . For the sake of simplicity, the fork f , its ACSILs 
2 2 2 
and the edge between e and e are not shown. The edges e and e 

3 4 3 4 

are extended till they meet at Q . e is parallel to e . e ’ is 

1 1 2 3 

the extension of e ’ till it meets the edge e at Q . e is the 

1 4 2 3” 

extension of e till it meets e at Q . P ’-0 is the NACSIL that 

1 4 3 1 

would be produced by edges e ’ and e . If C is the angle between 

12 

e and e ’ , C/2 is the angle between P ’-0 and P -0 as the NACSILs 
1 1 11 
are the angle bisectors of the edges. The rest of the angles are 

as shown in the figure. The relations between the angles which are 

shown, can be easily proved on the basis of elementary geometry. 

As P -0 is parall el to e ' , ^ ^ ^ ^ 

1 : ■ 1 

A - C / 2 - D + E = 180; 

in triangle Q-Q -A, 2 * ♦ is the exterior angle, hence, 

• 1 2 
' E = (j> - B ; ^ 
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Substituting in the first relation, 


A - C / 2 D - B / 2 + (j> = 180. 

For smallest A, C, D, B all should be zero. i.e. 
A = 180 - 4) f 


min " t 2.6.14.1). 

As A becomes smaller, <J» increases. At A = 150, ^ = 60. It will be 

proposed in section 2.7 that at around <|> = 60, the proximity and 

size qualitative predicates will be so adjusted to treat the fork 

pair f f as a large proximal pair, to be merged into one fork, 

X ^ 

or an uneven size pair, to be treated as rounded corner. Due to 
this treatment, the angle 4> must be quite small to generate a 
LL LL pair, which will lead to a sparse string. As observed 

above, only sparse strings lead to elongations. Hence, the angle 4> 
must be quite smaller than 60, i.e. A must be larger than 150. 
Thus A must be quite close to 180. 

At the terminal elongations also the same situation prevails as 
described above. If the angle A becomes too small, the fork pair 
is either merged or it is treated as rounded corner. Hence, even 
though the A may not be very close to 180, it can not be much 
smaller either. 


Elongation classification is developed in section 4. 

2 . 7 Exact Definitiona of the Qualitative Predicates ; 

The propositions given above obviously require a definition of the 
qualitative predicates which will keep the approximations to a 
minimum for the deacrlption process to give good results. The 
definitions should be such that the approximations are kept to a 
minimum even in the worst case. The critical predicates are 
exactly defined in this section. The exact definitions of the 
remaining predicates are not so critical and are largely self- 
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evident. they should be refined by experimentation on a larfie 
number of different shapes. 

The worst case for the propositions is that of an almost regular 
hexagon with three elongations on three alternate edges. This 
comes closest to a subtree of large and distant forks. It is 
illustrated in the Fig. 2. 6. 7. 2. 

A general elongation is illustrated in Fig. 2.7.1. 0-V -V ’-V ’-V 

12 3 4 

is the elongation of a polygon, corresponding to the pair of forks 

f and f at 0 and 0’ respectively. V -V is an edge the polygon 
12 2 3 

without an elongation. In the figure, 

1(0Q) = h = height of triangle P-V ’-V 

2 3 

1(00’ ) = d = distance between forks f and f ; 

1 2 

ICOV ) = R = radius of the skeletal circle at fork f ; 

1 1 
l(O’Q) = r = radius of the skeletal circle at fork f ; 

2 

From triangles P-Q-V ’ and Q-V ’-O’, 

2 2 

r = h tanC<J>) tanC45 + / 2 ) ; 

1(0P) = R / sin(<>); I(VIP) = R / tan(4>); 

Area of 0-V -P-V = A = R ** 2 / tanC<|>); 

1 4 

Area of triangle P-V ’-V ’ = A’ = h ** 2 * tan(<j>). 

2 3 

Let A = k * A’, k > 1, where k is the approximation area 
ratio. 

Hence, h = (1 / sqrt(k)) * R / tan(<l>); 

and r = (1 / sqrt(k)) * R * tan(45 + (^ / 2). 

Largeness ratio L = (1 / sqrt(k)) * tan(45 + <|> / 2); 

* r 

( 2.7.1) 

and proximity ratio P = d / sin((|>) - 1 / sqrtCk) tan(<|>)). 
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Fifi. 2 . 7 . 2 Larfienesa and Distance of Forks, Fuzzy Distributions 
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For the polygon in Fig. 2, 6. 7. 2, which is closest to a subtree of 
diatant large forks, il> - 30. We propose that the approximation 
area ratio k - 4 is a sufficiently close approximation in this 
case. For the elongation V^-V » -v '-v the position of edge V 

V about midway between an extreme elongation and a minor 

3 * 

elongation. It gives a critical value for L for the transition 

between •largeness’ and ' smallness’ of forks. 

The critical largeness ratio L = sqrt(3) / 2; 

rc 

and the corresponding proximity ratio P = 2 * (1 - sqrt(3) / 

r 

2 ). 

At L f the smallness and largeness should be at least 0.5 to 
rc 

ensure a good fuzzy spawning. The proximity size quality 
distributions is such that the LL LL situation is ruled out 

'I 

at 4> = 30. The desired relation between the size quality 

distributions for L and P against length along NACSIL 00’, for <|» 

r r 

= 30, are shown in Fig. 2.7.2Ca). It can be seen that this 

relation between L and P size quality distributions rule out 

r r 

LL LL at = 30. Thus, a subtree of large distant forks cati not 

exist in the skeleton. 

The exact size quality distributions of L and P are given in 

r r 

Fig. 2.7.2 (b) and (c). These distributions should be tested for a 
large number of situations for fine tuning and refinements. 

3 ■ Uni fork Classification ^ 

A feature analysis approach is adopted for Unifork classification. 

The features are the interior edge angles of the polygon. The 
angles, in conjunction with the radius of the unifork, completely 
represent the unlfork polygon. Hence a unifork description must 
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give the shape classification based the angles and the size in the 

form of the radius. The size may also be given in some other 

parameters such as the lengths of beaks and combinations thereof. 
As stated earlier, the classification is based on combinations of 
the ranges into which each feature is divided. The feature range 
division IS approximate and not exclusively exact. 

It may be noted that the classification based on the feature 

analysis is a fuzzy spawning process. 

lii Effect of Number of Sides in Classification 

The unifork polygons with large number of sides are very similar 
in shape to circles. Uhen the number of sides is small, the 

deviation from circular shape is too large. We propose to analyse 

triangles and quadrilaterals apart from the polygons with 
circular shape. Pentagons and hexagons are border line cases. A 
detailed exploratory spanning of shape space and concomitant 
analysis required for pentagons and hexagons is complicated. It is 
not attempted in the present work. They are treated as circular 
polygons instead. 

The number of sides is decided on the basis of the number of 

’not_low' intensity ACSILs in the unifork. The low intensity 

ACSILs produce very minor and gentle projections in the polygon. 

In Fig. 3.1.1, the low intensity ACSIL 0-V . produces a gentle 

2 

projections in what otherwise would be a triangle. It is isolated 
from other ACSILs and hence has not been approximated by applying 
rules developed in section 2.6.11. 

The decision of the number of sides is a fuzzy spawning process. 
In ambiguous cases, where the ACSIL intensity is not low and not 
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’not_low’, it generates two deciaion trees, one with the ACSIL 

being treated as a low intensity, and the other with the ACSIL 
being treated as a ’not__low’ intensity ACSIL. 

3.2 The Shape Features ^ 

Shape features are the distinctive angles which dominate the shape 
of polygon. Angles are directly related to the quench velocity 
(ref. section 2.2, chapter 2). The angle features can be detected 
by the quench velocities of the ACSILs of the unifork. 

The small angles in the unifork give rise to sharp projections, 
the 'beaks’ (ref. Fig. 3.2.1). The large obtuse angles approximate 
the circular curve. They do not affect the shape of the polygon. 
The middle level angles produce projections not as sharp as the 
beaks. They are called the 'bulges’. They are detected by the 
qualitative predicates ’Small_angles ’ and ’Middle_angles ’ . The 
ranges of the angles are defined as fuzzy ranges : 

beak : 0 < (J> < 60 ; bulge : 60 < < 100 ; 

where is the interior edge angle. 

The feature identification is a fuzzy spawning process. 

3 . 3 Properties of Convex Polygon 

% 

1) The sum of angles of a convex polygon with 'n' sides = 180 (n - 

2 ) . 

2) The sura of any three angles of a convex polygon is greater than 
or equal to 180. If the sura is equal to 180, the polygon is a 
triangle. 

The proofs of the above propositions are trivial and hence 
omitted. These propositions have been used extensively, explicitly 
as well as implicitly, in the classification. 
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3.4 Clasaif ication of Triain^ lo : 

Triangle classification is baaed only on beaks. Upto three beaks 

may be present in a triangle. Beaks are further subdivided into 

three ranges : 

Prominent beaks (PB) : 0 < <j> < 30 ; 

Middle level beaks (MB) : 30 < 4 . < 50 ; 

Incipient beaks (IB) : 50 < (j> < 60 ; 
where 4) is the interior edge angle. 

The combinatorial classes are : 


1) 

One 

Prominent 

beak 

(ref. Fig. 

3.4.1(a)) . 

2) 

One 

Middle beak (re 

£ . Fig. 3.4 

• Kb)). 

3) 

Two 

Prominent 

beaks 

(ref. Fig. 

3.4.1(c)) 

4) 

PB • 

- MB 

(ref . 

Fig. 

3.4.1(d)) . 


5) 

PB 

- IB 

(ref . 

Fig. 

3.4.1(e)). 


6) 

MB • 

- MB 

(ref . 

Fig. 

3.4.1(f)). 


7) 

MB 

- IB 

(ref . 

Fig. 

3.4.1(g)). 


8) 

IB 

- IB 

(ref. 

Fig . 

3.4.1(h)). 


9) 

IB 

- IB 

- IB 

(ref . 

Fig. 3.4.1 

CD). 


3 . 5 Claaaif ication of Quadrilateral t_ 

The approach employed in classification ia different from that for 
triangles. As a quadrilateral has four angles, the combinations 
would mushroom into a large unmanageable number. The angles used 
for classification are not wholly the edge angles of the polygon. 
The angles between the ACSILs are also used in the exploratory 
spanning of the quadrilateral shape space. 
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3.5.1 Shape Space for Quadrilaterals 

All the angles and their nomenclature are shown in the Fig. 



Fig. 3. 5. 1.1 Nomenclature for Quadrilateral Shape Feasible Space 


The constraints on the angles are, 

0 < © < 180 ; 0 < <t> < 180; i = 1, 2, 3, 4. 

i i 

(___ 

The governing equations for the angles are, 

Q 4 0 4 0 4 0 = 360 ; 

12 3 4 

(• 

From triangle V -0-V , 

1 2 

<j) /2 + <j> /2 4 0 ~ 180 ; 

12 1 

(— — 

From triangle V -0-V , 

2 3 

4>/2 + ij>/2 + e = 180 ; 

2 3 2 

c — 


3.5.1) 


3.5.2) 


3.5.3) 


3.5.4) 


From triangle V -O-V , 

3 4 

<j> /2 + <1> /2 + 0 = 180 ; 

3 4 3 



3.5.5) 
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From triangle V -0-V , 

4 1 

<I>/2 + 4‘/2 + 8 =180; 

4 14 


3.5,6) 

Equations 3.5.2 to 3.5.6 represent a system of linear equat- • 

'^^ions, 

with five equations and eight variables. Three variables, 

^ may 

assume random values, subject to the constraints expressa.^ 

by 

relations 3.5.1, to generate a solution i.e. it has three 

'degrees 

of freedom. Quadrilateral shape space exploration is don 

*ie by 

systematically spanning the solution space over three i. 

basic 

variables . 


© , 9 and (|>1 are selected as the basis variables. The rema,* 

1 4 "Gibing 

five variables can be expressed in terms of these variables . 


e = 180 - 0 ; 

2 4 
0 = 180 - 0 ; 

3 1 

<{> = 360 - 28 - <j> ; 

2 11 

4>=4f+e+e - 180 ; 

3 114 

4> = 360 - 28 - <j> ; (--■ 

4 4 1 

Expressing the constraints in 3.5.1 in terms of these 

variables : 


(• 

c- 

C- 

c- 


3.5 6) 

■"■3,5.7) 

^■5.8) 

~'3.5.9) 

3.5.10) 


bas 


IS 


20 + 4> < 

360 

1 

1 1 



20 + 4> > 

180 

f 

1 1 



4> /2 + 0 + 

0 

< 

1 1 

4 


<}► / 2 + 0 + 

0 

> 

1 1 

4 


20 + 4> < 

360 

> 

4 1 



20 + ^ > 

180 

11 

f 


■'^'5.11) 

"■^■5.13) 


-3.5, 


14) 


"■“ ——'3.5.15) 

4 1 — — — - — 3 

The feasible quadrilateral shape space boundaries are drawn 

• b in 

Fig. 3. 5, 1.2. The polyhedron ABCDEF is the feasible shape spac^ 
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Symmetry Cpnaiderationa 
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As the nomenclature is arbitrary i.e. any of the anelea could be 
the first angle and the nomenclature could be clockwise or 
anticlockwise, quadrilateral shapes are repeated in the feasible 
space under different nomenclature schemes. To pars down the 
feasible shape space to its exact size avoiding repetitions, the 
nomenclature scheme needs to be fixed uniquely. 

The nomenclature is fixed uniquely by specifying that 6 should be 

4 

the smallest of 0 and 0 should be the next smallest. This 

i 1 


implies that 


1) the feasible space polyhedron is cut by the plane 0 =8 and 

1 4 

the polyhedron with 0 > 0 forms the feasible space. 

1 4 

2) From this polyhedron, the part with 0 >0 i.e. 0 > 90 is 

13 1 


removed . 

The resulting pyramid ABFD is shown in the Fig. 3. 5. 1.3. 

3.5.2 Exploration of the Shape Space 

The pyramidal space AEFD has special quadrilateral categories at 
the vertices, edges and faces. They are illustrated in Fig. 3.5.2. 
In the figure, typical representative cases are drawn. In case of 
some of the categories, when an interior edge angle of the 
quadrilateral is zero, the angles shown are close to zero, about 
10 degrees, but not zero. 

1) Vertex A : 

In the vicinity of the vertex A, the quadrilateral has two very 

sharp beaks at angles and (ref. Fig. 3.5.2(a)). At the 

1 3 

vertex A, 0 = 0 = 0 =0 = 90 ; (t> ‘ == <1> = 0 ; «j> = 4> = 180. 

1 2 3 4 1' 3 2 4 
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(^C) 


F 



\J ■S-'X £■ 



Ve.T \^'f^. 



[pj Nftfe. eti^e. /^-D 


Fig. 3.5.2 Classes of Quadrilaterals 


(continued on next page) 
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eg) 










Fig. 3.5.2 Classes of Quadrilaterals (continued from previous 



2 ) Vertex E : 

In the vicinity of vertex E, the quadrilateral has a sharp beak at 

angle 0 (ref. Fig. 3.5.2Cb)). At the. vertex E, 0 =8 = 45; 9 = 

3 14 2 

0 = 135; 4> = 180; 4> = (j> = 90; (j> =0. 

3 i 2 4 3 

3 ) Vertex F : 

The quadrilateral has two sharp beaks at 4> and (ref. Fig. 

2 3 

3.5.2CC)). At F, 0 =0 = 90; 0 = 0; 9 = 180; ({> = <[> = 180; (}> 

1 3 4 2 1 4 2 

= <j> = 0 . 

3 

4) Vertex D : 

The quadrilateral has two sharp beaks at and c{> (ref. Fig. 

2 4 

3.5.2Cd)). AtD, 0 = 0 = 0 = 9 =90;<|> = 4> = 180; 4 . = <t> 

1 2 3 4 1 3 2 4 


0 . 

5 ) Edge AD : 

The quadrilateral is a rhombus (ref. Fig. 3.5.2(e). Along the edge 

,0 = 0 = 8 = 0 = 90;(|> = <1>;<|> = 4>.At the midpoint of the 

12 3 4 1324 

edge H,<j) = <i> = 4> = <1> =90. It is a square. 

12 3 4 

6 ) Edge AE : 


The quadrilateral is symmetrical about one of the diagonals, a 
shape resembling a kite, with a sharp beak at ^3 (ref. Fig. 


3.5.2Cf)). Along the edge, 0 = 0 ; 0 

14 2 

7) Edge Af : 


0;4) = = 0 . 

3 2 4 3 


The quadrilateral is a trapezoid with a sharp beak at 

Fig. 3.5.2(g)). Along the edge 8 = 0 =90;<t> = 0. 

13 3 


8 ) Edge EF : 


4 > 

3 


(ref . 


The quadrilateral has a large angle, almost 180 at <}> 1 , and a sharp 

beak at 4> (ref. Fig. 3 . 5 . 2 (h)). Along the edge, <j> = 180; 4> = 0. 

3 ^ 3 

9) Edge ED : 
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The quadrilateral is a kite with one large angle, airaost 180, at. 

(j) (ref. Fig. 3.5.2Ci)). Along the edge,, 4> = 180. 

1 2 . 

10) Edge FD : 

The quadrilateral is a trapezoid with a large angle and a sharp 
beak (ref. Fig. 3.5.2(j)). Along the edge, 0=0= 90; <i> = 180; 


2 

11) Face AEF : 

The characteristic of face AEF is that the angle 
quadrilateral is shown in Fig. 3.5.2Ck). 

12) Face EFD : 


4 ^ 

3 


0. A typical 


Face EFD is a degenerate face representing triangles as the angle 

(J) = 180. A typical quadrilateral lying very close to the face is 

1 

shown in Fig. 3. 5. 2(1). The quadrilaterals lying very close to the 
face and near the center will be approximated as triangles. 

13) Face AED : 

The quadrilaterals on this face are 'kites’ (ref. Fig. 3.5.2Cm). 

14) Face AFD ; 

For this face, 0 =0 =90 which implies that the quadrilaterals 

1 3 

on this face are trapezoids. The proof of this statement is 
trivial. A typical quadrilateral lying close to the center of the 
face is shown in Fig. 3.5.2(n). 

3.5.3 Regions in Shape Space and Classification 

The regions are those close to the vertices and point H 
representing square, and the edges and faces. The regions are 
defined as qualitative predicate values of the specifications of 
the points, edges and faces. The predicates are defined for faces 
first. The edges are an 'AND' conjunction of the faces which 
intersect to produce the edges. Similarly, the vertices are 
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defined as conjunctions of the edges. Point H is taken as the 
laidpoint of edge AD. The region between point H and vertex A is 
the same as the region between H and D for membership grade for 
this class; otherwise they differ in their relationship with other 
classes . 


The central regions of edges are defined as conjunction of the 

edge predicate and negations of vertices predicates. For example, 

the central region of edge AE is defined as (<{> = 0) AND (<}> = 

3 1 

180) AND (0 /= 90) AMD (0 /= 0 ). All the comparisons and the 

1 4 1 

conjunction are fuzzy. The fuzzy ranges of the angle values are 

shown in Fig. 3.5.3. The fuzzy membership functions shown are for 

the particular quality such as 0 =90. The angles observed in the 

1 

polygon say, 0 , has its membership function as shown in part(d) 

1 

for 0 and 0 . The application of the predicate is the conjunction 
1 4 

of fuzzy sets. 


The classes and their predic 



o eb “SO <5|5 


(a) 0 = 90 



o 165 no 

(c) d) = 180 
1 

Fig. 3.5.3 Fuzzy membership 


es are listed in table 3.5.3. 



o fo 15 

(b) C|) = 0 



(d) e =0 

1 4 

Functions for Quadrilaterals 
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Table 3.5.3 


Quadrilateral classes and their 

predicates : 



Angle Quality : 

01 = 90 

01 = ©4 

01 

03 

1) Interior 

NOT 

NOT 

/= 0 

/= 180 

2 ) Vertex A 

YES 

YES 

= 0 

= 0 

3) Vertex F 

YES 

NOT 

= 180 

= 0 

4) Vertex E 

NOT 

YES 

= 180 

= 0 

5) Vertex D 

YES 

YES 

= 180 

= 180 

6) Edge FD 

YES 

YES 

= 180 

/= 0 

7) Edge AF 

YES 

NOT 

/= 180 

/= 0 

8) Edge ED 

NOT 

YES 

= 180 

/= 0 

9) Edge EF 

NOT 

NOT 

= 180 

= 0 

10) Edge AE 

NOT 

YES 

/= 180 

/= 0 

11) Edge AD : 





a) Point H : 
02 = 04 

YES 

YES 

/= 180 

/= 0 

b) Zone HA : 
02 /= ©4 

YES 

YES 

/= 180 

/= 0 

c) Zone HD : 
02 /= ©4 

YES 

YES 

/= 180 

/= 0 

12 ) Face AFD 

YES, 

NOT 

/= 180 

/= 0 

13) Face AFE 

NOT 

NOT 

/= 180 

= 0 

14) Face EFD 

NOT 

NOT 

= 180 

/= 0 

15) Face AED 

NOT 

YES 

/= 180 

/= 0 
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.3.1^ — — M.ith Bj^her Mtiab Ml £1. Sidea j_ 

The clasji i f i cat ion of convex polygons with raore than four sides 

(note that the number of sides is the one left after various 
approximations) is based on the beak and bulge features and their 
geometrical relations with each other. The definitions of beaks 
and bulges is the same as that given in section 3.2. 

3.6.1 Classification Basis j_ 

The number of beaks and bulges in a convex polygon are constrained 

by the angles subtended by them at the center of the unifork. In 

Fig. 3.6.1, the lines OP , i = l,...,n , where n is the number of 

i 

sides, are perpendicular to the edges E . Let 8' be the angle 

i i 

between OP and OP . Then, 
i i + 1 

0’ +...+9' +...+0’ = 360; 

1 in 

( 3.6.1) 

0- + ^ = 180; ( 3.6.2) 

i i 

where <|> is the angle between Ei and Ei + 1. Thus 0’ and ij> are 
i i i 

directly related. Beaks and the bulges can be detected using 

angles 0’ . The beak and bulge definitions in terms of 8' are, 

i . i 

Beak : 180 >0' > 120 ; C 3.6.3) 

1 

Bulge : 120 >0’ >90 . ( 3.6.4) 

i 

The number of beaks and bulges is constrained by the relation 
3.6.1. This limits the number of possible combinations of beaks 
and bulges and allows combinatorial approach to remain within 
manageable limits. 
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3.6.2 Combinatorial Clasaea 
The possible combinations are, 

1) One beak 

2) Two beaks 

3) One bulge 

4) Two bulges 

5) three bulges 

6) One beak and one bulge 

7) One beak and two bulges 

The cases of three beaks, four bulges, two beaks and one bulge are 
not possible, as they entail very low intensity ACSILs separating 
the features; these cases would be treated as triangles or 
quadrilaterals, as the case may be, with very gentle minor 
projections . 

A brief discussion and illustrations of the classes are given 
below. In the illustrations, the fork point of the unifork, its 
skeletal circle and the perpendiculars from the fork point onto 
the sides are drawn. The skeleton is not drawn unlike the previous 
f igures . 

1) One beak : 

It is illustrated in Fig. 3. 6. 2.1. As there are no other features, 
geometrical relations are absent. 

2) Two beaks : 

There are two classes of polygons with two beaks, depending on the 
geometrical positions of the beaks. In Fig. 3 . 6 . 2 . 2 (a) , the beaks 
are opposite each other , the minimum angle C ’ of section 3.2.1) 
between them being 60. In Fig. 3. 6. 2. 2 (b), the angle between the 
beaks is close to 0. 
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ere are no other features^ 


3 ) One bulge 

It is illustrated in Fig. 3. 6. 2. 3. As th 
geometrical relations are absent. 

4 ) Two bulges 

There are three classes of polygons with two bulges, based on the 
geometrical positions of the bulges. In Fig. 3.6.2.4Ca), the 
bulges are opposite each other, the minimum angle CO ’ of section 
3.2.1) between them being 90. In Fig. 3.6. 2.4 Cb), the angl e 
between the bulges is close to 45. In part (c), the angle is 0. 

5) three bulges 

Uith three bulges, atleast 270 degrees of the total angle are in 
the bulges. The remaining 90 degrees may be distributed as 
30,30,30 in between the bulges, or 0,45,45 or 0,0,90. 
Correspondingly, there are three categories of polygons with three 
bulges. They are illustrated in Fig. 3. 6. 2. 5. 

6) One beak and one bulge 

There are three categories of a polygon with one beak and one 
bulge. The beak and bulge may be close to one another with angle 
between them almost 0 (ref. Fig. 3.6.2.6(a)). The angle may be 
around 35 (ref. Fig. 3.6.2.6(b)), or the angle may be around 75, 
with the beak and bulge exactly opposite each other (ref. Fig. 
3.6.2.6(c)). 

7) One beak and two bulges : 

There are totally five categories of the class. The total hi’ 
contribution by the beak and bulges is around 300. By relation 
3.6.1, the remaining angles must sum to around 60. These 60 
degrees are to be distributed in three angles between the beaks 
and bulges. Hence the categories are. 
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(C) 


Fifi. 3. 6. 2. 6 Unifork with On© Beak and One Bulge 
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a) angle distribution : 0 - 0 - 60; two subcategories due to 

the asymmetry. (ref. Fig. 3. 6. 2, 7, parts (a) anb (d)). 

b) angle distribution : 0 - 30 - 30; two subcategories due to 

asyro^^® Fig. 3. 6. 2. 7, parts Cel and C d ) ) . 

c) angle distribution : 20 - 20 - 20; (ref. Fig. 3.6.2.7Ce)). 

4 . Elongation Classification ^ 

Elongation is a string of long or middle length NACSILs, split 
into two by the unifork point. The two substrings thus formed, 
henceforth referred to as the elongation branches, are almost 
collinear. Successive NACSILs in the branches, have a very large 
angle, almost 180 degrees, between them. The angles in the 
branches, henceforth referred to as bends, result in deviation of 
the elongation from rectilinear form. The branches, may be almost 
rectilinear, with successive bends in opposite directions, or they 
may curve in one direction. These features in addition to the 

lengths of the branches, are used to classify elongations. 

The classification scheme is hierarchical. The number of 

elongation branches, and lengths of the elongation branches give 
the first level of classification. The bend features are then 
combined with the first level classes to give final 
classification. 

Before proceeding with classification, a number of new concepts 
representing some elongation properties and entities need to be 
defined. 

4.1 Definitions of Elongation Froperties 

1) Hean elongation line : Mean elongation line is defined for each 
of the elongation branches. It is defined as the line passing 

through the origin and with inclination which is the weighted 
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average of the NACSILs in th#. 4-u • 

in Tne Dranch, the weights being the 

intensities of the NACSILs. Thp i <=,««+.%. 4-v t j .- 

me length of the mean elongation line 

is the rectilinear distance between the unifork point and the tip 
of the elongation branch. In Fig. 4.1.1. the elongation branch 0- 

"the origin which is the unifork point. The 
line 0-P is the mean elongation line. It is used as a rough 
approximation of the elongation branch. It gives a reasonably good 
approximation as the elongation branch can not have sharp bends 
and does not contain short lines. 

An overall mean elongation line is defined for the whole of 
elongation along exactly the same lines as for the individual mean 
elongation lines for the elongation branches. 

2) Elongation deviation : It is the standard deviation of the 
inclinations of the NACSILs from the inclination of the mean 
elongation line, weighed by the intensities. It is an indicator of 
the closeness of the approximation in the mean elongation line. 

4 . 2 The Number and Lengths of Branches j_ 

The number of branches may be one or two. The length of branches 
is the length of the mean deviation lines of the branches. It is 
divided into fuzzy long and short categories, the normaliser and 
comparator for the qualitative division being the unifork radius. 
The combinations of these give the classes. 

4 . 3 Benda in Elongation Branches 

The bends in the branches are analysed for curve trends. They are 
classified as with a curve trend and without curve trend. 

A curve is indicated by raonotonic inclinations of the NACSILs. 
Local deviations from the curve trends are approximated by cutting 
across the deviations. In Fig. 4.3.1 the NACSILs forming an 
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elongation branch are drawn as the lines OPl, P1-P2, and so on. 

The local deviation from the curve trend, P -P -P , ,is 

2 3 4 

approximated by substituting the segment P -P -p by the line P - 

2 3 4 2 

p . All local deviations are successively approximated. As the 
4 

elongation branches are not complicated, this simplistic method 
can give good results for most of the cases. For the cases where 
it is complicated, for example in an extremely elongated polygon 
having a large number of NACSILs in its elongation branches which 
are quite randomly disposed, the elongation branch will be 
considered to be a line with a large deviations. 

The value of the curvature of the resulting elongation branch 

would indicate the curve trend. The local curvature of a ’curve' 

formed by a set of line segments joined end to end, can be defined 

by considering the line segments to be approximating a curve. 

Consider two line segments, L and L , have lengths S , S ; if the 

1_2_ 1__2_ 

unit vectors along the lines are t and t . (Note that "t , t are 

12 12 

the unit tangent vectors.) If L and L are located end to end, 

12 

then the local curvature could be defined as the length of the 

normal vector N, given by 

N = (t* - t ) / CCS +S )/2) ; 

12 12 

( _____ 4.3.1) 

where the denominator is the length of the region of the line 
segments associated with the 'local curve’ which gives the local 
curvature. 

The average of the local curvatures is assumed to give the curve 
trend. This assumption is justified as the inclination change, and 
hence the range of local curvatures is quite small in an 
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elongation branch. The standard deviation of curvature, weighed by 
the intensities of l.he NACSILs, gives f.he mean deviation from the 
curve trend and thus it is a measure of the confidence value of 
the curvature. 

The average curvature and its standard deviation together define 
the two bend classes, curve trond and rectilinearity . The curve 
is defined as ’large’ curvature and rectilinearity is 
defined as ’low’ curvature. 

4 . 4 Elongation Classes 

As stated ear’l l cn-. the elongation classes arc combinations of the 
features given above. In addition to the combinations given above, 
the dii'-oc'tion of curvature of the two branches, gives cla.oses of 
’curved towards each other’ and ’curved away form each other'. 
Table 4.4 lists tlx! elongation classes. 

Table 4 . 4 Elongation classes 


Symbols : 


L long, S - short, R rectilinear, C curved, 
T curved towards one another 

A curv<;d away from one another. 


One branch : LR, LC, SR, SC. 
Two branches : 


LR-LR 

LC-LR 

SR-LR 

SC-LR 


LK-LC 

LC-LC-T , LC-LC-A 
SR-LC 

SC-LC-T , SC-LC-A 


LR-SR LR-SC 

LC-SR LC-SC-T , LC-SC-A 

SR-SR SR-SC 

SC-SR -SC’Se-T , SC-SC-A 
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c lasiii i f 1 cation , 


_5 . Convex Polygon Deacription ; 

Convex poly£on description is based on the shape 
size parameters of the polygon, and identification of the 
•favoured directions' in the polygon. Shape classification has 
beoTi (ilaborated upon in the discussion so far. The size parameters 
have boon briefly mentioned in the section 1. The concept of 
favoured directions in the polygon has nol. yet been introduced. As 
the size paramo tiers are dependent on favoured directions, the 
concept of favoured directions wi, ] :i bo expounded first. 

5 . 1 Favoured Directions i_ 

Favoured directions may be said to be the directions along which a 

nesting allocation is likely. Favoured directions may be the 

approximate sjymmeti'y and stretch operation directions of the 
polygon. They are utilised in the entire profile description for 

the rfilations of the constituent convex polygons which are the 

result of the polygon decomposition. 

Favoured directions are the atietch axes, e.g. the mean elongation 

linos, the beak and bulge AC,SIL3, and the axes of symmetry. The 

stretch directions are directly available from the earlier 

analysts. Symmetry axes lie in between the adjacent stretcli axes. 

Favoured axes arc geometrical lines located in the plane of the 

polygon . Favoured directions are illustrated in Fig. 5.1.1. In 

part (a) the CSIL 0 -0 is a favoured axis. In part (b). line OP 

1 2 

is a symmetry axis, and hence a favoured direction. 

For the stretch axes, the measure of ’favouredness’ may he the 
intensity of the CSl i. along the stretch axis. For the mean 
elongation 1 ines, it is the sum of the intensities of the flACSILs 
forming the elongation or elongation branches. 
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A favoured directj,.,, ,„ay be forced by more than one features of 
etretch. In this oaoe, the favouredness is the sum of l:he 

contributions by the fcal«rfi<a tc « it ^ , a. ^ a. 

vtiiajres. If squally favoured two feature 

axeb h.^ve almost equal inclinat i oa.s then the angle biooctor of the 

two ajces i.. .il-aj a favoured axis, called an int ermed.5 ate feature 

axis. If the t\.;o feature axes are parallel, then the line 

equidistant to them is th© intermediate feature axis. The 

f avouredn ess of the intermediate feature axis is the average of 

the favourednes.s of the two feature axes. In Kig. 5.1.2, CSILs 0- 

V and 0-V are almost collinear. Line OP is the intermediate 
1 4 

feature axi.s. 


The symmetry 

axes generally lie either 

along a nf. retch feature 

or 

between two 

stretch features. Appro x 

imixte symmetry axis 

is 

the 

angle b.is(;(^ 

:l.(>r of two equally favoured feature axes. 

The 

f avou r ed n es s 

of the symmetry axis 

is the average 

of 

the 

favour edner.n 

of the two feature axes. 

In Pig. 5.1.3, the 

f eature 


axes O-V and 0-V generate the approximate symmetry axis OP. 

1 2 

Intermediate axes and symmetry axes are determined only for highly 
favourtMl axes. 

5. 2 Size Parameters of Shapes ^ 

In the proposed shape elaa.si f ication scheme, the unifork radius, 
the elongations and the information on whtcli of the sides of the 
uni fork form the MACSILs of the elongations, are suff icient to 
determine the approximated skeleton and hence the approximated 
polygon. Thus, from completeness of the doHcription issue, just 
the unifork radius and the elongation branefuis, in addition to the 
classes, are sufficient. 

XJe propose to take the unt fork radius and the length of the 
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^long its favoured directions as the size parameters. As 
favour*--* <i tJiroctions are associated with classos, with each class 
for * * ‘-’l-i the polygon has an acc<!i>l.able class membership grade, 
atleast- one size parameter will be associated. It plays an 
jriipoT-l. .i.ri-t r-ole in matching of profiles. 
hi ^-11 Classification j. 

The (sv-er-a.!! classification is the combination of the unifork 
£ic'a.tion and the elongation classification. It is a 

hierar oil i. ca.1 classification scheme. The classification decision 
tree is £ oj nuid by attaching Idse whole of elongation decision i.i'oe 
to eael* of the leaf nodes of the unifork decision tree. It i« a 
fuzzy wining decision tree. 

The final classes which occur at the leaf nodes of the decision 
tree, mcvy be treated as a de.scriptive string rathf'.T l.han a 
homogeneous entity. For example, the class ’Two beaked circular 
polygon with one rectilinear long elongation branch’ is a 
description of the polygon. The class is a fuzzy conjunction of 
the qua.1 ital. ive predicates forming the description cl.rings. The 
membership grade of the polygon for the whole of the class is the 
application of fuzzy conjunction of the qualitative predica i. , 
The menilic'.rship grades for each of the predicates forming the 

string stire individually preserved. The dominance of any of the 


features repi'csented by the pr 


cdicate is given by the product of 


J IK/I in-t-ensitv of the feature. The relative 
I'.he grstGi© and l-he intensi y 

>,.K?r^a+ton of the predicates with respect to the 
doiainances of a combination ot xnc pi eu, 

at-e define.! Hkewiee. Not ell the possible cosibinatiom. are 
all-o«d. only those .ein'eaenting a continuous path in the decision 
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tr6© cine! 
r el evaxit - 


fo.vraiiig a stejj in the classification process arc; 


Uith i.ho fiiial cia.3si f ications , the area approximations that were 
done Kuiing the approximation process to produce the 

1 f i (.&tions are preserved as the ratio of original area to the 
a.pproxinia'fced area. Different approximation areas ai-c possible as 
some of tin; steps in the approximation process are also fuzzy 
spawning steps, producing different approximations. 

5 . 4 The Description 

Convex polygon descripl.ion consists of the overall shape 
classification, the area approximation measures, size porainet ers 
and the favoured ax<>.‘i. 

The overall shape classification generates multiple class 

membership for the polygon. All the classes must be prenorvsd in 
the dor.f-r Iption in toto. None of the classes by it.Mtil f would be a 
complete shape description. 

The convex polygon description is a complex Hyol.em of multiple 
classes, their memberships, !;h6ir internal dominance relations and 
the size parametars and the favoured axes. 


6 . ConvaK Polygon Description Algorithm 

The convex polygon description process is based on the 
propo.sitions given in section 2. T t is a fuzzy spawning process, 
producing multiple descriptions. The simplifying approximation'! 
are done before l.ho two pronged feature, analysis of unifork 
analysis and elongation analysis. It is structured into 
functionally different steps, Ea<^h step completely deletes a 
part i <vu 1 ric type of manifestation of minor stretch operations, by 
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1 Elongation Curve Approximations 1 

i ! 

/ 

! 

/ \ 

! I 

I Fork Pair Approximations | 

1 1 

\ / 

I 

/ . — \ 

i ! 

I Global ‘imall Fork Approx it, ration I 

I 

1 

/ \ 

I I 

I Fork Simplifications ! 

1 I 

\ / , 


1 Sejieration of Elongations and Uni fork 




'r- ■ ;"'j 

I Unifork Classification 1 

( ) 


E] ongation Class! f icat ion 


I Favoiirad Axes Analysis 


— / 


I Convex Polygon Desci - 1 p i-ion ! 


--/ 


Fig . 


A.l Rlock Diagram of convex polygon Mgorithm. 
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appt oxiinat ing it. The steps from 1 to 4 are appr-oximation steps. 
51 ,rps 3 to 7 are clavssif icat von srteps. Steps 8 and 9 give the 
descr i p ovi - The algorithm is given in the fovin of a block d iagr .im 
in Fig. 6.1. All the sl.ivps are later elaboral.nd upon and presented 
in the form of procedural algorithms in sections 7 to 14. 

The description process is sequentially broken down into the 
following steps : 


1) Elorigation curve approximations : The curves in the NAr.'>ILs 
are detected and replaced. It is based on the observation in 
section . 6 . 8 . The fork - WACSIL tree is d<;void of dense strings 
al. idvG end of this step. 

2) Fork pair approximations : The 1. 1.- -LS, LL Ls and LL--LL 

fork pairs are replaced in this ntep. It is based on unctions 
2.6.1 through to 2.6.7. It is fuzzy spawning process. The LL- LL 
are the small elongations and LL LS are the rounded corners. 


ITi o. r. ta a 1 1 e 1 o n g a t i o n s a i • cr 

.•approximated 

by merging the 

f 0 1 ’ k i 

t and 

adjusting the 

length of 

ACSILs in the 

1 forks. 

It is 

based on 

observations 

in flections 

2.6.5 and 

2.6.7. 

Al 1 -1 

he 

small 


elongations are deleted from the fork - MACSIL tree in this step. 
The rounded corn<n-u are approximated by extensions of appropriate 


edfie.s of the polygon. 

This step deletes all the small elongations and the locally small 
forks. 

3) Global Small fork apprv>i; » mation : The globally small forks are 
app’roxiiaat ed . It. is based on observations in sections 2.6,1 
through to 2.6.5. Al I the globally small forks are deleted in this 
step . 

4) Fork iiiTiipliflcations : The corners rounded by circular curves 
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in fuJ K-. -iT V detected and replaced. It is based on obsci'v.-il.ion in 
section 2 . 6 . j. 1 . A.g a result of tiiiu step, in conjunction wil.h step 
2 , all the rounded corners are deleted fv<,m the fork - WACSIL 

l.r G e . 

5 ) separation of elongations and unifork : Tlie major fork is 
located. All tne other forks are merged into it and the MAC'lfl.s 
are separ,-ited from the resulting unifork. 

6) Unifork clasn i f ication : The beak and bulge features of the 
uni fork are identified and a classification based on the nuinbcn- of 
sides and the features is accomplished. The classification has 
alri'.ady been proposed in section 4. 

7) Elongation classification : The elongation line diagram is 
classified on the of the. .-mgular distribution and the lengths of 
the elongations. The elongation classification has ali-e.ady been 
proposed in section 5. 

8) Favoured axes analysis : The favoured axes are determined on 
the basis of the propos i I. i ons given in section 5.1. 

9) Description ; The final classification is done as a combination 
of unifork classification and elongation classification. The 
favoured directions are noted and the complete <3cscription is 
«t;ored. 


Control of the description process j. 

The process is modeled as a fuszy spawning process ft of. chapter 
1, section 6). The steps given above, are diffeiout sequential 
parts of the process, which may be viewed as independent 


subproeesses . 
The f uazy quali 


lal. i ve predicates foi- each of the subproces,ses are 
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elaborat <'•<! i vi the respective discussions. The pruninfi of bhe 
process is bancsd on the membership grade returned by the fuzzy 
qual'ita. ! i e predicates. The goal state evaluation should be done 
by u iii.ipping function of the quality grades. The suggested, mapping 
funct ion i -u the weighted average of the fuzzy quality gradon, the 
weights being the intensities of the CSILs or profile boundary 
edges involved in the quality evalualion. 

In the artificial intelligence paradigm, convex polygon 
classification is a heuri.si. i rally guided generative space search. 
The sol u I- ion generation is done by the application of the fuzzy 
qualitative predicates and the rules for the action to be taken 
thereof. Backtracking is controlled by the mapping function which 
is viewed as a heuristic function guiding the backtracking. 


7^ Procedure o± Elongation Cur v^ i- 

Ih. <:onal^ts of a uolquo ideotl f .1 o;. , i ,.n of curve baoed 

on a string and the string doUmitation; followed by its 
replacement. This process is not a fuss, ..-pawning process. It is a 

direct irrevocable algorithm. 

7 . 1 Elonaation Curve Tdentif ication j. 

•vsniitr 1on;i. dense strings of locally 

1) Identify all the numerically long, 

, ra i-ifv the curves accociated with it by the 

large forks. Identify the 

Tnrv<=i"fc of llHo fOirliCS 

low intensity ACSILs eraanatmti f ^ 

pr6seBCbi o f: low id.*- « ii«>*** 

. The curve is delimited by the extreme low 

• and ending on, the .curve- 

intensity ACS I Lis . 

onrlosed by the string and its associated 
2) Calculate the arcan enclos y 

I M Ha that a string may have two curved 
boundary curve segroen''- Not 

+WO areas will be calculated in such 

boundary segments and hence 
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cases. If one of the curve has alreadv Ixuvn replaced, it will not 
be coTinidored in further processing. 

3) v->oi I. f.fie strings in ascending order of areas. For' ni.T'ings with 
two uriT' <‘::o 1 ved curves, the smaller of tluj Lwo areas will be 
considered for sorting. 

4) Replace the curve associated wi l.h the first string. If two 
unresolved curves are associated with l:he string, the smaller area 
curve will be replaced. Hark lh<; boundary as resolved. 

5) If all the curves are replaced, the process is !:oiiii)lete; 
otherwise go to step 3. 

7 . 2 Curve Replacement ^ 

Curve r epl ac oniont is by external envelopes. It in an iterative 
process. 

1) K>:!. find the two extreme edges till they meet (ref. Fig. 7.2). 

2) Det ©rrainc: ihc additional area, the shaded area in the 7.2. 

If the add i 1 . 1 onal area is .small compared to the area enclosed 
within the string and the curve, then the curve replacement is 
complete. The smallnefm of additional area is a qualitative 
predicate for ’good apprsjx 5 wation’ . 

3) If the additional area i .s not small, then split the curve* into 
two. Replace the two curves. For the replacement of these two 
curves, the good approximation criterion should be applied to the 
total area rather than individual application to each of the 
curvet! . 

4) Curve splitting : The curve is spl 5 I- at the point on the 
boundary which is at the maximum distance from the nhord (the 
dotted line in Fig. 7.2). 
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5) The nu^e-rically short curves', are approximated as rounded 
corner.^ as in the step 7. of the section 8. 

7 . 3 Closeness of Approximation 

ThJ.£3 criterion is adjusted as per the requiT-tuiitvn i:s of section 10. 

A circle can at worst be reduced l.o a hexagon. In Fig. 7.3, a 

closenes-s appi ejx i inat ion would be the ratio of the shaded area to 

the area of the sector 0-V -V . The ratio R , is given by, 

12 g 

r**2 / sqrt(3) II * r**2 / 6 

£ 

g 

r **2 / 6 


= 2 * sqrtC3) -II =0.3225; 

Renee the closeness can be defined as ’less_than 0.3225 
8 . Procedure of Fork Pair Approximation j_ 

The LL--LS, LL Ls fcjrks are detected as rounded corners . and 

LL — LL forks as ama ! 1 eloi^gations . The rounded corners are 
by extension of the polygon. The small elongations 
are approximated by merg’ffi the forks. It is a fuzzy sp.iwi> i ng 
process. The fork pairs which have a reasonably good confidence 
value of both the classes spawn two children processes. It is also 
sequential process in the sense that the fork.-, are examined one 

by one. 


1 ) 

If 


Examine 
i !: i. s 


a fork for LL--LS, LL 
LL — LS or LL Ls then 


--LS or LL--LL 
go to step 2. 


categorisation . 

If it is LL — LL 


go to step .3 . 

2) LL-LS and LL Ls rounded corner replace.,™ i . Extend the 

ed«es the NACSIL in the pair till they .»oet . Delete the 

snail fork fron tho fork - NACSIL tree and the skeletal tree. In 


the I'ifi. 8 


1 the small fork £ in deleted as show., hsr deleting the 


edgesl 6^ a^^ e 

12 34 


till they Jiieet at new vertex V 
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3) LL LL sinotll elongation merging : Merge the two forka into one 
fork at the point on the NACSIL dividing the NACSIL into a ratio 
of into.ruii l.iea of the two forks. The CSILa in l.he two forks also 
deviate accordinely (see Fig. 2.6.12.1 and Fig. 2.6.12.2 for ihe 
exact deviation in the Cr>!l,,s). 

li. Procedure for Global Small Fork Approximation i. 

It is a direct sequential process in which the globally small 
forks are detected and approximated as rounded (mrners. 

1) Examine a fork for- global smallness. 

2) Doberraine the larger of the two neighbouring forks. Delete idle 
fork uTuier consideration by extending the edges forming the NACSIL 
betwoon the fork and the largest neighbouring fork. The e.xliension 
is the .same as that in secttion 5. (ref. Fig. 8.1) 

10 ■ Procedure for Fork Simplifications 2 . 

The circular cur-vtui boundary sections in forks may be rounded 
corners or curves. The curves are replaccul first and then the 
rounded corners. 

The curves are detected by application of tlie qualitative 
predicate stated in section 2.6.11. They are approxima I; ed as in 
section 7.2. The criterion for the acc<'i)!able approximation 
cTnnoness is adjusted so bhat a circular polygon with large number 
of sides does not get reduced l.o a quadrilateral, but remains a 
circular polygon. It is a r.iwple sequential process. 

The rounded corners ident.if leal. ion and replacement is also a 
simple sequent ia.l process. - 

1) Rounded corner detection ; Examine a fork for a circular 
rounded corners, by applying Khe qualitative predicate for 


176 



cl €? 1- >>- c I t 1 Oil. cii:i s “t b ,. t cd. in b v . cfelon 2 6 11 

2) Replacement ; II. is the same as that in section 5. 
iiju i£o .£. gJM , re for. aeparation of. Elongation and Unifork ± 

The sep.-iT-al. 1 on of elongation and tmifork is basically a tree 
traversal prob1<;m. The largest fork is located firsl. and it is 
treated .j.s (.he root of the tree for tree traversal. 

1) Location of tlio largest fork : The fork - NACSIL tree may have 
two forks with the same largest radius. In sucki a ca-se a new fork 
is created at a point on the NACSIL between the two forks, 
dividing the NACSIL in the ratio of the respective fork 
intensities. This fork is treated as the root for .separation, 

2) The tree is scanned in a ((Rj)th first traversal. The WACSII.s are 
roni.-j ded to form the elongation tree. The ACSIL.S ar-c; shifted to 
the root node to f<irm the uni fork. 

12. Procedure for Unifork classification ; 


The uni fork classification is a fu^szy spawning process in all its 
.steps. The rniiaber of sides is determined fitiit. Then the 
approprtal:Q feature detection and classification is followed. 
Feature detection is coiiimon for triangles and polygons with mor’c; 
than four sides. For a quadrilateral, the classification is based 
on the angler, between the ACSILs. 

12 . 1 Determination of the Humber of. Sides a Pol..y:a.Qn i. 

The number of sides in the polygon - number of ’not low’ intensity 
ACSILs. As ejcplained in section 10, if an ACSIL i.s on (.he border 
of ’not low’ intenn.i l.y , two branches of the decision tree are 
creatod, one treating the ACSIL as ’nol; low’ intensity, and the 

it as,, .low intensi,t,y. 
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12.2 Feature Identification j_ 

The beak and bulgn detection is done on the basis of the quench 
velocities of tho ACSILs. The quench velocity range division is 
fuzzy. 

Beak : V > 1 / sin(30) ; 

Q 

Prominent beak : V >= 1 / sinClS) ; 

Middle level beak : 1 / sin(15) < V <= 1 / sin(2S) ; 

q 

Incipient beak : 1 / sin(25) < V <= 1 / sin(30) ; 

q 

Bulge : 1 / sin(30) < V <= 1 / sinC45) . 

q 

The proposed fuzzy membership distribution functions, in terms of 
the polygon interior edge angle (p are given in Fig. 12.2. 

Feature identification for quadrilaterals is the comparison of 
various angles with constants and among themselves. The fuzzy 
membership functions for these features have already been proposed 
in section 3 . 5 .3 . 

The feature detection process' is sequential i.e. the features are 
detecte'd one by one. It is also a fuzzy spawning process. 
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to be hi erar cl) f cfi.'! . 


The 


l£jlLf.9,£lS Cj3.aai f icat i nn . 

Uni fork claaaon ha ^ 

way De considered 

classification r>r 00^530 10-^ r 

a fuzzy spawning tra.vnrnal of a path in 

the decision tree of the hierarchy. The decision tree of the 
hierarchy is very simple to construct. The first decision is based 
on the number of ACSILs which is the same as the numho. of sides. 

unifork is a triangle, then the further decision tree is 
based on the classification given in section 3.4. If it is a 

quadrilateral , i;he fut'+h.ar. ^ ... 

turther decision tree is based on the 

classification given in section 3.5. For a polygon with higher 

Kuxibf.r of sides, the tree is based on the classification in 
section 3.6. 


111. . g U g -C g teLg. for Elongation Classification j_ 

Elongation classification is also a hierarchical scheme. The first 
level of the hierarchy is the number and length combination 
decision. Ihen the curve trend analysis is done and I. he resulting 
combinations with number and length classes are the final classes 
of elongation. 

Evalual. ion of the mean elongation line and the mean deviation 
pT-cnrede the classification process. 

13.1 Number and Lengths of Eloncations 2 . 

Elongation bf.anchea may be on or two in number. The length of 
elongation is the length of the mean elongation line. It is 
compared with the unifork radius for relations of ’much larger’ 
and ’ comparable' length. These two are the length clasmfsn, called 
respec: i.l vely ’very long’ and ’long’ .The combinal. i ons of the 
number and length clas.ueii are the first level classes . In addition 
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to thes-, a category of two very long elongation branches wJtich 
are at an - Tglo with one another and the category of two very long 
collinear elongal. icm branches are also formed. The mean cldngation 
line is considered for the purpose of angle evaluation. 

13 . 2 Curve Trend 

Morioi.onic inclinations represent a curve trend. The average 
curvature value and a low inclination deviation for curvature 
give; the curve trend. 

The inclination trend is followed from the unifork origin till the 
tip of the elorigation branch. At each local rever.sal in trend of 
inclination, incl i nations may be increasing or decreasing, the 
deviation-^ are smoothened out, as stated in section 4.3 (ref. Fig. 
4.3.1). This is continued till a uniform increasing or dec.reasing 
trend is obtained through out the elongation In aneh. The average 
curvature and ' the mean inclination deviation are calculated. A 
curve is detected by a high curvature and low mean inclination 
deviation. The sign of curvature indicates whether it curves 
upwards or downwards. 

The presence or absence of a curve, and its direction of cuT vature 
in relation to the other curve if it is present, form the basis of 
the classes . 

14 . Pro cedure for Favoured Direct ion Analysis 
1) Stretch a-xns determination : 

Stretc;h axes are the ACSIL lines of beaks and bulges, and the 
NACSILs of elongations. The mean elongation line in also taken to 
be a str etch axis . The measures of f avouredn(*.ss of all the stretch 
axes are calculated. 
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2) Inlnr-mediate featurA ^ ^ . 

feature ax8s aetermination : 

PaXjT'Q of ixiffnl.y favoii'Pofl « -1 4- 

feature axes with close inclin.i.1. i ons are 

d c fc e c t e d » 1 1 x e x ix ter iti erfiA*f*A 

te f eatuiT. axes and their favouredness are 

calculated as iil.ah.ed in section 5.1. 

3) Symmetry axes determination : 

Pairs of highly favoured axes with equal favouredness measure, and 
which do not have close inclinations, are detected. The Symmetry 

axes and thexr favouredness are calculated as stated in section 
5.1. 


15 . Procedure for Description ^ 

The classifications are based on the features detected and their 
geometrical relations. The propo.sed classification is hierarchical 
in na.ture. The process of classification can be modnlod as a 
decis.iovi tree with fuzzy spawning. The decisions are based on the 
detected features and tho ijooinetf'icai interrelations among i.h«m. 
ThoJKi have already been expounded earlier. 

The decision trees are different for each of uni fork 
c3 aciii f ication of triangles, quadrilaterals and polygons with 
higher number of sides, and the elongation classification. The 
exact codification of the decision trees and the process of their 
comb 1 iiat ions to form the cuwplote decision tree is self evident, 
from l.he proposed classification scheine. 

The complete convex polygon do-scription consists of the multiple 
classes, their membership grades and the area approximation 
measxtres, their internal dominance relations, favxHirod axes and 
the SI i z 0 parameters. The multiple classes and the memberships are 
generated by the classification proeodure given in sections 7 to 
12. The internal dominance relations within a class are generated 
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in 


the 


sarae process. The favoured axes generation procedure is 
given in section 13. The size parameters are det <u-ini nod along with 
the favoured axes determination. 

The iusZu* .■-■{.’•.iv/iiing can be modeled as backtracking upto a marked 
decisiOii foi' the puT'pose of progra.ia design. Hence essentially the 
process model tor program design is the same as that for polygon 
decompos i t i on . 


£i .Convex Polygon Description x 


The examples 

of 

description 

are taken 

from 

the first 

decomposition, 

of 

the profile 

given in 

the 

example of 


(5 octtKsiposition algori I'.hm (ref. chapter 3, section 9. Fig. 9.3). All 


the !.1 u o 8 CCFa of the profile, C ,, C and C , are described here. 

12 3 

1) CCP C is .shown in Fig. 16,1. It is a triangle with angles 0 = 

1 1 
32, 0 = 56 3.nd 0 =92. As ther<^ is only one fork, no 

2 3 

approximations are needed. It can lx; directly classified. 

Beak analysis (ref. section 3.4) : 

0 : Prominent Beak , 0.45 membership, and Middle Ueak 0.55 

1 

membership . 

0 : Incipient Beak i.O membership. 

2 

Hence the given CCP C is a triangle, of class 

1 

a) PB-IB with 0.45 membership, 

b) Mii-IB with 0.55 membership. 

The feature axes are along tho C.OILs emanating from the fork. They 
are , . 

a) C511L F-V i slKS 21 mm, 

’ * 1 

b) CSIL F~V ; sige 13 mm, 

2 

c) CSIL F-V : siae 8 mm. 

3 
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Feature axis in a) is the most important as 
has maximum intensity. 

2) The CCP C has two forks F and F (ref. 
pair Is LL LL. Hence there are no 

separatee} tin {fork and elongations are shown 
polygon is a quadrilateral. 

Unifork Analysis (rof. section 3.5) : 


the corresponding 

Fig . 16.2). The 
approx i mat i ons . 
in Fig. 16 . 2 (b) . 


CSIL 


fork 

The 

The 


0 = 87, 0 = 109, 0 = 93, 0 = 71; 

12 3 4 

(J> = 130, = 56, Ip = 88, <p =86. 

1 2 3 4 

Appl t i on of the rules for classification : 

a = 90, e /= e , 4> /= o, 0 /= iso. 

1 1 4 14 

Hence from table 3.5.3, the convex polygon is of unifork type 12, 
with unifork radius of 12.6 mm. 

Elongation analysis : 

The polygon has only one elongation line (line F - F in Fig. 16.2 

1 2 

(a)). It’s length is comparable to the unifork radius. Hence it is 
’long’ elongation (ref. section 13.1). There me no curves in 
elongation. 

Feature axes : 

The feature axes are, in descending order of importance, F -F , 

12 

F -V , F -V , F -V , F -V . The sizes along the feature axes are, 

1 2 2 3 2 4 1 1 

F 'F : 18 ram, F -V : 13.5 mm, F -V : 25.5ram, F -V : 19 mm, F - 

1 2 1 1 1 2 2 3 2 

V : 19mm. 

4 

3) The CCP C (ref. Fig. 16.3) is a convex polygon with large 
3 

number sides. It has no elamgation curves, as there are no long 
strings of local ly large prox iiiial forks (ref . section 7.1). T l. has 
the following fork pairs. 
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F -F : Largeness ratio L = 0.934, 

12 r 

Proximity ratio P = 0.23, LL- Ll.; 

r 

F -F : Largoiiass ratio L =0.961, 

2 3 r . 

Proximity ratio P = 0.205, LL — LL; 

r 

F -F : Largeness ratio L = 0.738, 

3 4 r 

Proximity ratio 1' =0.57, LL LS; 

r 

F -F : Larg<iTioss ratio L = 0.748, 

4 5 r 

Proximity ratio P =0.4, LL LS. 

r 

Fork pair approximation C ref. section 8) : 

F -F and F -F are approximated by rounding of the corners, i.e. 
3 4 4 5 

by extending the edgen V -V and V -V till they meet at V 

94 12 11 13 

(ref. Fig. 16.3 (b)). F -F and F -F are now negligibly nmall. 

12 2 3 

They are merged in to the CSlI. F -V 

3 13 

There are no more approximations. The next nl.ep is separation. The 
unifork radius is 15 mm. The WAC’.SILs are very minor, and they have 
already been merged in to tht! neighbouring CSILs. Hence 
elongat j onn are absent. 

Unifork Classification : 

There is only one beak at V , <j) = 54, Incipient Beak. 

13 

Hence the given CCP is a polygon with large number of sides and 
with one incipient beak. 


Feature axes : 


One important feature axis along CSIL F -V , size 32 mm. 

3 13 

The low intensity axes are : 

F -V : size 15.5 mm, F -V : size 15.5 mm, F -V : size 15.5 mm, 

1 4 , 13 12 

F -V : six<! 17 mm, F -V : size 15mm. 

2 1 3 12 

The profile and CCPs after approximations are shown in Fig. 16.4. 


184 




Fig. 16.1 Description of a Triangle 
Vs 



F 

- V : 

132 , 

15.5 mm ; 

F - V 

: 130, 

15.5 mm ; 

1 

4 



1 3 



F 

- V : 

H 30 , 

15.5 mm ; 

F - V 

: 125, 

1 7 ffna ; 

1 

2 



2 1 



F 

- V 

: 150 

, 11 mm ; 

F - V 

: 153 

, 8 mm ; 

3 

12 



4 11 



F 

- V 

: 153 

, 8 mm ; 

F - V 

: 107, 

10 iin»; , 

5 

10 



5 9 



F 

- F : 

80, 

.5 mm ; 

F - F 

; 54, 

8 mm; 

5 

4 



4 3 



F 

- F : 

54, 

8 mm ; 

F - F 

: 33, 

4 mm . 

3 

2 



2 1 



F 

- F : 

I. 0 

.934, P 0.23; 

LL— LL 



1 

2 

r 

r 




F 

- F : 

L 0 

.961, P 0.205; 

LL- LL 



2 

3 

r 

r 




F 

- F : 

L 0 

.738, P 0,522; 

LL 

LS 


3 

4 

r 

r 




F 

- F : 

L 0 

.75, P 0.4; LL 

LS 



4 

5 

r 

r 





Fig. 16.3 Description of a Polygon with High Number- of Sides 
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Chapter 5 

PROFILE DES CRIPTION 

1 . Introduction 

The approach of the present work towards nesting is based on a 
shape and size matching of profiles. A defii:i- iption of the profiles 
to be allocated is to be matched against the description of the 
compl imen l.ary area of the allocation sheet fni-med by previous 
allocations, to maximise nesting efficiency. Profile description 
is the crux of the approach. 

Proposals for decomposition of profile in to the ’Constituent 
Convex Pol. ygons ’ , henceforth abbreviated as ('CPs, and description 
of the CCPs have been outlined in the previous chaptern. After the 
conclusion of these analyses, the profile model is in l;he form of 
a union of the CCPs, and their descriptions. The proposed profile 
des«ir i p l:ion is based on the following ; 

1) A further analysis of the concavities of the profile, 
henenf <»!• l;h referred to as bays. 

2) Idenl.i f i cot ion of shape features likely to be ustsful in 
nesting. Only the projections of the profile are considered in 
the present work. Boundary features such as long straight edges 
for edge justification in jnx I aposition of profiles has not been 
considered. 

3) A study of the geometrical and topological relations among the 
projections and the bays of the profile. 

Description of the profiles for matching and concomitant nesting 
is enuraerative in nature rather than a classification as was done 
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for the convex polygons. A classification is not attempted i ti the 
present work as it would be extremely complex. 

ZllS ^ 9 HQquir eroents be satisfied by prof i 1 e description ^ 

A profile may be rotated, displaced along the X and Y directions 
in th<i nesting plane, or it may flipped for allocation. Hence the 
description of a profile must be independent of rotation, 

displacement and mirror imaging transformations in the nesting 
plane. 

2 . Bay Analysis j_ 

The aim of the description of the shape of a profile is that the 

description should be useful in the nesting allocation of the 

profile, A nesting allocation needs an identification of the 

'projections’ and the 'bays' of the profile which may be fitted 

snug against the bays and projections respectively of the profiles 

already allocated. The projections and bays of a profile could be 

matched against the projections and bays of the 'complimentary 

profile’. In the Fig. 2.1, the profile P , is being allocated. The 

f 

shaded area is the region of the sheet on which the earlier 

allocations of profiles have been done. The unshaded area of the 

sheet is the complimentary profile with which the profile P is to 

f 

be matched, and on which P is to be allocated. Tt can be seen 

f 

that not only the projections of P and the complimentary profile 

f 

may be matched, but also the concave parts of the two profiles may 

be matched for nesting efficiency. P has projections P , P , P 

f 1 2 3 

and bays B , B lying in between the projections. The projection 
12 

P is matched with a projection of the complimentary profile. The 
1 

projections P and P do not match exactly with the projections of 
2 3 

the complimentary profile but the bay B matches with a bay of the 

.2 
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complimentary profile. 

A d escf i pt. ion of a profile can not be truly complete without a 
consideration of the bays of the profile. The bays give the 
int fractions of the profile boundary across the exterior (ref. 
section 3, chapter 2), 

2 . 1 Discussion on bay analysis 

Bays for-iii an integral part of the shape of a profile. Bays are 
regions of the exterior of the profile delineated V>y the 
projections of the profile. Alternatively, bays can be thought of 
as the projections of the exterior of the profile or projections 
of the complimentary area of the profile, in the plane of the 
profile. 

Bays could be detected by a decomposition of the complimentary 

area of the profile along the same lines as the decomposition of 

the profile, A sintp>lfr approach is proposed in the present work- 

Bays may be tiofined as the complimentary regions of a profile in 

the convex hull of the profile. In Fig. 2.1.1, the convex hull of 

the profile P is shown in dashed lines. B and B are the 

f 12 

irampLimeat.acy regions of the profile in tha CQrv.ej£ hiilL.. They are 

the bays of the profile. This definitiors sixactly delineates bays 

as two dimoTis ional profiles which are exactly located in the plane 

with respect to the original profile. The problem of identifying 

the bays is that of obtaining the set of cycles in the geometrical 

graph formed by the boundary of the profile and the convex hull of 

the profile. Due to the definition of the convex hull, it is 

guaranteed that the bays, i.e. the cycles in the aforementioned 

graph', are mutually exclusive. The problem of idantification of 



such cycles is a standard problem and standard solutions exist .in 
hhe literature [12]. 

As bays thenisclves are profiles, they can be analysed in the same 
way .i-s the profilos. The integration of the two analyses may be 
used as a description of the profile. The shape analysis of a bay 
differs from the shape analysis of a profile only in the details 
of the approximations done in the analysis process. The 
approximations in the analysis of a profile are of the form of 
external envelopes; whereas the approximations in the analysis of 
bays are, what may correspondingly be called as the ’internal 
envelopes’ of the polygon. The exact meaning of this term will be 
clarified in the next section where the ’internal envelope 
approximations’ are proposed. 

The difference between approximation method arises because of the 
basic nature of nesting. The most basic constraint on allocation 
of pr-ofiles is that there should be no overlap of the profiles. 
External envelope approximations in shape analysis of profile 
guarantee that if the approximated profile is allocated, checking 
the constraint satisfaction, there would be no overlap of the 
exact, profiles. Similarly, internal envelope approximations in bay 
analysis ensure that the approximate bays and their descriptions 
may be used in allocation. 

2 . 2 Bay analysis algorithm vis a vis profile shape analysis 
As stated above, bay analysis is the same as the profil<s shape 
analysis, except for the approximation methods, it suffices to 
just give the ’internal envelope approximations’ of the bay, to 
lay down the proctuJural algorithm for bay analysis, after the bays 
have been identified. The algorithm is as follows : 
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2.2.1 Identif ication of bays j_ 

1) Detorinine the convex hull of the approxima t. (>.<1 profile. These 
approximations are the external envelope approximations outlined 
so far in chapters on profile decomposition and convex polygon 
description, 

2) Bays are the complimentary areas of the approximated profile 
with its convex hull. 

3) Bays are analysed as independent shapes upto and including the 
neighbourhood relationship analysis given in section 3. 

2.2.2 The approximations in shape analysis of bays 

1) Approximations in decomposition of bays into constituent convex 
polygons : 

The ’straightening up’ approximations (ref. section 6.5.1, chapter 
3) are the same as 'cuts’ (ref. section 6.5.2, chapter 3). A cut 
divides the profile in to two or three constituent polygons. The 
criterion to distinguish between a cut and a straightening up 
approximation is that one of the polygons resulting from the cut 
should have a small area in comparison with the area of the other 
polygons. The polygon with small area is neglected altogether. 

2) Curve approximation in convex polygon description : 

Curve approximat i ciius (ref. section 7, chapter 4) for bays are 

illustrated in Fig. 2. 2. 2.1. The curve V -V ’ is approximated by the 

11 

pair of lines V -A and A-V ’. The point A is the point on the 

1 1 

curve V -V ' at the maximum distance from the chord V -V ’. The 

11 ... ^ ^ 
procedure for approximations and the definition of closeness of 

approximation remain the same. 
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Curve V - V ’ approximated by 
1 1 

line scgiiicnts V - A V ' 

1 1 

Point A is at maximum distance 

from V - V ' 

1 1 

2. 2. 2.1 Curve Approximation in Bays 



2. 2.2.2 Fork Pair Approximation in Bays 



3) Fork Pair Rounded Corner Approximation : 

It is done for unevnn size forks, as in convex polygon description 
(ref. section 8, chapter 4). In Fig, 2. 2.2.2, f is the large fork 

and f i.s the small fork of the part of polygon V -V -V -V -V -V . 

^ . . 1 2 3 4 5 6 

The approximation linos e ’, e e ' ' ' from the upper s i de of 

111 

thn polygon and e ’ , e ’ ’ , e ’ ’ ’ from the lower side of polygon, 

6 6 6 

form the possible approx J tna Lions . Of all the approximations, the 

approx 1 i ons with roim’inum area approximation and with i:he 

approximated ama equally divided on either sides, am .nolected. 

For the approximated area being equally divided on either sides, 

in Fig. 2. 2. 2. 2, if the approximation e ’-e ’’’ has minimum area 

1 6 

then, the areas V -V -V and V V -V -V , should be approximately 

1 2 3 6 5 4 3 

equal . 

'4) Rounded corner at a fork ; 

It is the approximation given in section ID of the chapter 4. It 
is exactly the same as that given in the pari, two of the current 
section. 

2 . 3 Niche of bays and bay analysis in profile description ^ 

After the identification of bays and the analysis of the bays up 
to the convex polygon description stage as expounded in the 
previous chapters, the shape knowledge extracted from the profile 
is in the following state. 

The profile model is in the form of a convex hull. The convex hull 
is a union of the approximated original profile and the bays of 
the profile. The bays are decomposed in to their constituent 
convex polygons (CCPs). The descriptions of the CCPs of the bays 
and the profile are present in the profile model. The bays have 
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diverse geometrical and topological relations with the CCPs of the 
profile. The geometrical and topological relations among the bays 
and the profile are important for nesting. The formulation and 
analysis of these relationships are given in the n<ix.t section. 

The bay identification and analysis is performed after the 
projection relation analysis, proposed in the section 3.2 of the 
present chapter. The bay analysis, as .stated earlier, differs very 
little from the profile analysis. It consists of the profile 
analysis , with the differences stated in section 2.2.1, up to and 
inclusive of projection relation analysis. 

2 . 4 Bay Merging _£. 

In the course of nesting, some of tlie identified bays may not be 

useful. It may happen that, a strong match of some other 

projections and bay features may preclude the utilisation of some 

of the bays. Such a -situation is illustrated in Fig. 2.4.1. A very 

good match of the dominant projections P , P and ttie bay B 

1 2 1 
forces an allocation which wastes tho area of the bay B . Such a 

2 

situation may also arise for highly complicated shapes. It may 
also happen that a matching bay is not found at all in the 
complimentary profiles formed during the allocation. To obtain a 
shape ma(;ch match in such situations, it is required to merge some 
of the bays into the profile. After merging a bay in to the 
profile, the new profile, referred to as ’Bay Merged Prof i 1 e' (BMP) 
henceforth, is subjected to shape analysis to obtain its 
description. 

A model of the profile in which a number of BMP versions of the 
profile are formed by merging combinations of bays in to the 
profile, can generate descriptions to satisfactorily handle the 
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cases described above. 



Description o£ a profile, as stated earlier, is an enumeration. 
The geometrical and the topological relations among the CCPs and 
bays, along with the CCP description, const i l.u l.<! the profile 
description. All the CCPs do not play an equal role in 
description. The CCPs, which project out of the profile and form 
bays between them are more important for nesting. These CCPs, 
called projections, are iderit. ified on the basis of their 
topological relations. 
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3 • 1 Topolo£ical relations &mon& the CCPs : 

3 ■ 1 .1 Ba-sic Idea of topological relations j_ 

Topological relatiorii-i are the connectivity relations among the 

CCPs, indicating the neighbourhood of the CCPs. Classically, these 

relations would be represented by an adjacency graph of the CCPs; 

where in the individual CCPs are the nodes of the graph and the 

existence of a shared edge betworn two CCPs producing an edge of 

the graph b<; tween the two nodes. Such a graph, does not 

satisfactorily model the relationships among the CCPs. In Fig. 

3. 1.1.1, CP and CP arc two CCPs. In part (a), they share an 
1 2 

edge; and hence they would be depicted as neighbours in the 

adjacency graph. In part C^) they do not share an edge. They are 

separated by a thin wedge. Here, they would not havo any 

neighbourhood rnlaUionship betwrion them. This does not discern 

between the neighbourhood relation of CP with CP and some othft 

1 2 

distant CCP. Also, an adjaooncy graph does not distinguish between 
the two oiluations depicted in parts (a) and (b) of Fig. 3. 1.1. 2. 
These two cases are distinctly different, and should be 
represented distinctly from one another. 

Different types of neighbourhood relations are defined in the 
present work, whidi overcome these drawbacks. They I'epresent 
actual connectivity relations more accurately than the .'wi jacency 
graph. 
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^ of. the neighbourhood relatione : 


The neighbourhood relations should fticorporate the folloving 
conceptji : 

1) The importance of a shared edge, from tiu; perspective of both 

the CCPs in the relation. The importance of the shared edge may be 

different in the them. In Fig. 3.1.1.2(a), the shared edge is more 

important in the CP than in CP . The neighbourhood of CP is more 

12 2 
important to CP that that of CP to CP . 

1 12 

2) Tunneling effect : Two CCPs may not share an tedge, but they may 

be separated by a thin wedgeCref. Fig. 3.1.1. 1). The CCPs CP and 

1 

CP in part (b) of the figure should be considered as neighbours. 
2 

It can be viewed as a tunneling of the neighbourhood relatiorifdiij! 
across a barrier. This maintains a continuity of the neighbourhood 
relation as the two CCPs move a.nart . The quantity representing the 
neighbourhood effect should reduce as the separ. tion between the 
CCPs increases. 

3) Gcoiiiotrical location of the neighbourhood relation : A CCP may 

share more than one edge with other CCPs. In Fig. .3. 1.2.1, CP is 

1 

a neighbour of CP and CP . The neighbourhood relations of CP 

2 3 1 

with CP and CP are exactly located in the plane of the figure. 
2 3 

They can be said to act along the lines L and L . L and L are 

2 3 2 3 

the perpendicular bisectors of the respective shared edges. 
Definition of neighbourhood relations, henceforth referred to as N 
relations, incorporating the concepts stated above, is based on 
the following discussion. The profile model consists of the 
unapproximated CCPs, which are the result of the decomposition. 
The skeletons of the CCPs form a geometrical tree. The edgor. of 
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the CCPs in combination with the skeletons forii! a geometrical 

graph, henceforth referred to as the CCP - r.knleton graph. A CCP 

skfileton graph is illustrated in Fig. 3. 3. 2. 2. The cycle C -C - 

1 2 

C -C -C -C -C -C -C , i>,-i.osGS through tho CCPs CP , CP , CP . The 
3456 7 8 1 123 

shaded area is the area enclosed with in the cycle. The area C - 

1 

C -C is a subarea of the cycle. It is the area of the cycle 
2 8 

associated with the CCP CP . The subarea C -C -C is associated 

1 3 4 5 

with the CCP CP . The cycle establishes a neighbourhood relation 

2 

between CP and CP , tunneling through CP . In CP , the area C - 
12 3 1 1 

C -C is associated with the cycle, out of the total area of the 
2 8 

polygon. The area C -C -C -C -C -C -C -C is the subarea of the 

12345678 

cycle, linking CP with CP . There may exist more than one cycles 

1 2 

connecting two given CCPs. The cycle with l.lie minimum enclosed 

area is the cycle on which the neighbourhood relation 

definition is based. Let the cycle C -C -C -C -C -C -C -C be the 

12345678 

minitiium area cycle for the CCPs CP and CP in the Fig. 3. 1.2. 2. 

12 

Let the arnas in the Fig. 3. 1.2. 2 be denoted by, 

A = area of the cycle C -C -C -C -C -C -C -C -C ; 

c 123456 781 

A = area of the polygon C -C -C ; 

1 12 8 

A = area of the polygon C -C -C ; 

2 3 4 5 

A = area of the polygon C -C -C -C -C -C -C ; 

3 1235678 
A = area of the polygon C -C -C -C -C - C ; 

4 2 3 4 5 6 7 

A = area of CP ; A = area of CP . 

5 1 6 2 
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Fig, 3. 1.2. 2 Cycle Defining Neighbourhood Relations 


3-1.3 Type one neighbourhood relations ^ 

A neighbourhood relation between two CCPs , CP and CP of Fig. 

1 2 

3. 1.2. 2, can be defined from the perspective of CP , to indicate 

1 

how important is CP as a neighbour to CP . It is called ’Type One 

2 1 

Neighbourhood Relation’ between CP and CP , hnneeforth 

1 2 

abbreviated as N1 relation between CP and CP , represented 

1 2 

symbolically as N (CP ,CP ). 

' 112 
N (CP ,CP ) is defined as 
112 

N (CP ,CP ) = (A / A ) * (A /A ); 

112 1 3 1 5 

acting along a line of action joining the. centers of the edges of 

CCPs, C -C of CP and C -C of CP .It is shown in Fig. 3 .1,2 .2, 
2 : 8 1 3 5 3 ^ 

as the tine L. 
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It can be seen that N is not commutative, i.e. H (CP , CP ) = 

^ ^ ) = (A / A ) * (A / A ); acting along the 

1 2 1 1 2 1 2 4 2 6 
line of action L. 

Typs tvo neighbourhood relations j_ 

A neighbourhood relation between the CCPs CP and CP of the Fig. 

1 2 

3. 1.2. 2 can be defined so as to indicate a global perspeci. ive of 

the relation, rather than the perspective from one of the CCPs as 

in M . It is called the ’Type Two Neighbourhood Relation’, 
1 

henceforth abbreviated as N relation, r-epresented symbolically by 

2 

N (CP ,CP ). 

2 12 

N (CP ,CP ) is defined as 
2 1 2 

N (CP ,CP ) = (A + A ) / A ; 

2 12 12 c 

It is taken to be acting along the 1 ine of action L, which joins 

the centers of the two interacting regions of the two edges 

involved in the N relation. If the said regions are collinear, 

2 

then .it i r. the same as that of N (CP ,CP j. 

112 

N is commutative, i.e. N (CP ,CP ) = N (CP ,CP ). 

2 212221 
3.1.5 Overall neighbourhood rel at ions j_ 

The overall neighbourhood relat.ion between the two CCPs CP and 

1 

CP of the Fig. 3. 1.2.2, is a triad acting along the line of 
2 

action L i . e . , 

N(CP ,CP 3 = {N (CF ,CP 3, N (CP ,CP 3, N (CP ,CP 3)- 
1 2 112 12 12 12 
The W relation between two CCPs has a unique line of action. Two 

CCPs may share a t most one edge. It can be seen in Fig. 3 . 1 . .5 . 1 , 

that the two CCPs CP and CP must interact through the shared 

1 2 

edgtv only. The cycle C , is the minimum area cycle. Any other 

i- 

cycle, such as C , would have a larger area than cycle C , tho 

2 1 
cycle through the shared edge. 
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An edge of a CCP , which it shares with another CCP , may be 

contained in the edge of the second CCP. In Fig. 3.1.5.?., the edge 

C -C of CP is contained in the edge C -C of CP . The edge C -C 
12 2 131 13 
of CP thus is shared with more than one CCP or with exterior. The 
1 

edge C -C is divided by introducing a vertex with 180 degrees 
1 3 

interior angle .it C . It gives rise to a zero inteiu-i i l.y skeletal 

2 

line, shown dashed in the figure. Thus the N relations of CP with 

1 

CP and the other CCPs or the exterior interacting through edge 
2 

C -C , are appropriately determined. 

1 3 

Two CCPs in an N relation, may share a vertex, rather than an 
edge. In such ,i case, the line of action is defined as the line 
passing through the common vertex, with inclination equal to the 
average of the inclinations of the skeletal lines of the two CCPs 
emanating from the vertex. 

The N relation may be; extended to cover the interact ieni between a 

CCP and the exterior of the profile. In .absence of a bay analysis, 

the exterior may be treated as a polygon with very large area. The 

intoraction of the CCP CP and the exterior is shown in Fig. 

3.i.5.3. Only a definition of the N relation is meaningful for C 

1 P 

and the exter'it)!-. As the skeleton of the exterior is not defined, 

definition of N is not feasible. 

2 

In Fig. 3. 1.5. 3, five different N relations can be defined for CP 

1 

and exterior. They are betuoon edge e - exterior, e - exterior, 

1 2 

o - exterior, e - exterior and e - exterior. 

3 4 ' . 5 

N C® , exterior) = Area(C -C --C ~C ) / Area (CP); along L ; 

1 1 1 2 3 4 1 

N (e .exterior) = (AreaCC -C ~C -C ) / Aroa(C ~C -C -C -C )) 

1 2 3 4 5 6. 3 4 5 6 7 

* (area (C -C -C -C ) / Area(CP)); along L . 

: . 3 ' 4. ■ 5 6 

N relations of e , e and e with the exterior are similar top 

• 1 - ■ 3 ' 4 ■ '5 ■ ■ 
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that of e . 

1 

?. i , ? , Projections of the profile ^ 

All the CCPs of a profile may not be equally important for 
nesting. The CCPs which project, out of the profile are the 
.iiiiporLant CCPs for noKfing. Such CCPs are called the projections 


of the profile. 


3 .2.1 Pro j ect ion identification j_ 

The projection identi f i cat. i on is based on the N relations of a CCf 
in a profile. A Projection would typically' liave a large peripheral 
track exposed to the exterior of the profile. In Fig. 3.2.1, the 
CCP CP is a projection of the profile. It has multiple K 

1 

relations wi !:h the exterior, K , i - 1 , . . . ,n; where n is the 

li 

rmiiiber K relations between CP and exl.<5rior. 


1 

For a projection, 'y N is close to one. Projections are 

i li 

identified on the applical.ion of the corresponding fusssy 

qualitative predicate to all the CCPs in hhe profile. 

3.2.2 Interrelations of projections 

For further analysis, only the projections are considered rather 

than all the CCPs. The projections, linked together by the N 

relations form a geometrical graph. The grapli and its constituent 
ti’tios and strings, are henceforth referred to nr. projection graph, 
projection trees and projection strings respectively. 

The important N relation for the identification of strings, irees 
and graphs is the N relation. An N relation of value almost 

. . 2 : 2 ■■ f 

equal l:o one indi cartes -a-N relation linking two projections in to 
a string. ^ 
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Various projection N .interrelations ar-o illustrated in Fig. 

3. 2. 2.1. In part (a) the projections P , P and P form a 

12 3 

projection string. In part (b) P , P , P and P form a tree. In 

12 3 4 

part Cc) P , P , P and P for'm a cycle i.e. a graph . 

12 3 4 

The distinction between projection strings and tr-oes may not be as 

clear as in Fig. 3. 2. 2.1. The Fig. 3. 2. 2. 2 shows an ambiguous 

case. Tho projection P is only linked to P . But this is 

2 1 

definitely not a string but a tree. Such cases are identifiable by 

the presences of two bays B and B', between P and P . Th' >/ may be 

1 2 

treated as trees by dividing l.he projection P in to two at one of 

1 

the concave vertiems of the profile on P . 

1 

3.3 Integration of bay analysis and the profile : 

The integration of bays in to the profile model is done after 
projection .-inalysis. The bays are analysed up to and im-lusive of 
projection relations. Since the bays are also treated as profiles, 
they may have projections, called the bay-projections. The.se bay- 
projections are in contact with the pro jectionn of the profile. 
The first step in Integration of bays is the formation of N 
relations of the bays and projections. The N relations forge links 
between various CCPs of profile and its bays. Again a geometric 
grapti, called the Bay - Profile N relation graph, is formed by the 
K relations. The process of forming the N relations between the 
bay-projections and the projections (of the prof.il c), is the same 
as that of the N relation formation between l.ho CCPs of the 
profile as given earlier. 

The model of the profile now consists of th«; projections, bays, 
and the bay - profile N relation graph. For a complex profile, 
with a large number of intermingled l»ays and projections, the 
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grrijj)) ill extremely complex. A thorough analysis of such a is 

beyond I. ho scope of the present work. Some of the oompl exiti es of 
a pT’ofile may have been removed by th<; approximations in the 
profile; decomposition. In spite of it, the bay - profile N 

relation graph may remain very complex. 

The model of the ptofile, at this stage is that of a goometrical 
graph i.e. a graph in which the nodes and the links; ai-e 
geometrical entities which bear- srhrict positional relations in 
addition t.o the topological relations embodi ec3 by the graph. The 
gooinetrical relations are elaborated in the; next section. A shape 
matching based on the dcnoi- iption of a profile at this stage 
involves an inexact partial matching of a graph, and then the 
matching of the geomotrical aspects of the rela-tiono. Inexact 
matching of a graph has been explored by Shapiro and Haralicfc 
[33]. The concepts developed by them could be explored and 

extended for a .shape matching algorithm. It is not attempted in 
the present work. 

The complexities are sought to be reduced by approximations. The 

secondary features such as the bay-projections and the low 

intensity projections on the projections are deleted by merging 

the bays between the secondary projections. Thus the bay- 

proj<;ctions are approximated out of the profile model. The profile 

model at this stage, consists of simple strings of bays and 

projections. Subgraphs, if they exist, are among the members of a 

string. An example of such a string relation and a subgraph is 

given in Fig. 3.3.1. In part (a), the projections P , P and P 

12 3 

form a string with bays B and B . In part (b), a subgraph is 

1 2 
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formed among the bay B and the pro jectionr, P and P . Thus, in 

i 1 • ■ 2 

the further work, the bayu and projections will be deemed to be 

only first ordnr features, and an absence of any second order 

featurc^ii wFli be assuraed. 

Ar. l>ie shape matching for nesting is likely to be inexact and as 
it is to a containment matching, i.e. the profile is to be 
contained in the complimentary profile, some of the bays may not 
be utilised at all. To facilitate a match urider such 

circumstances, it is desirable to merge bays in to the profile in 
combinations and analyse the shape of the resulting profile (ref. 
section. 2.3). A multitude of approximattul models of the profile 
are formed by merging the bays in to profile in all possible 
combinations. The profile is now reprnu-.ented as a number of 
approximated models. 

3 . 4 Geometrical relations ^ 

The N relations are fixed in the plane of the figure. They act 
along their particular lines of art-.ion. The geometrical relaf.tons 
encode the geometrically located nature of the K relations and the 
relative location and the orientations of the CCPs or projections 
and bays. The geometrical nature of the N relations has already 
been elaborated earlier. The relative locations and relative 
orientations of the CCPs are quite easy to formulate. It is done 
by eo l.ablishing a local co-ordinate syste.m on each of the CCPs and 
storiTig the transformations between the local c.o or-dinate systems. 
The co-ordinate systems are polar co-ordinate system.*?. This 
selection helps in ensuring tho independence of the description 
with respect to the mirror imaging transformation. Thcs .axes of the 
co-ordinate systems are the feature axes. The important feature 
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axes are selected to form th 


e zero ar)£] e direction of the polar 

CO ordinate system. The zero angle dir(i(:l. i ons are taken to b.e 

those which poiT)t. towards the exterior of the profile and the 

convex hull of the profile. Thus multiple co-ord i tioS: i- systems are 

associated with every CCP , In Fig. 3^4.1, a simple profile is 

shown. It is composed of the CCPs C , which i 13 a rectangle, C 

1 2 
which is a triangle and C which is a square. The feal-ur'o axes for 

3 

all the thr-ci <3 are drawn dotted in the figure. Each of these 

feature axes can function as the co-ordinate axis direction. In 

tho figure, only one co-ordinate system has been drawn for- each of 

the CCPs. Ray 0 -0 ’ in the co-ordinate axis for Cl. Ray 0 -0 ’ is 
11 2 2 
thn co-ordinate axis for C . Ray 0 -0 ’ is the co-ordinate axis 

2 3 3 

for C .The transformations between the pairs of local co-ordinate 
3 

systems can be written down trivially. The lines of action of the 
N relations have also been lihown in the figure. 

The comp! ate polygon description constitutes of the K-relation 
graph, the descriptions of the CCPs, and the geometrical 
rel al. i ons . 

4 . Example of Profile Description 

The profile, given in chapter 3, section 16 (ref. Fig. 16.4, 
chapter 3) is again taken as the example for profile description. 
1) The N-relations of the CCPs are taken. They are very simple as 
there; are only three CCPs. The N-rel at ! ens of the CCPs, their 
linos of action and the cycles of the N-relations, ar-e .shown in 
Fig. 4 . 1 . As all the throe CCPs have strong N-relations with the 
exterior, all are termed as projections. 
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The profile, after identification of bays and the N-relations of 

all hhe project! onn and bays, is shown in Fig. 4.2. Thnre are 

three projections, P , P and P (corresponding to CCPs C , C and 
,12 3 12 

C respectively); and the bays B , B and B are interspersed 
3 12 3 

among tin- pi’o jections . 

The resulting bay-projection strings are : 

P -B -P , P - B r , P -B -P , B -P -B , B -P -B , B ' P -B ; 

11 2 223 331 112 223 311 

P -B -P -B , P -B -P -B , P -B -P -B , B -P -B ■ J> , B -P -B -P , 
1122 22 33 3311 1223 2331 

B -P -B -P ; 

3 112 

P -B -P -B -P -B . 

1 1 2 2 3 3 

The sti-ings P -P and other such strings are not considered as 
1 2 

they are not bay-pro j or. 1. i on strings. They are not peripheral and 
hence do not influnnce nesting allocation. 

The local co-ordinate systems for the bays and projections also 

have been shown in Fig. 4.2 in the form of arrows. The origins of 

the local co-ord i tial.t; systems are the locations of the unifork 

center, indicated by 0 in the figure. Projection P has only one 

i 1 

important feature axes. The z.r.ro angle direction of the local co- 
ord inal.e system (polar) is directed along the feature axis 0 -V , 

I 8 

and directed towards V from 0 , pointing towards the exterior of 

8 1 

the profile. P has one important feature axis 0 -0 ’ . P has 

2 2 2 3 

mull- i pie feature axes of comparable importance; F -V (reversed 

3 13 

to be directed from V to F ) is the most important, F -V , F - 

13 3 122 

V , F -V and F -V are the otVier less important feature axes. 
1 ’ 13 3 12 

The approximate symmetry axes formed by combinations of feature 

axes F -V , F -V , F -V and F -V , i.e. 0 -0 ’ is the most 
1 2 * 2 1 1 3 3 12 3 3 

important feature axis, hence the best candidate for local co- 
ordinate direction. For the. bays , the symmetry axes formed by 
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combinations of the featijin ,j.xes form important feature axes. Thus 

for Bl, Ob -Ob ’ is the local co-ordinate direction. Similarly, 

1 1 

for- T\ and B , Ob -Ob ’ and Ob -Ob ’ form the local co-ordinate 
2 3 2 2 3 3 

directions . 

The N-relationship graph, the geometrical a.nix'.cts of the N- 
relations, the -ometrical relations between the local co cjr-dinate 
systems, along with the CCP descriptions given in llie earlier 
chapters, constitute the description of one ’vci ci.on’ of the 
profile. The bays are then mcrgcti in to the profile in 
combinations and hhe same exercise is carried out. The bay merged 
profiles (BMPs) are illustrated in Fig. 4.3. 

The nliape descriptions of all the BMPs forms l.he complete 
dcncription of the profile. 




Chapter 6 


CONCUJSIONS 

1 ■ Summary j_ 

"Iho aim of the thesis has been a full design automation, rather 
tlian an implementation of an interactive pro£si;unme for nesting 
allocation. Shape analysis has been recognised to be crucial to 
any fully automatic nesting solution. Hence the basic thrust: of 
the work is on shape analysis consistent with human perception of 
shape. 

The problem of analysis of an arbitrary two dimensional profile is 
solved by decomposing the. profile in to a set of cDuvex polygons. 
The descT-iption of the constituent convex polygons, their inter- 
relations, anti l:he boundary features, such as projections and hays 
generated thereof, form the d<;scription of the shape of the 
profile. Thus the shape deuen- i pt ion problem has been uulxHvided in 
to subproblems which have been solved independently. The shape 
description generated thereby will form the basis for fuiM.her 
development of shap<! matching and nesting algorithms. 

The aubproblems and their solutions 

Polygon Decomposition : Pattern Recognition and computational 
geometry. In profile decomptasi l:ion, which is a standard 
Computational Geometry / Pattern Analysis problem, the aim of the 
decompos i t i t>Ti has been modified from a decoroptis i I; ion in to minimum 
nvniilxsT' of polygons in to a decomposition optimising a heurisH. m' 
function which reflects cons it; tuncy with hu'man perception. The 
skeleton and skeletal interaclions of the concave vertices guide 
the decowpot; i tion cuts . A decomposition algorithm, I'.tnod on the 
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skeleton, and approximations, has been proposed. It is a 
heur ist ically guided si>ace search algorithm. The heuristics have 
been developed as fuzzy meitihcrnhip functions. Hence it also fits 
in to a fuzzy paradigm. 

Convex Polygon Description : A scheme has been developed for a 
meaningful approximate dciscr iption of an arbitrary convex polygon. 
It is based again on the special fealur-os of the skeleton of the 
convex polygon, such as the convexi Ly of the radius of the 
skeletal circle as a function of the distance along th(! skeleton. 
The classification of tint convex polygons, which is the basis of 
the de-scr i p I. i on, is a fuzzy classification. 

I’T’ofile Description : The proposed profile dciscription is based on 
the gcoiiiotrical and topological inter-relations of the const itur.n I. 
convex polygons. Features of the profile useful for nesting, such 
as pro j oc Lions , bays and their strings are recognised on the basis 
of these relations. An enumeration of these features along wi l;h 
the inter-relations, and the description of the constituent convex 
polygons is the profile description. The proposed algorithm 
gene.Tal.es multiple approximate descriptions. The furl.her work on 
shape matching and nesting algorithms should be based on the 
profile description. 

Shape matching and nest it ig algorithms have not been developed due. 
to the time consLraint on the thesis work. A thorougli study of 
nesting algorithms has not been done; rather, the shape analysis 
approach has bottn pursued fully. The emphasis in the present work 
has been on the (levelopment of a coHiplel.e conceptual framework for 
shape analysis and concomit-n> I- shape matching- The proposed 
framework can completely handle any non-self intersecting polygon 
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wa ll-i competence. The algorithms in the work have not been worked 
down 1.0 last details at places; but the problems have been reduced 

to Jitandard problems in literature, for wli 5 ch complete solutions 
are known. 

The present work is seen to be a fundamental work on two 
dimensional tseometry. The approach developed in this work can be 
applied to diverse two dimensional geometry problems. As an 
illustration, an auxiliary probltnn in manufacturing automation, 
ban been solved. The conceptual algori l.iim giving a complete 
solution has been given in the appendix. 

2 . Scope for Further Uork in the Probl em j 

Since the present work is a conceptual framework development 
rather than a development of an algorithm, some of tiie standard 
problems to which l.iio solution has been reduced, need to be 
stuciied further, for the purpose of an selection of the exact 
algoril.hm for the solution of the standard problem. 

For the fuzzy- logical operators and fuzzy inferencing, it is 
suggested that the generalised mixlus ponens developed by Didier 
and Prade [14 J, be followed. In addition, the fuzzy logical 
op<iT-ators may be fine tuned by experimentation. 

The heuri.sl.ic functions / the fuzzy qualitative predicates 
proposed in the present work also should be improved and fine 
tuned by experimentation over a wide range of shapes. 
Par1. i cn.larly , the qualitative predicates for the proximity and 
size of the fork;;, in a convex polygon, need extra atlention. 

In profile decomposition, apart from cuts based on the skeletal 
interaction, the cuts generated by extending the edges meeting at 
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the concave vertex should bo explored. 

In the profile dnscription, the projection-bay graphs have been 
aii.ilysed only after itiergiing the bays retjuir-ed to reduce the 

complex c-r-.-iph to a set of strings. For a complete filiapo analysis, 
a tlioTough explorative study of the graphs may be done for a 
possible classification, as don<! in convex polygon description. 

The shape description model is strictly additive in nal.uio; i.e. 
the constituent convex polygons are only added to form the final 
profile. The concept of developing an additive as well as a 

subtra<Ctive model, wherein the enveloping approximations madn for 
the sake of .“i i mplicity , may be considered to bo ( he approximation 
minus the differential area bel.wnen the approximation and the 
original chape, seems to bo a good avenue for further work. 

The whole shape descripl. ion has been divided in to mutually 
exclusive set of subproblems. It has been assumed that the 

i-iubprobl ems are mutually exclusive, but in reality it is not so. 

The subproblems do affect ono another; e.g. In giving the CCP 

iloncription, the exposure to the exterior should bo considered 
while deciding the dominating features of the CCP. These cross 
effects should be studied, and explored further. For shape 
matching algor i l.hms , the concepts developed by Shapiro and 
llaralick [33] may be found useful. Before proposing a nesting 
algorithm, a thorough slucJy of nesting stral.ogies, vis a vis the 
range of shapes, sisses, number of profiles, and the d i s l.r ibution 
of shape and irixe on the number of profiles must be done. For 

n< -il ing based on shape matching, clustering is an important 

nestii:ig strategy. The potential of the shape matching approach 

would be more apparent in an implementation which allows 
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clustering of the profiles in to simpler shapes and allocation of 
the clusters. 

The posioiblo cross-effects of the shapes of other profiles have 
not been con-gidered in the present work; e.g. If a profile has a 
feature such as a bay, which i 0 (luite small and there are no 
projections which may poowibly fit this bay in the other profiles, 
then it would be better to merge this bay in to the profile, 
rather than considering it to be a bay and there by cmplicating 
the .'j.rifj.lysia. Also, it would be better to approximate the features 
of a small profile, as the increase in nesting efficiency would 
not amount to mufb; and the complexity of the problem would be 
reduced. Such cross-effects, if considered, would greatly enhance 
the usefulness of any nesting imp! tuiientation . 

3 . Generalisation of the work j_ ^ 

Polygons with holes can be hantil cd through the clustering 
approach, by allocating small profiles in the holes; otherwise, 
holes do not play any pari, in nesting. After allocating small 
profiles in the holes, ti»c rmnaining area should be merged in to 
th<! profile ai'id shape analysis should be conl.inued for the merged 
profile. 

As menti(MM'.d earlier, curved boundary handling in the present work 
is not completely satisfactory. For a fully complete development 
capable of handling any two dimensional shape, curves should be 
handled as curves and ttofc as linear approximations. The pi<;cont 
work eij-n not be directly extended to profiler with curved 
boundaries as l;he profile decomposition depends on concave vertex 
elimiTval.ion; and also, in convex polygon description, the 'number 
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of aides' play an important ro.l e in the analysis- An investigation 
in to the curvaturci i>1on^ the boundary of the profile mi£ 5 ht: load 
to a inolh(.u;l to extend the present approach to curved profiles, as 
the <'h.u»£o of sign of curvature and any diacont i tuj 1 ly in curvature 
i rul i ca. 1. the presence of a concavity. 

Extension of the shape analysis ho 3~D solids is extremely 
diffictail- It is suggested that the skeleton could again play a 
vital role in 3-D analysis, as it does in 2-D analysis. A 3-D 
shape analysis using skeleton of 3-D solids has been mooLcd in 
[26]. ri could be used as a starting poini: for the study of 3-D 
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APPENDIX 


problem solved by shape analysis approach 

1 . Statement of Problem j_ 

The generation of the machining scqtujnce for a prismatic solid, of 
polygonal cross-section to be machined on a mining center, by 
end, peripheral or face milling, or by trepanning operations, from 
a bar of some standard cross-sections such as rectangular, 
hexagonal or circular. 

2 . Discussion on the approach to the problem 2 . 

This problem can be considered to be a two d i mensional geometrical 
problem. The considerations of feasibility of the machining 
sequence from clamping of the woi-k piece, the cutting forcen 
minimisation and the optimality of the machining sequence for 
minimum time etc. can be imposed later, on the solution to the 
geometric problem. The geometrical problem must be completely 
solved first indicating th<s areas to be removed, the tools 
required to remove the arean, and the sequence of the removal of 
tho area. 

The basic approach to the problem will follow the following 
sequence. 

1) Generate a minimum area envelope of the cross- see Lion of a 
standard shape. 

2) Identify the areas to be removed, by taking the differential 
area bfstween the cross-section and the envelope. Disjoint 
components of the differential area are taken as different 
irregular profiles to be removed. These areas am referred to as 
profiles henceforth, to maintain consistency with the terminology 
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used in the thesis. 


3) It may be noted that a the milling cutters considered in the 
probltm, in one pass or in multiple, overlapping or non- 
overlapping patuKui, cut and remove rectangular areas. Hence the 
profiltsn -•should be decomposed in to rectangular conn I. i tuents which 
may overlap. The const i l.u on ts may also be triangular in shape. The 
machining o£ the triangular primitives is explained later. Each of 
the constituent primitives is cut by a milling cutter in one 
machining operation. To obtain a minimum number of machining 
operations , the number of primitive const i l.uents should be 
minimum. 

4) To obtain a feanihio machining sequence, the rectangular 
constituents should be ordered on the basis of their exposure to 
the exterior, or their external visibility. The exterior referred 
here is the exterior of the polygon composed of the cross-section 
and the parl.n of the profiles yet to be removed, it is referred to 
as the intermediate. cross-section henceforth. 

5) The feasibility of eacli of the machining operations should be 
checked for the accessibility of the primitive at the time of 
removal, and the basic feasibility of machining such a primitive. 
For example, a pocket may not be accessible as the tool holder may 
interfere with the work-piece, a triangular notch can not be 
nincbined by a milling operation. A machining sequence which has 
been checked for the feasibility of all the machining operations 
in it, is referred to as a verified machining sequence. 

6) A number of machining sequences may be possible. Some of which 
may turn ou t to b e inf eas ibl e f rora cl amp i ng and other 
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considerations. A completely automated proceosi planning should 
then consider the force analysis also. That is out of the purview 
of the present work which is strictly geometrical in nature. The 
geometrically feasible verified machining sequences may be all 
displayed l.o l.he user and the user may be aske<} to select one the 
machining sequence which is deemed to be the best for the t:l;unping 
and other considerations. 

Each of these sLeps is elaborated upon in the following sections. 

3 . Standard Envelope Generation 

This problem is not a very difficult one. A convex hull of the 
cross-section can be taken and the envelope may be built around 
it, on the basis of the available standard cross-sections and 
their sizes. Identification of the disjoint differential areas 
between the cross-section and the envelope is again a simple 
problem. 

4 . Decomposition j. 

The profile (the profile is the area to be machined away) may be 
decomposed in to convex polygons. Each of the convex polygons may 
then be decomposed in to rectangles and triangles quite trivially. 
This decomposition may not be optimal but it can form the basis of 
a .solution methodology which can completely solve the problem and 
exactly identify and indicate the polygonal cross-sections which 
can not be machined by milling. The feasibility of the machining 
s<!<iu«>nces which may be obtained from the decompositions, are 
discussed in the next section. A decomposition which can give a 
foaiilble machining sequence is called a feasible decomposition. 
Multii>le feasible decompositions may be for-med, and an optimal 
decoraposit ion may be selected from the set of the feasible 
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decomposit iorjK . The heuristics for guiding tht; decomposition 
procnas should be based on the cutting cona i (Jorations . 

^ Feasibility of a Machining Sequence j. 

The feasibility of a machining sequence is the feasibility of the 
constituent machining operations. A machining operation in fully 
identified when the milling cutter required to remove the 
primitive under consideration is identified. At a stage in the 
determination of the machining sequence, a primitives may be in the 
situations shown in the Fig. 1. 

A rcictangular primitive, with atlrsast one side fully exposed to 
the exterior of the intermediate cross-section (ref. step 4 of 
section 2) may be milled by approaching the rectangle from that 
sides, as shown in the part (a) of figure 1. Tenl.alively select a 
milling cutter, from those available, for th<s rectangle. Determine 
the area swept by the cutter as it approaches the rectangle. If 
the swept area does not intersect with the Intermediate cross- 
section, then the machining operatiots is feasible. 

A recl.angle with only a part of a side exposed to the. exterior of 

intermediate crosn flection, necessarily ii<s<sds a milling operation 

followed by a trepanning operation. The feasibility of this 
operation again depends upon the int erf orence of the swept area of 
the cutter with the intermediate cross-section and also upon the 
position of the exposed boundary segment (Fig. ICb)). 

A triangle which has two sides exposed, may be matthined as shown 
in part (c) of the figure 1. A triangle with only one side exposed 
to tlif! (ixterior, can not be milled. The rest of the feasibility 
considerat i f»ns are the same as the ones given earlier. 
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Fig. 1. Area Primitives and Hilling Operation 
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6 . Machining Sequence j_ 

The n e i £l'ibourhood relations of the priroitiven form a graph 
Csimilar to the CCP—N relation graph given in section 3 of the 
chapter 5) The neighbourhood relations may be the adjacency 
relations i.e. the existende of a shared edge between two 
primitives, ins’tea<] of the N-relations developed in the thesis. 
Determination of a machining sequence is the determ? nation of a 
hierarchy in the graph, in teriiu; of the necessity of the prior 
removal of a primitive with respect to the other primitives. The 
directed graph formed by the imposition of the hierarchy, and 
after an exhaustive machining operation feasibility check 
(verification), given the set of feasible machining sequeneen. The 
verified and directed primitive neighbourhood relation graph is 
henceforth abbreviated an VDPN graph. 

The primitives which are directly exposed to the exterior of the 
envelope may be removed first and hence they form the roots of the 
graph. In Fig. 2, primitives A and B are directly exposed to the 
exterior of the envelope. They foT«i the roots of the graph. Any of 
them can be removed first. If A is removed, then B is exposed to 
the cvxl fjrior of the intermediate cross-sec; I. i on. Note that C is not 
exposed l.o the exterior for a feasible machining operation with 
just the removal of A. After the removal of B any of C and D can 
be roiiioved in any order. If C is removed first, then D is removed, 
and finally B is removed by a trepanning operation. F.imilarly, if 
B had been removed in the first step rather than A, the graph 
could again have been follow<^d similarly. 
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At each of these steps, the feasibility of a inachinin£ operation 
must be checked before deciding upon the set of the possible 
candidates for removal. A VDJ’N graph must contain all the 
primitive nodes present in the primitive N-rolation graph; 
otherwise it indicates the infeasibility of the production of the 
object by a set of only milling operations, for the given set of 
cutters. It may be noted that a sot of cutters is associated with 
every primitive, ond thus every requisite machining operation is 
identified. The exact assignment of a unique cutter to every 
machining operation may be done later, on the basis of the 
optimality of the machining operation for simultarHsi ty of 
machining operations and minimisation of the overall machining 
time. The VDPN graph for the object in figure 2(a), obtained after 
the hierarchy imposition and feasibility check is illustrated in 
the part (b) of the Fig. 2. 

The paths from the roots to the leaf nodes of the VDPN graph, give 
the priority constrain t.q for machining operations. A set of paths 
which t:ompletely spans the primitive nodes gives a net of feasible 
machining nequences. The set of all the sets of such spanning 
paths, give tlie complete set of feasible machining operation 
sequences. The machining sequences ABCDE, ABDEC, ABDCE, BACDE, 
BADEC, BACDE are the feasible machining sequences for the object 
in f igure 2 - 

The set of the feasible machining sequences may be processed 
further for the determination of an optimal machining sequence. 
The possibility of simultaneous machining operations is indicated 
by the subpaths, within tb« net of paths giving a machining 


229 



sequence, which do not intersect i.e. which do not have any cojiitnon 
primitive nodes. 

7 . Conclusions 

The basic concepts developed in the thesis, of decomposin£ (.he 
profile in to a set of constituent primitives and forming a 
neighbourhood relation graph to structure the shape of the 
profile, have been successfully applied in a diverse geometrical 
problem of automatic process planning. Thus tlus universality of 
the methodology has been demonstrated. It may be concluded that 
the method developed in the thesis has a very good potential for 
the solution of a lot of diverse geometrical problems encounfored 
in ninr.i»anical engineering. 
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